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VORTEX-INDUCED  TRANSVERSE  LOADING  ON  AN  ARTICULATED 

TOWER 

By 

AHMED  FAHMY  OMAR 
August  1992 

Chairman:  Prof.  D.  Max  Sheppard 

Major  Department:  Coastal  and  Oceanographic  Engineering 
This  research  is  an  experimental  investigation  of  vortex-induced  transverse  loading 
on  fixed  and  compliant  structures  in  nonplanar  oscillatory  flow.  The  effects  of  flow 
nonuniformity  and  transverse  motion  on  the  transverse  force,  vortex-shedding  fre- 
quency, and  vortex  span-wise  coherence  length  have  been  investigated.  To  quantify 
such  effects,  two  types  of  experiments  were  performed.  Each  type  of  experiment  con- 
sisted of  a series  of  tests.  All  tests  were  performed  by  oscillating  an  instrumented  10  ft 
long,  6 inch  diameter  aluminum,  articulated  cylinder  in  a 30  ft  diameter  cylindrical 
tank  with  a water  depth  of  8.5  ft.  In  the  first  type  of  experiment,  transverse  motion 
was  not  allowed,  while  in  the  second  type,  the  cylinder  was  allowed  to  move  freely 
in  the  transverse  direction.  Reynolds  numbers,  Re,  between  0.61  x 104  and  1.3  x 105 
and  Keulegan-Carpenter  numbers,  KC,  between  2.4  and  9.35  were  obtained.  A new 
method,  involving  the  use  of  miniature,  quick  response  thermistors  was  employed  for 
measuring  frequency  of  vortex-shedding. 

For  the  range  of  Re  and  KC  tested,  transverse  force  amplitudes  as  large  as  70% 
of  the  in-line  force  were  measured.  Transverse  forces  occurred  at  frequencies  that 
were  multiples  of  the  driving  frequency.  The  dominant  vortex-shedding  frequencies 
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clustered  around  one  of  the  harmonics  of  the  driving  frequency.  The  specific  harmonic 
depended  on  the  value  of  KC.  A mathematical  model  for  computing  transverse  forces 
taking  into  account  the  dependency  of  the  lift  coefficient  on  both  Re  and  KC  as 
well  as  the  fact  that  these  forces  have  multiple  frequency  components  has  been  also 
proposed.  Using  this  model,  the  magnitude  of  the  transverse  force  was  found  to  be 
deterministic  and  repeatable.  However,  the  phase  of  this  force  was  random  in  nature. 
The  constrained  tower  results  for  in-line  and  transverse  force  coefficients  were  found 
to  agree  well  with  data  obtained  by  others  for  waves  acting  on  vertical  cylinders. 

Transverse  motion  was  found  to  have  a significant  effect  on  vortex-induced  load- 
ing. The  effect  of  transverse  motion  was  also  found  to  increase  the  vortex  correlation 
length,  increase  the  lift  coefficient  (by  at  least  two  and  one-half  times)  and  alter  the 
nature  of  the  vortex-shedding.  For  the  range  of  KC  tested,  the  results  showed  that, 
the  larger  the  value  of  KC,  the  larger  the  amplitude  of  maximum  transverse  motion 
and  the  larger  the  value  of  reduced  velocity,  Vr,  at  which  this  maximum  occurs. 


xvm 


CHAPTER  1 
INTRODUCTION 


1.1  Problem  Statement 

With  the  depletion  of  near  shore  oil  reservoirs,  exploration  and  production  is 
moving  to  deeper  water  and  more  remote  locations.  In  most  cases,  this  means  more 
severe  environmental  conditions  and  more  stringent  design  and  safety  regulations.  At 
the  same  time,  there  is  a great  need  to  reduce  the  cost  of  producing  hydrocarbons. 
One  way  to  cut  the  cost  is  to  optimize  structural  designs.  An  essential  component  of 
optimization  is  being  able  to  predict  accurately  the  loadings  to  be  experienced  by  the 
structure.  This  trend  has  recently  led  to  a massive  research  effort  into  the  design  and 
assessment  of  the  short  and  long  term  reliability  of  offshore  structures.  The  design 
usually  involves  three  major  steps,  first,  the  long  term  prediction  of  environmental 
conditions,  second,  the  estimation  of  the  forces  associated  with  these  environmental 
conditions  and  third,  the  determination  of  the  effects  of  these  forces  on  an  intended 
structure  and  its  ability  to  survive  the  expected  extreme  environment.  It  is  interesting 
to  note  that  by  1970,  over  10%  out  of  200  drilling  rigs  had  collapsed  and  a further 
20%  had  suffered  severe  fatigue  failure  of  structural  members  due  to  flow-induced 
vibrations  (King,  1974).  This  is  in  addition  to  frequent  shut-down  of  operation  for 
days  due  to  the  large  oscillations  induced  by  fluid-induced  forces;  see  Lewis  et  al. 
(1991)  and  Koch  et  al.  (1991). 

In  the  past  two  decades,  a significant  amount  of  research  on  fluid-structure  interac- 
tion problems  has  been  conducted.  This  research  can  be  divided  into  two  categories, 
those  concerned  with  i)  the  fluid-induced  forces  and  ii)  the  fluid-induced  response 
and/or  vibration.  The  latter  category  has  received  the  most  attention. 
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The  fluid-induced  forces  are  comprised  of  in-line  forces  and  vortex-induced  trans- 
verse (lift)  forces.  Many  studies  have  been  done  to  predict  the  fluid  in-line  forces 
acting  on  a cylinder  in  steady  and  oscillating  flow.  The  majority  of  those  studies 
have  been  based  on  the  Morison  equation  (Morison  et  ai,  1950),  where  the  two  in- 
line force  components  of  drag  and  inertia  are  identified.  The  equation  expresses  these 
forces  in  terms  of  the  velocity  and  acceleration  of  the  fluid  particles  at  the  location 
of  the  cylinder.  In  comparison  to  in-line  forces,  vortex-induced  transverse  forces  have 
received  little  attention.  This  is  most  likely  due  to  the  complexity  of  the  vortex- 
shedding  phenomena  and  the  lack  of  a clear  cut  methodology  for  analyzing  the  data. 

Structural  elements  are  constantly  subjected  to  loading  due  to  wind  and/or  ocean 
currents  and  waves.  Most  flow  situations  encountered  in  nature  are  turbulent,  non- 
planar  (nonuniform),  and  unsteady.  To  further  complicate  matters,  the  structural 
element  of  interest  is  often  in  close  proximity  to  other  members,  compliant,  and 
perhaps  partially  covered  with  biofouling.  To  predict  the  vortex-induced  transverse 
loading  and  response  of  such  a structural  element  under  such  complex  flow  conditions 
is  truly  a challenging  task,  but  nevertheless,  one  that  is  routinely  faced  by  design 
engineers  in  the  offshore  industry. 

The  problems  of  in-line  and  transverse  loading  and  the  corresponding  structural 
response  are  all  related  and  interdependent.  However,  the  quantities  involved  are 
difficult  to  isolate  and  to  measure;  and,  therefore  some  of  the  reported  data  have  been 
affected  by  the  techniques,  apparatus  and/or  instrumentation  used.  In  addition,  the 
diversity  of  the  investigators’  backgrounds,  ranging  from  fluid  to  structural  mechanics, 
resulted  in  differences  in  approaches  taken  and  differences  in  the  manner  in  which 
results  were  presented.  This  diversity  of  backgrounds  has  its  advantages,  but  the  lack 
of  uniformity  in  data  reporting  makes  it  difficult  to  compare  and  extend  data  sets. 

The  focus  of  this  research  was  on  the  investigation  of  some  of  the  fluid-structure 
interaction  problems,  in  particular  on  quantifying  the  effects  of  flow  nonuniformity 
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and  transverse  motion  on  vortex-induced  transverse  loading  and  on  developing  a 
mathematical  model  for  predicting  transverse  forces. 

It  is  extremely  difficult  to  obtain  a theoretical  solution  for  vortex-induced  trans- 
verse forces  especially  at  high  Reynolds  numbers.  The  difficulty  is  partly  due  to 
incomplete  knowledge  of  the  flow  field  around  the  structure,  and  to  problems  as- 
sociated with  the  coupling  of  structural  motion  and  fluid  flow.  Consequently,  the 
approach  for  obtaining  the  solution  to  this  problem  has  been  the  same  as  taken  here, 
i.e.  experimental. 


1.2  Research  Objectives 
The  overall  objectives  of  this  study  were 

1.  To  design  and  conduct  experiments  to  establish  the  dependence  of  the  vortex- 
shedding  process,  vortex  span-wise  coherence  length,  magnitude  of  the  vortex- 
induced  transverse  forces  and  transverse  motion  on  the  in-line  flow  and  structure 
parameters.  The  physical  model  must  be  sufficiently  large  to  cover  meaningful 
ranges  of  Reynolds  ( Re ) and  Keulegan- Carpenter  ( KC ) numbers  and 

2.  To  develop  an  improved  mathematical  model  for  predicting  vortex-induced 
transverse  forces. 

To  achieve  these  objectives,  the  following  tasks  were  performed: 

1.  Design,  construct  and  instrument  an  articulated  tower  system.  Two  designs 
were  necessary,  one  with  and  one  without  transverse  motion. 

2.  Perform  two  types  of  experiments  by  sinusoidally  oscillating  the  articulated 
tower  in  an  otherwise  still-water  tank.  In  the  first  type  the  tower  transverse 
motion  was  constrained.  During  the  series  of  tests  performed  for  this  type,  the 
following  quantities  were  measured 


4 


• the  in-line  and  transverse  forces. 

• the  frequencies  of  vortex  shedding  along  the  physical  model. 

• the  in-line  motion. 

During  the  series  of  tests  for  the  second  type,  transverse  motion  was  allowed 
and  the  following  quantities  were  measured 

• the  in-line  forces. 

• the  in-line  and  transverse  motion. 

• the  frequencies  of  vortex  shedding  along  the  physical  model. 

3.  Develop  a mathematical  model  for  predicting  the  vortex-induced  transverse 
forces. 

4.  Quantify  the  effects  of  flow  nonuniformity  and  structure’s  transverse  motion  on 
the  lift  coefficients,  frequency  of  vortex- shedding  and  spanwise  coherence  length 
of  vortices. 


CHAPTER  2 
LITERATURE  REVIEW 


In  this  chapter,  a literature  review  and  discussion  of  topics  relevant  to  the  work 
undertaken  are  presented.  The  review  is  limited  to  the  vortex-shedding  process  and 
its  induced  transverse  forces  and  concentrates  on  smooth,  rigid  cylinders  in  1)  steady, 
planar,  2)  oscillatory,  planar,  3)  oscillatory,  nonplanar,  and  4)  wave  flows.  The  effects 
of  transverse  motion  on  the  vortex-shedding  process,  span-wise  coherence  length  and 
transverse  loading  are  also  discussed  for  those  cases  where  sufficient  data  exist.  It  is 
the  author’s  opinion  that  the  above  classification  of  quantities  affecting  vortex-induced 
loading  and  response  will  be  beneficial  in  analyzing  data  from  various  investigators 
and  in  understanding  the  processes  involved. 

Analytical  and  computational  approaches  to  the  transverse  force  problem  have 
been  hampered  by  the  complexity  of  the  processes.  Thus,  most  of  the  work  to  date 
has  been  experimental,  guided  by  dimensional  analysis  techniques.  For  this  reason, 
only  experimental  work  will  be  discussed  in  this  review. 

2.1  Overview 

When  a viscous  fluid  such  as  water  or  air  flows  past  a bluff  body  with  sufficient 
velocity,  flow  separation  occurs  and  a wake  region  is  formed.  Over  a wide  range  of  flow 
and  structure  parameters  of  interest,  vortices  are  observed  to  form  near  the  points 
of  flow  separation.  For  symmetric  structure  shapes,  void  of  sharp  edges,  such  as 
right  circular  cylinders,  vortices  are  formed  on  both  sides  of  the  body.  Under  certain 
conditions  these  vortices  remain  attached  to  the  body  while  under  other  conditions 
they  are  shed  from  the  body  in  or  out  of  phase  with  each  other.  The  net  effect  of  this 
phenomenon  is  a fluctuation  in  the  points  of  flow  separation,  which  in  turn  causes 
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a time  varying  distribution  of  normal  and  tangential  stresses  over  the  body.  This 
results  in  time  dependent  in-line  and  transverse  loads  on  the  structure,  even  when 
the  flow  is  steady  and  planar  (uniform). 

The  processes  associated  with  flow  separation  are  complex  and  difficult  to  predict. 
Yet  minor  changes  in  the  separation  point  can  result  in  relatively  large  changes  in  both 
the  in-line  and  the  transverse  forces  on  the  structure.  This  flow  instability  problem 
is  sensitive  to  perturbations  such  as  those  introduced  by  surface  roughness,  motion 
of  the  body,  free  stream  turbulence,  flow  orientation  relative  to  the  structure,  flow 
around  the  ends  of  the  structure,  etc.  In  an  attempt  to  understand  and  model  this 
phenomenon,  researchers  have  isolated  various  aspects  of  the  problem  starting  with 
the  (seemingly)  simplest  case  of  uniform,  steady  flow  and  moving  toward  the  more 
complex  flow  and  structure  situations.  The  processes  are  of  course  nonlinear  and  thus 
their  individual  effects  cannot  be  simply  superimposed  to  obtain  the  combined  effect. 
However,  much  can  be  learned  about  the  mechanisms  involved  and  some  guidance  for 
the  design  engineer  can  be  achieved  by  such  a process. 

Vortex-induced  loads  are  of  interest  in  a number  of  engineering  disciplines  and 
of  particular  importance  in  the  design  of  offshore  structures.  Structural  elements  of 
interest  are  often  compliant  and  subjected  to  complex  flows  (turbulent,  oscillatory  and 
nonuniform).  For  such  a situation,  the  main  danger  from  the  vortex-induced  loading 
arises  from  the  possibility  of  resonance  created  by  the  vortex-shedding  frequency  being 
close  to  the  system  natural  frequency  (or  one  of  its  multiples).  Large  and  damaging 
amplitudes  of  oscillation  can  result  from  the  complicated  and  pernicious  mechanism 
of  resonance  which  can  occur  over  a considerable  range  of  conditions.  The  interaction 
between  the  flow  and  the  structure’s  motion  causes  the  frequency  of  vortex-shedding 
to  be  controlled  by  the  response.  This  can  result  in  what  is  known  as  “lock-in”. 
Although  the  fluctuating  pressure  that  causes  the  transverse  force  is  predominantly 
at  right  angles  to  the  direction  of  the  approaching  flow,  it  can  also  produce  dynamic 
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forces  in-line  with  the  flow.  The  induced  transverse  motion  can  also  increase  the 
time- average  drag  significantly. 

The  problems  of  in-line  and  transverse  loading  and  the  corresponding  structural 
response  are  all  related  and  interdependent  as  pointed  out  in  several  excellent  reviews 
on  this  subject  (e.g.  Bearman,  1984;  Sarpkaya,  1979;  Chen,  1987;  Griffin,  1984).  On 
the  other  hand,  transverse  forces  depend  on  the  nature  of  the  flow,  the  structure’s 
geometric  parameters,  and  (when  it  exists)  the  structure  motion.  This  is  illustrated 
in  the  diagram  in  Fig.  2.1.  As  stated  above,  the  flows  are  very  complex.  Perhaps  the 
most  surprising  thing  about  this  phenomenon  is  that  it  displays  some  degree  of  order 
and  repeatability  and  thus  predictability. 


Figure  2.1:  Flow  chart  for  cylinder-loading  response. 


The  quantities  involved  are  difficult  to  isolate  or  to  measure  and  therefore  some 
of  the  reported  data  have  been  affected  by  the  techniques,  instrumentation  and/or 
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apparatus  used.  In  spite  of  these  difficulties,  researchers  have  managed  to  acquire 
at  least  a qualitative  understanding  of  most  of  the  processes  and  how  the  various 
geometric  and  flow  parameters  affect  the  transverse  force.  The  matrix  presented  in 
Table  2.1  is  an  attempt  to  classify  the  existing  work  from  the  point  of  view  of  the 
more  important  geometric  and  flow  parameters  affecting  the  vortex-shedding  process 
and  the  resulting  transverse  forces.  Reference  numbers  for  some  of  the  more  recent 
and  historically  important  papers  on  the  various  subjects  are  given  in  the  matrix. 
For  a more  complete  list  of  references  on  these  subjects  the  reader  is  referred  to  the 
review  articles  by  Bearman  (1984),  King  (1977)  and  Sarpkaya  (1979).  As  one  can 
see  from  Table  2.1,  some  areas  have  received  more  attention  than  others  (e.g.  steady 
and  oscillatory  planar  flows  have  attracted  most  of  the  interest  while  other  flows  like 
oscillatory  nonplanar  and  waves  have  received  less  attention).  The  author  is  aware 
that  important  proprietary  research  has  also  been  conducted  in  this  field,  but  these 
results  are,  of  course,  not  available  for  review. 


Table  2.1:  Categorization  of  flow  and  structural  parameters  influencing  vor- 

tex-induced transverse  loadings  and  response  (numbers  refer  to  references). 


Flow 

Structure  Parameter 
(Rigid  Cylinders) 

Smooth 

Rough 

Fixed 

Inline 

Motion 

Transverse 

Motion 

L/D 

BR 

End 

Effects 

Steady,  Planar 

1,  18,  34,  63, 
71,  76,  77,  55, 
19,  75,  67,  33, 17 

1, 18 

1,  76,  77,  55, 
19,  75,  67, 17 

63,  33 

18,  34,  63, 
71,  75,  33 

77 

77 

19 

Oscillatory, 

Planar 

60,  61,  2,  40,  46, 
71,  74,  76,  80,  43, 
50,  32,  3,  45,  42, 
69,  47,  78,  56,  70,  17 

62,  60,  61, 
71,  45,  56, 
70 

62,  60,  61,  76,  80, 
43,  50,  32,  3,  45, 
42,  69,  47,  78, 17 

2,  40,  46,  71, 
74,  42,  56,  70 

74,80 

43 

Oscillatory, 

Non-planar 

13, 11 

13, 11 

Regular  Waves 

7,  30,  65,  12,  11, 
44,  76,  9,  66 

11 

7,  30, 12, 11, 
44,  76,  9,  66 

65, 11, 
44,  76 

7,65, 
44,  76 

Irregular  Waves 

73,  23 

73.23 
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Some  of  the  earliest  studies  were  on  steady,  planar,  relative  flow,  perpendicu- 
lar to  rigid,  fixed,  circular  cylinders  (see  King,  1977;  Fleischmann  and  Sallet,  1981; 
Blevins,  1990).  The  term  relative  flow  is  used  since  in  some  cases  the  cylinder  was 
towed  through  still  water  while  in  other  cases  the  fluid  was  forced  to  move  around 
a fixed  cylinder.  This  was  followed  by  experiments  with  oscillatory,  planar,  relative 
flow  around  circular  cylinders.  More  recent  experiments  with  oscillatory,  nonplanar, 
relative  flows  have  been  performed  both  by  moving  articulated  cylinders  in  still-water 
(Chakrabarti  et  al.,  1983)  and  by  subjecting  cylinders  to  surface  waves  (see  Bearman, 
1988a;  Isaascson  et  al.,  1976,  1977;  Sawaragi  et  al.,  1977;  Chakrabarti  et  al.,  1976). 
While  the  above  work  was  proceeding,  other  aspects  of  the  problem  such  as  the  effects 
of  surface  roughness  (e.g.  Achenbach  and  Heinecke,  1981;  Bearman,  1988a;  Sarpkaya, 
1976a,  1976b,  1990b;  Wolfram  et  al.,  1989),  in-line  and/or  transverse  motion  (e.g. 
Chakrabarti  et  al.,  1984;  Maull  and  Kaye,  1988;  Bearman  and  Hall,  1987;  Donazzi 
et  al.,  1981;  King,  1974;  Laird,  1962;  McConnell  and  Park,  1982a;  Sarpkaya,  1978; 
Sumer  and  Fredspe,  1988;  Tprum  and  Anand,  1985;  Verley,  1980),  structure  aspect 
ratio  (i.e.,  cylinder  length  to  diameter  ratio)  (e.g.  West  and  Apelt,  1982),  flow  block- 
age (e.g.  Kozakiewicz  et  al.,  1991;  Tprum  and  Anand,  1985;  Ramamurthy  and  Ng, 
1973;  Yamamoto  and  Nath,  1976),  free  stream  turbulence  (e.g.  Tprum  and  Anand, 
1985),  end  effects  (e.g.  Tprum  and  Anand,  1985;  Matten  et  al.,  1978),  etc.  were  also 
being  investigated.  One  of  the  factors  affecting  the  overall  lift  coefficient  is  the  span- 
wise  coherence  length  of  the  vortices  (i.e.  the  length  along  the  axis  of  the  cylinder 
where  the  vortices  are  being  shed  in  unison).  Even  though  this  parameter  length  has 
been  the  subject  of  much  discussion,  very  few  measurements  with  circular  cylinders 
in  water  have  been  reported  in  the  literature  (see  King,  1977;  Wolfram  et  al.,  1989; 
Obasaju  et  al.,  1988;  Kozakiewicz  et  al.,  1991).  It  is  thought  that  this  length  de- 
pends on:  cylinder  aspect  ratio,  end  effects,  free  stream  turbulence,  nonuniformity  of 
the  in-line  flow,  roughness  and  roughness  gradients  along  the  cylinder,  two  and  three 
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dimensionality  of  the  flow  (such  as  would  be  generated  by  uni-  and  multi-directional 
surface  waves,  respectively),  nonaligned  currents  and  waves,  etc. 

During  the  past  few  decades,  much  has  been  learned  about  the  vortex-induced 
transverse  force  problem;  yet  it  is  safe  to  say  that  more  work  is  needed  before  reli- 
able information  required  by  design  engineers  is  available.  In  the  following  sections  a 
review  of  the  vortex-induced  transverse  forces  (in  particular,  the  lift  coefficients  and 
frequency  of  vortex-shedding)  and  span-wise  coherence  length  of  vortices  on  smooth, 
fixed  and  transversely  unconstrained  cylinders  in  1)  steady,  planar,  2)  oscillatory,  pla- 
nar, 3)  oscillatory,  nonplanar,  and  4)  wave  flows  is  presented.  Lock-in  conditions  and 
their  relation  to  the  vortex-shedding  frequency  and  the  system  natural  frequency  are 
also  presented  where  sufficient  data  exist.  Flexible  cylinders  will  not  be  discussed  in 
this  review.  This  is  because  it  is  difficult  to  devise  experiments  with  flexible  cylin- 
ders that  provide  quantitative  results  specifically  on  the  effects  of  flow  nonuniformity. 
Results  obtained  from  the  more  easily  controlled  experiments  discussed  here  could; 
however,  be  helpful  in  analyzing  the  loading  and  response  of  flexible  cylinders. 

2.2  Steady,  Planar  Flow 

2.2.1  Transversely  Constrained  Cylinder  in  Steady,  Planar  Flow 

Vortex-induced  transverse  forces  on  a smooth,  circular  cylinders  in  steady,  planar, 
flows  have  received  considerable  attention  during  this  century  (see  e.g.  King,  1977; 
Fleischmann  and  Sallet,  1981).  For  this  type  of  flow,  it  was  found  that  Reynolds 
number,  Re  (relative  magnitudes  of  inertia  and  viscous  forces),  was  the  most  im- 
portant parameter  to  characterize  the  flow  around  a circular  cylinder.  For  steady, 
turbulent  flows  around  right  circular  cylinders,  the  in-line  loading  has  been  formu- 
lated in  terms  of  a drag  force  that  is  proportional  to  the  square  of  the  relative  speed, 
the  mass  density  of  the  fluid,  and  a projected  area.  The  constant  of  proportionality 
is  one-half  the  drag  coefficient  (CD/ 2).  Plots  of  CD  versus  Re  for  this  flow  situation 
can  be  found  in  numerous  publications  (e.g.  Sarpkaya  and  Isaacson,  1981).  In  the 
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absence  of  a better  formulation,  a similar  model  for  the  transverse  or  lift  force  in  terms 
of  the  in-line  relative  speed  squared,  fluid  mass  density,  and  projected  area  has  also 
been  used.  In  this  case,  the  constant  of  proportionality  is  one-half  the  lift  coefficient 
(Cl/ 2).  In  both  cases,  the  coefficient  is  a catch-all  term  whose  value  depends  to  some 
degree  on  the  flow  and  most,  if  not  all,  of  the  quantities  given  in  Table  2.1. 

Sarpkaya  and  Isaacson  (1981)  showed  that  there  is  considerably  more  scatter  in 
Cl  data  for  this  flow  than  say  for  Co-  In  an  attempt  to  show  both  the  consistency 
and  the  scatter  in  the  available  experimental  data,  C/,(rmj)  (=  2 Fyrma/ pwDU2)  versus 
Re(=  UD/u ) has  been  plotted  in  Fig.  2.2  for  twelve  different  investigators.  Most 
of  these  data  were  taken  from  Sarpkaya  and  Isaacson  (1981),  but  the  results  of  two 
more  recent  studies  (Dronkers  and  Massie,  1978;  van  der  Vegt  and  van  Walree,  1987) 
have  been  added.  However,  it  should  be  pointed  out  that  many  of  these  data  are  for 
experiments  conducted  some  years  ago  in  air.  In  addition,  it  is  possible  that  at  least 
some  of  the  cylinders  experienced  transverse  vibration  during  the  tests.  Dronkers  and 
Massie  (1978)  had  a fixed  vertical  cylinder  in  a circulating  flume  while  van  der  Vegt 
and  van  Walree  (1987)  towed  a horizontal  cylinder  in  still- water.  These  data  appear 
to  only  add  to  the  scatter.  This  lack  of  agreement  in  the  reported  data  has  been 
attributed  to  differences  in  free  stream  turbulence,  aspect  ratios,  flow  uniformity,  etc. 

In  structural  design,  the  frequency  components  of  the  transverse  force  are  perhaps 
as  important  as  the  magnitude.  The  basic  phenomenon  of  vortex-shedding  is  illus- 
trated in  Fig  2.3,  where  the  major  regimes  of  vortex-shedding  from  a right  circular 
cylinder  in  steady,  planar  flow  are  sketched,  based  upon  the  observations  of  various 
investigators. 

Recently,  a flow  visualization  by  van  der  Vegt  and  Walree  (1987)  has  confirmed 
some  of  the  patterns  by  which  the  vortex-shedding  process  occurs  as  presented  in 
Fig.  2.3.  In  all  the  investigations  carried  out  in  this  category,  the  frequency  with 
which  the  individual  vortices  are  shed  was  found  to  be  proportional  to  the  ratio  U/ D, 
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Figure  2.2:  Lift  coefficient  versus  Reynolds  number  for  steady,  planar  flow  around  a 
smooth,  fixed  cylinder. 

so  that  fv  = StU/D,  where  the  constant  of  proportionality  is  called  Strouhal  number 
and  fv  is  the  frequency  of  vortex-shedding.  Later  St  was  presented  versus  Re  in  an 
envelope,  within  ±10%  accuracy  over  a large  Reynolds  number  range,  see  Fig.  2.4. 
Transverse  force  spectra  (in  this  case,  spectral  density  versus  Strouhal  Number)  for 
various  values  of  Re  are  also  shown  in  Fig.  2.4.  These  spectra  were  computed  by 
Schewe  (1983),  using  data  from  steady  flow  around  a fixed,  rigid  cylinder  in  a pres- 
surized wind  tunnel.  It  is  important  to  note  the  relative  energy  levels  in  the  various 
flow  regimes.  Single  spike  spectra  similar  to  Schewe ’s  spectra  for  Re  = 1.3  x 105 
(see  Fig.  2.4)  were  recently  obtained  for  water  flows  by  van  der  Vegt  and  van  Walree 
(1987).  Thus  for  all  practical  purposes,  it  appears  that  vortices  are  shed  at  a single 
frequency  for  60  < Re  < 2 x 10s  and  well  defined  by  a Strouhal  number  of  about  0.2 
for  smooth  cylinders. 


n«  < 5 REGIME  of  unseparated  flow 


5 TO  15  < R.  < 40  A FIXED  PAIR  OF 
VORTICES  IN  WAKE 


40  < Rt  < 90  ANO  90  < Rt  < 150 
TWO  REGIMES  IN  WHICH  VORTEX 
STREET  IS  LAMINAR 


150  < Rt  < 300  TRANSITION  RANGE  TO  TURBU- 
LENCE IN  VORTEX 

300  < R»  ?:  3X  10s  VORTEX  STREET  IS  FULLY 
TURBULENT 


3 X 10s  ^ Rt  < 3i  X 10S 

LAMINAR  BOUNDARY  LAYER  HAS  UNDERGONE 
TURBULENT  TRANSITION  ANO  WAKE  IS 
NARROWER  ANO  DISORGANIZED 


3.5  X 10s  < R» 

RE-ESTABLISHMENT  OF  TURBU- 
LENT VORTEX  STREET 


Figure  2.3:  Regimes  of  steady,  planar  flow  across  a smooth,  fixed  circular  cylinder, 
(Ref.  14). 


Fatal  (1989)  and  others  presented  the  relationship  of  Fig.  2.4  in  the  form  of 
empirical  equations: 

fv  = (StU/D)  Re  < 60  (2.1) 

fv  = St  (1  + 19.7 / Re)(U/ D)  60  <Re<  2 x10s  (2.2) 

In  the  turbulent  regime  (2  x 10s  < Re  < 7 x 106),  the  Strouhal  number  varies 
between  0.15  and  0.4  depending  on  the  intensity  of  the  free  stream  turbulence.  In 
this  region,  the  spectrum  is  broad  banded,  reduced  in  magnitude  and  very  sensitive 
to  flow  disturbance.  In  the  supercritical  regime  ( Re  > 7 x 106),  the  spectrum  becomes 
narrow  banded  once  again. 
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Figure  2.4:  Strouhal- Reynolds  numbers  relationship  with  transverse  force  spectra  for 
steady,  planar  flows  around  a smooth,  fixed  cylinder  (Refs.  14,  67). 


The  vortices  shed  along  a cylinder  can  be  in  or  out  of  phase  with  each  other. 
The  total  transverse  force  on  the  cylinder  is  very  sensitive  to  these  phase  angles.  It 
is  important  to  understand  what  quantities  control  or  influence  the  shedding  pro- 
cess. A parameter  known  as  the  span-wise  coherence  length  (a  length  over  which  the 
vortices  are  considered  well  correlated)  has  been  used  as  a measure  of  the  coherence 
between  the  vortices  shed  along  a cylinder.  In  general,  the  greater  the  coherence 
length  the  larger  the  total  transverse  force.  This  correlation  length  has  been  observed 
to  vary  with  Reynolds  number,  surface  roughness  and  free  stream  turbulence  (King, 
1977).  The  only  data  found  that  demonstrate  the  parameters  influencing  the  coher- 
ence length  support  the  importance  of  the  Reynolds  number  (Re).  These  data  (by 
Scruton,  1967;  taken  from  Overvik,  1982)  are  shown  in  Fig.  2.5  and  generally  indicate 
a reduction  in  the  span-wise  coherence  length  with  increased  Reynolds  number. 
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Figure  2.5:  Span-wise  coherence  length  versus  Reynolds  number  for  steady,  planar 
flow  around  a smooth,  fixed  cylinder  (Ref.  68). 

2.2.2  Transversely  Unconstrained  Cylinder  in  Steady,  Planar  Flow 

When  the  structure  is  allowed  to  move  in  the  transverse  direction  and  its  support 
is  such  that  a restoring  force  exists  (resulting  in  a system  natural  frequency),  there  can 
be  strong  interaction  between  the  transverse  response  and  the  transverse  loading  as 
indicated  in  Fig.  2.1.  Extreme  caution,  therefore,  should  be  exercised  when  using  lift 
coefficient  data  obtained  for  fixed  cylinders  when  the  structural  element  is  compliant. 
When  the  vortex-shedding  frequency  is  in  the  vicinity  of  the  natural  frequency  of 
the  cylinder-support  system,  large  transverse  excursions  can  occur  which  in  turn 
result  in  further  increases  in  the  transverse  force.  This  large  amplitude  motion  at 
the  natural  frequency  can  change  the  frequency  of  vortex-shedding  to  that  of  the 
oscillation  frequency  and  “lock-in”  occurs. 

Even  though  a number  of  experiments  have  been  performed  with  cylinders  free 
to  move  in  the  transverse  direction,  only  a few  investigators  (e.g.  King  et  al.,  1973; 
King,  1974;  Griffin  and  Koopman,  1977;  Tprum  and  Anand,  1985)  have  reported 
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information  on  the  response  of  the  cylinder.  The  spectra  of  the  transverse  force 
and  corresponding  response  by  Tprum  and  Anand  (1985)  were  the  only  data  found 
that  show  the  relationship  between  the  vortex-shedding  frequency  (/„)  and  response 
frequency  (/r).  Tprum  and  Anand’s  investigation  was  primarily  to  study  wall  effects, 
however,  their  results  for  the  largest  cylinder  gap  to  diameter  ratio  (i.e.,  G/ D = 3) 
shown  in  Fig.  2.6  should  be  very  similar  to  an  unobstructed  cylinder.  For  these  data, 
the  authors  did  not  report  the  values  of  Reynolds  numbers.  A value  of  v ~ 10-6 
m2/ sec  was  assumed  in  the  computations  of  Re  shown  in  the  figure. 


Figure  2.6:  Schematic  transverse  force  and  corresponding  response  power  spectra  for 
steady,  planar  flow  around  a smooth,  transversely  unconstrained  cylinder  (Ref.  74). 

The  spectra  given  by  the  authors  were  only  for  the  conditions  after  lock-in.  For 
those  conditions,  Fig.  2.6  shows  that  the  vortex-shedding  frequency  displays  a depen- 
dency on  Reynolds  number  as  for  the  fixed  cylinder.  As  the  vortex-shedding  frequency 
moves  further  from  the  system  natural  frequency  the  response  at  the  natural  frequency 
decreases  while  the  response  at  the  vortex-shedding  frequency  increases.  In  a plot  of 
vortex-shedding  frequency  versus  velocity  their  results  show  that  the  vortex-shedding 
frequency  follows  the  St  = 0.2  relationship  except  in  the  neighborhood  of  the  natural 
frequency  where  lock-in  occurs. 
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It  is  known  that  the  transverse  motion  of  the  cylinder  has  a significant  effect  on  the 
vortex  correlation  and  consequently  on  the  transverse  force.  It  has  been  demonstrated 
by  Koopman  (1967)  and  Toebes  (1969)  that  the  vortex  span-wise  coherence  length  for 
a transversely  unconstrained  cylinder  exposed  to  a steady,  planar  current  increases 
drastically  with  increasing  amplitudes  of  oscillations.  However,  the  author  has  not 
been  able  to  locate  information  on  lift  coefficients  for  this  flow  situation. 

2.3  Oscillatory,  Planar  Flow 

2.3.1  Transversely  Constrained  Cylinder  in  Oscillatory,  Planar  Flow 

In  oscillatory,  planar  flow  around  rigid,  stationary  cylinders  both  the  acceler- 
ation and  relative  velocity  of  the  free  stream  are  constantly  changing  with  time, 
but  are  uniform  along  the  cylinder.  The  dimensionless  groups  that  characterize  this 
flow  situation  are  the  Reynolds  Number  ( Re  = UmaxD/v)  and  Keulegan- Carpenter 
Number  (KC  = UmaxT/D  = 2nA/D).  Combinations  of  these  groups  such  as 
/?  = Re/I<C  have  also  been  used  to  correlate  and  present  experimental  data  with 
varying  degrees  of  success.  Several  different  lift  coefficients  have  been  used  ( CLmax  = 
2 Fymax/pDUnax,  C Lrms  ~ ^Fyrms / pDU^^,  f-'Lrma  ‘^•FyTTna/ pDUrms,  etc.)  to 
present  results  (often  without  specifying  which  was  used).  Since  the  transverse  force 
had  been  found  in  nature  to  exhibit  some  degree  of  irregularity,  many  investigators 
have  commonly  used  an  rms  lift  coefficient  (Cxrm*)  and  so  is  this  review. 

The  investigations  for  this  flow  were  conducted  by  oscillating  the  flow  past  a fixed 
cylinder  or  by  oscillating  a rigid  cylinder  in  still-water.  Some  very  good  work  has 
been  done  in  this  area  resulting  in  the  largest  data  set  for  Cl  of  any  of  the  categories 
outlined  in  Table  2.1.  In  an  attempt  to  compare  the  results,  the  author  compiled 
rms  lift  coefficient  data  for  smooth  cylinders  from  a number  of  these  investigations 
and  plotted  them  versus  KC.  Some  of  the  data  were  taken  holding  Re  constant 
while  other  investigators  maintained  /3  constant.  These  data  are  plotted  separately 
in  Figs.  2.7  and  2.8.  Those  data  not  reported  in  terms  of  Cxrm4  = 2Fyrm3/ pDU^ax 
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were  converted  to  this  definition  prior  to  plotting.  The  data  for  Chaplin  (taken 
from  Bearman,  1988b)  and  Justesen  (1989)  in  Fig.  2.7  had  to  be  converted.  Other 
investigators’  data  (such  as  Bearman  et  al .,  1984;  Maull  and  Milliner,  1987)  are  not 
shown  since  they  presented  their  data  in  a manner  that  would  be  difficult,  if  not 
impossible,  to  convert  to  the  coordinates  used  in  these  figures.  In  Fig.  2.7  recent 
data  from  Longoria  et  al.  (1991)  are  presented  along  with  those  for  Sarpkaya  (1990a) 
and  Skomedal  (1989).  Both  plots  show  that  the  lift  coefficient  is  a maximum  at  KC 
between  10  and  12.  There  is  surprisingly  good  agreement  among  the  data  within 
each  of  the  two  plots.  The  fact  that  both  plots  have  the  same  shape  and  magnitude 
means  that  for  oscillatory,  planar  flow  the  lift  coefficient  depends  primarily  on  the 
Keulegan- Carpenter  number. 


2 - 


(0 
E 

ft. 

O i _ 


Symbol 

Ref# 

Re 

o 

25 

(0.3-1)  x 104 

▲ 

60 

(2-20)  x 104 

T 

49 

7 x 104  - 105 

• 

61 

(1-2)  x 105 

■ 

42 

2.5  xIO5 

A 

42 

5 x 10s 

* 

42 

7.5x105 

□ 

42 

1.0  x 106 

* * 

□ 

■TS 


7 

I AA flA  I 


T • 


J_L 


± 


o 

8 


10 


1 


100 


KC 


Figure  2.7:  Lift  coefficient  versus  Keulegan- Carpenter  number  for  oscillatory,  pla- 
nar flow  around  a smooth,  transversely  constrained  cylinder  (for  constant  Reynolds 
number,  Re). 


The  basic  nature  of  the  flow  in  this  category  depends  on  the  period  of  time  for 
which  the  flow  continues  in  one  direction  before  it  reverses.  If  the  period  is  very  short, 
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Figure  2.8:  Lift  coefficient  versus  Keulegan- Carpenter  number  for  oscillatory,  pla- 
nar flow  around  a smooth,  transversely  constrained  cylinder  (for  constant  frequency 
parameter,  /?). 


there  will  not  be  sufficient  time  for  the  vortices  to  form  before  the  flow  reverses.  If 
the  period  is  very  long,  the  flow  will  be  quasi-steady  and  will  have  the  character  of 
streaming  flow,  first  in  one  direction  and  then  in  the  other.  For  flow  periods  between 
these  extremes  some  downstream  vortex  effects  will  occur.  In  nondimensional  terms, 
it  was  found  that,  the  parameter  that  best  determines  the  general  character  of  vortex- 
shedding  from  a circular  cylinder  in  oscillatory  flows  is  the  ratio  of  how  far  a fluid 
particle  moves  in  one  half-cycle  to  the  characteristic  cylinder  dimension  that  the 
particle  flows  past.  This  nondimensional  term  is  what  is  known  by  the  Keulegan- 
Carpenter  number  aKC  = UmaxT/D” . 

Since  the  early  seventies  many  vortex-shedding  flow  visualization  studies  have 
been  conducted  for  this  flow  situation.  Most  of  these  investigations  were  carried  out 
to  further  the  understanding  of  the  mechanisms  and  thus  help  in  the  development 
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of  numerical  models  for  predicting  transverse  force.  Detailed  descriptions  of  these 
processes  have  been  given  by  Bearman  (1988b),  Sarpkaya  (1976a),  Skomedal  et  al. 
(1989)  and  Williamson  (1985).  Even  though  there  is  general  agreement  on  the  nature 
of  the  vortex  shedding  processes  there  are  differences  among  researchers  regarding 
the  details.  The  subjectivity  in  the  flow  visualization  techniques  used  to  measure  the 
conditions  under  which  vortices  are  shed,  the  number  shed,  etc.  probably  accounts 
for  many  of  the  differences. 

Experimentally,  the  frequency  of  the  vortex-shedding  is  computed  from  the  spec- 
tral analysis  of  the  transverse  force.  It  is  surprising  that  more  investigators  have  not 
reported  transverse  force  spectra.  The  spectra  that  have  been  reported  (Justesen, 
1989,  KC  = 1.7  to  15.6;  Bearman  and  Hall,  1987,  I<C  = 36.13;  McConnell  and  Park, 
1982b,  KC  = 37.7),  are  shown,  schematically,  in  Fig.  2.9.  The  intent  of  this  plot 
is  to  illustrate  the  behavior  of  the  frequency  and  not  the  magnitude  of  the  spectral 
density  of  the  transverse  force.  The  results  are  by  no  means  conclusive,  but  the  vor- 
tex energy  for  this  flow  appears  to  cluster  around  one  of  the  harmonics  of  the  in-line 
driving  (or  flow  oscillation)  frequency  with  moderate  energy  in  the  surrounding  har- 
monics, depending,  primarily,  on  the  value  of  KC.  As  KC  increases  the  frequency  of 
vortex-shedding  increases. 

Correlation  measurements  have  also  been  made  for  a fixed  cylinder  in  oscillatory, 
planar  flow,  but  unfortunately,  only  those  results  for  correlation  length  measured  by 
Obasaju  et  al.  (1988)  have  been  reported.  They  measured  the  span-wise  correlation 
of  vortex-shedding  for  a range  of  KC  from  4 to  55.  Their  results  show  that  the 
correlation  length  does  not  decrease  monotonically  with  increasing  KC.  The  highest 
correlation  length  was  obtained  at  I<C  = 10  (see  Fig  2.10)  and  at  this  value  it  was 
approximated  by  4.5 D.  They  also  found  that  for  KC  >30  the  correlation  is  no  longer 


sensitive  to  KC. 
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Figure  2.9:  Schematic  transverse  force  power  spectra  for  oscillatory,  planar  flow 
around  a smooth,  transversely  constrained  cylinder  for  various  Keulegan-Carpenter 
number,  KC  (Refs.  2,  32,  47). 


22 


2.3.2  Transversely  Unconstrained  Cylinder  in  Oscillatory,  Planar  Flow 

The  work  conducted  for  this  case  can  be  put  into  two  categories;  1)  that  which  is 
concerned  with  the  effect  of  the  transverse  motion  on  the  transverse  force  (McConnell 
and  Park,  1982a,  b;  Sarpkaya  and  Rajabi,  1979)  and  2)  that  which  is  concerned  with 
determining  the  conditions  under  which  transverse  motion  can  be  excited  (Bearman 
and  Hall,  1987;  Sarpkaya  and  Rajabi,  1979;  Sumer  and  Fredspe,  1988,  1989).  In 
general,  when  the  cylinder  is  allowed  to  have  transverse  motion,  the  vortex-shedding 
frequency  as  well  as  the  strength  of  the  vortices  is  modified.  For  this  situation  KC 
alone  is  no  longer  adequate  to  characterize  the  transverse  loading.  This  was  clearly 
demonstrated  by  Sumer  and  Fredsqe  (1988,  1989)  in  their  experiments  that  covered 
a wide  range  of  KC  ( KC  = 10  to  100)  and  a large  range  of  reduced  velocity,  VT  = 
Umax/Dfn  = 2n(A/ D)(fw  / fn)  = KC(fw/fn).  Their  results  demonstrated  that,  at 
least  for  the  range  of  KC  tested,  both  Vr  and  KC  are  necessary  to  describe  the 
behavior  of  the  transverse  response.  Their  plots  of  transverse  motion  due  to  vortex- 
shedding,  i.e.,  yrmj  versus  VT  for  constant  values  of  KC  show  that  for  lower  KC  (up 
to  « 10)  a single  spike  exists.  As  KC  increases,  the  number  of  spikes  increases  until 
at  KC  ~ 100  the  response  versus  Vr  is  flat  for  6 < VT  < 11.  The  reader  is  referred  to 
their  paper  for  a detailed  interpretation  of  this  behavior. 

To  the  author’s  knowledge,  McConnell  and  Park  (1982a)  are  the  only  investigators 
reporting  lift  coefficients  for  this  situation.  Their  results  showed  that  the  lift  coeffi- 
cient increased  up  to  twice  that  for  a fixed  cylinder.  Others  presented  only  structure 
response  information  such  as  yrms/D  or  ymax/D  versus  Vr  or  fw/fn  ratio. 

As  shown  in  Fig.  2.1,  when  transverse  motion  exists  there  can  be  an  interac- 
tion between  the  response  and  the  transverse  force.  The  level  of  this  interaction  is 
very  much  dependent  on  how  close  the  stationary  vortex-shedding  frequency  is  to 
the  system  s natural  frequency.  In  order  to  study  this  phenomenon  the  transverse 
force  and  response  should  be  measured  simultaneously.  Of  the  literature  reviewed, 
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only  McConnell  and  Park  (1982b)  reported  spectra  for  both  transverse  force  and 
the  corresponding  response.  Schematized  versions  of  these  spectra  are  presented  in 
Fig.  2.11.  The  frequency  for  the  force  and  response  spectra  has  been  normalized  by 
the  cylinder  natural  frequency  in  still-water.  Their  results  showed  for  4.4  < Vr  < 6.6 
and  fn/fd  = 6.22  lock-in  occurred,  i.e.,  /r//„  = fv/fn  = 1.  This  show  that  it  is 
not  necessarily  for  /„/ fj  to  be  an  integer  to  have  a lock-in.  Their  results  also  show 
(at  least  qualitatively)  the  influence  of  the  transverse  motion  on  the  frequency  of 
vortex-shedding. 


Figure  2.11:  Schematic  transverse  force  and  corresponding  response  power  spectra 
for  oscillatory,  planar  flow  around  a smooth,  transversely  unconstrained  cylinder  for 
various  reduced  velocities,  Vr  (Ref.  47). 

Data  or  information  on  the  vortex  span-wise  coherence  length  for  this  category 
was  not  available  until  recently  when  Kozakiewicz  et  al.  (1991)  measured  and  reported 
the  correlation  length  for  a vibrating  cylinder  near  a wall.  The  study  is  not  directly 
related  to  the  work  undertaken  here  as  it  was  mainly  concerned  with  the  wall  effect 
on  the  vortex  correlation;  however,  some  insight  can  be  gained  from  the  results. 
Their  results  for  the  largest  gap/diameter  ratio  (that  is  G/ D = 2.3)  are  reproduced 
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in  Fig.  2.12.  As  shown  in  this  figure  the  correlation  between  vortices  is  largest  at 
KC  = 6 (as  compared  to  KC  = 10  for  a fixed  cylinder  [Obasaju  et  al.,  1988])  and 
like  the  case  of  fixed  cylinder  is  mainly  dependent  on  KC.  The  results  of  their  study 
also  showed  that  for  a fixed  KC  the  larger  the  amplitude  of  transverse  motion  the 
larger  the  correlation. 


Figure  2.12:  Span-wise  coherence  length  for  transversely  unconstrained  cylinder  near 
a wall  in  oscillatory,  planar  flow  (Ref.  39). 

2.4  Oscillatory,  Nonplanar  Flow 

The  final  flow  configuration  to  be  considered  in  this  review  is  that  of  nonplanar, 
oscillatory  relative  flow  around  rigid,  circular  cylinders.  Very  little  work  has  been  done 
with  these  flow  conditions  as  can  be  seen  from  Table  2.1.  The  need  for  information  on 
the  effects  of  flow  nonuniformity  and  transverse  motion  on  transverse  force  and  vortex 
span-wise  coherence  length  was  the  motivation  for  the  work  of  this  dissertation. 

Both  transversely  constrained  and  unconstrained  cylinders  will  be  discussed  in 
this  section.  Surface  waves  acting  on  vertical  cylinders  are  also  included  in  this 
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category,  but  will  be  treated  separately  due  to  the  effects  of  the  vertical  component 
of  water  particle  velocity  on  the  transverse  force. 

The  question  is,  what  effect  does  flow  nonuniformity  have  on  the  transverse  force? 
More  specifically,  what  effect  does  it  have  on  the  span-wise  coherence  of  the  vortices? 
Does  the  coherence  length  of  the  vortices  diminish  with  nonuniformity  as  a result  of 
KC  gradients  along  the  cylinder  or  do  the  higher  energy  vortices,  associated  with  the 
regions  of  higher  velocity  and  vorticity,  dominate  the  shedding  process?  Also,  what 
effect  does  the  transverse  motion  of  the  structure  have  on  the  vortex-shedding  process 
and  transverse  force? 

Chakrabarti  et  al.  (1983)  and  Chakrabarti  and  Cotter  (1984)  conducted  exper- 
iments on  an  articulated  cylinder  where  the  top  was  oscillated  sinusoidally  while 
constraining  the  transverse  motion.  They  measured  local  and  total  forces  on  the 
cylinder  and,  among  other  things,  investigated  the  frequency  and  magnitude  of  the 
transverse  force  for  a range  of  KC  values.  They  used  a five  term  Fourier  series  (sim- 
ilar to  Mercier,  1973)  to  represent  the  measured  transverse  force,  defining  a different 
lift  coefficient  for  each  term  in  the  series.  These  lift  coefficients  are  reproduced  in 
Fig.  2.13  but  will  be  discussed  in  the  following  section  on  waves.  They  did  not  calcu- 
late the  phase  angles  associated  with  the  different  harmonics.  However,  they  advised 
the  use  of  a random  phase  angle  to  calculate  the  transverse  force.  Their  results  also 
show  that  the  dominant  frequencies  of  vortex-shedding  are  dependent  on  KC  and 
cluster  around  the  driving  frequency  for  the  lower  values  of  KC.  As  KC  increases, 
the  frequencies  cluster  around  a value  twice  the  driving  frequency  but  became  less 
organized  at  higher  I\C.  The  trend  continued  but  the  basic  assumption  of  clustering 
to  one  of  the  multiples  of  the  driving  frequency  seemed  to  break  down.  The  author 
was  unable  to  find  any  published  data  on  vortex  span-wise  coherence  length  for  this 
flow  category.  Likewise,  no  information  was  found  on  the  effect  of  transverse  motion 
on  transverse  force  and  vortex  span-wise  coherence  length. 
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Figure  2.13:  Lift  coefficient  harmonics  versus  Keulegan- Carpenter  number,  KC, 
— curve  fit  of  data  from  harmonically  oscillated  articulated  cylinder  in  still  water, 
• data  from  waves  impinging  on  a smooth,  fixed  vertical  cylinder  (Refs.  12,  13). 


2.5  Wave  Flows 

Flows  generated  by  regular,  unidirectional  surface  waves  impinging  on  fixed  ver- 
tical cylinders  are  discussed  separately  due  to  the  three  dimensionality  of  the  particle 
motion  and  its  potential  effects  on  the  transverse  force.  Several  researchers  (Bearman 
et  al .,  1985;  Bearman,  1988b;  Bidde,  1971;  Chakrabarti  et  al. , 1976,  1983;  Chakrabarti 
and  Cotter,  1984;  Isaacson  and  Maull,  1976,  1977;  Sawaragi  et  al.,  1976,  1977;  Tprum 
and  Reed,  1982)  have  conducted  experiments  with  regular  waves  and  have  reported 
lift  coefficient  data  as  a function  of  KC.  Others  (e.g.  Maull  and  Kaye,  1988)  have 
observed  the  effect  of  transverse  motion  on  the  in-line  motion  and  transverse  forces  in 
regular  waves.  Random  waves  have  also  been  investigated  (see  Bearman  et  al.,  1985; 
Graham,  1987). 

Chakrabarti  et  al.  (1976)  measured  vortex-induced  transverse  forces  on  local 
sections  of  a fixed  vertical  cylinder  (same  articulated  cylinder  they  oscillated  in  still 
water)  in  regular  waves.  They  used  the  same  series  technique  as  outlined  in  the  previ- 
ous section  to  obtain  the  lift  coefficients.  Data  reported  for  the  first  three  harmonics 
of  lift  coefficient  are  plotted  in  Fig.  2.13  along  with  curves  representing  (i.e.  poly- 
nomial curve  fits)  data  for  the  articulated  tower  discussed  above  under  oscillatory, 
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nonplanar  flows.  Note  that  Cm)  and  Ci( 3)  from  the  articulated  tower  experiments 
are  lower  than  the  corresponding  values  for  waves  while  Cl( 2),  which  is  the  larger 
of  the  coefficients,  is  about  the  same  for  waves  and  the  oscillated  articulated  tower. 
More  data  and  comparisons  are  needed  before  conclusions  can  be  drawn  for  this  flow 
situation. 

The  experiments  performed  in  this  category  (wave  flows)  covered  a wide  range 
of  apparatus  sizes  ranging  from  the  25  mm  (1  in)  diameter,  51  mm  (0.167  ft)  long 
cylinder  (used  by  Sawaragi  et  al.,  1976)  to  the  0.5  m (19.7  in)  diameter,  10  m (32.8  ft) 
long  cylinder  (Bearman,  1988b).  Since  the  transverse  force  was  found  to  exhibit  some 
degree  of  irregularity,  many  of  the  investigators  presented  their  data  in  terms  of  an 
rms  lift  coefficient.  As  in  the  case  of  planar  flows  discussed  above  there  are  differences 
in  Cxrma  used  by  the  different  investigators.  Once  again  the  author  compiled  available 
data  for  Cltuis  and  plotted  them  versus  KC  (see  Fig.  2.14).  Those  for  Bidde  (1971) 
were  not  included  since  they  are  defined  in  terms  of  C/,max.  In  this  figure,  Bearman’s 
data  are  based  on  an  rms  velocity  from  measurements  near  the  location  of  his  in- 
strumented segment.  It  is  not  clear  what  definition  was  used  by  Isaacson  and  Maull 
(1976)  but  the  author  assumed  it  to  be  based  on  the  maximum  value  of  the  velocity 
at  the  instrumented  segment  as  computed  using  linear  wave  theory.  Sawaragi  et  al. 
(1976)  did  not  define  their  Cz,rm4  but  since  they  compared  their  results  with  Isaacson 
and  Maull  the  same  assumption  made  for  Isaacson  and  Maull  was  applied  to  them. 
As  would  be  expected  there  is  more  scatter  in  these  data  than  for  the  less  complex 
flows.  This  is  particularly  true  for  the  data  of  Sawaragi  et  al.  for  larger  values  of  KC 
and  kh  between  1.2  and  2.4.  Possible  reasons  for  such  scatter  include  the  large  differ- 
ences in  model  sizes,  differences  in  H/D,  and  differences  in  ellipticity  of  the  particle 
motion  (Bearman  et  al .,  1985  found  this  was  not  a factor  for  the  conditions  of  their 
experiments).  However,  it  is  interesting  to  note  that  the  lift  coefficient  is  a maximum 
at  about  the  same  KC  values  as  for  oscillatory,  planar  flow  {KC  between  10  and  12). 
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Figure  2.14:  Lift  coefficient  versus  Keulegan- Carpenter  number  for  regular  waves 
around  a smooth,  fixed  vertical  cylinder. 


Some  flow  visualization  tests  have  also  been  conducted  in  this  flow  category  to 
study  the  vortex-shedding  process  and  its  frequency.  Bearman  (1988b),  Isaacson  and 
Maull  (1976)  and  Sawargi  et  al.  (1976)  carried  out  visualization  tests  to  characterize 
the  frequency  of  the  vortex-shedding  process.  However,  due  to  the  subjectivity  in 
the  visualization  techniques  the  behavior  of  vortex-shedding  frequency  is  examined 
from  the  transverse  force  power  spectra.  The  spectra  obtained  by  Chakrabarti  et  al. 
(1976)  and  by  T0rum  and  Reed  (1982)  were  thus  compiled  and  plotted  in  Fig  2.15. 
From  these  spectra,  one  can  see  the  existence  of  multiple  frequencies  at  all  values 
of  KC,  but  the  dominant  frequencies  have  the  same  trend  as  in  oscillatory,  planar 
flows.  That  is,  at  lower  values  of  KC  the  dominant  frequencies  cluster  around  the 
wave  frequency.  As  KC  increases  to  ss  9 the  frequencies  cluster  at  twice  the  wave 
frequency,  then  at  three  times  the  wave  frequency  a s KC  increases  further,  etc. 


29 


The  author  was  unable  to  find  reported  data  on  span-wise  coherence  length  or 
transverse  motion  for  this  flow  category  (waves).  However,  as  mentioned  before,  some 
investigators  allowed  transverse  motion  in  their  studies.  These  investigators  were 
mainly  concerned  with  determining  the  conditions  under  which  transverse  motion 
can  be  excited  since  this  motion  can  cause  shut-down  of  operations  and  fatigue  of 
structural  members  (such  as  marine  risers).  In  an  attempt  to  quantify  such  conditions, 
the  maximum  values  of  yTma/D  (defined  as  Yrma/D)  from  several  investigators’  data 
were  plotted  versus  several  of  the  pertinent  parameters.  A plot  of  Yrma/ D versus  VT 
(see  Fig.  2.16)  was  found  to  be  a meaningful.  This  plot  shows  that  the  maximum 
values  of  the  response  {Yrmal D)  for  different  flow  conditions  fell  in  a relatively  narrow 
range  of  reduced  velocity,  VT. 

In  summary,  the  existing  data  on  the  vortex-induced  transverse  loading  and  re- 
sponse for  fixed  and  transversely  unconstrained,  smooth,  cylinders  were  reviewed  for 
the  different  flow  situations.  As  a result  of  this  review  one  can  say  there  have  been 
significant  advances  in  the  understanding  of  vortex-induced  forces  on  bluff  bodies 
for  some  types  of  flow,  especially  steady  and  oscillatory,  planar  flows.  In  spite  of 
these  advances,  however,  more  work  is  needed  before  the  effects  of  many  of  the  flow 
and  geometric  parameters  on  the  transverse  force  and  response  can  be  quantified  for 
the  design  engineer.  Lack  of  consistency  in  data  reporting  has  hindered  progress  in 
some  areas  while  the  difficulty  of  measuring  some  of  the  important  quantities  (such 
as,  span-wise  coherence  length  of  vortices)  has  limited  the  data  available  for  analy- 
sis. With  an  increasing  awareness  of  the  problems  associated  with  inconsistent  data 
reporting  and  with  rapid  advances  in  instrumentation  technology,  hopefully,  these 
problems  can  be  eliminated.  This  study  focuses  on  the  effects  of  flow  nonuniformity 
and  transverse  motion  on  vortex-induced  transverse  forces  (in  particular,  lift  coef- 
ficients and  frequency  of  vortex-shedding)  and  on  the  vortex  span-wise  correlation 
length. 
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Figure  2.15:  Schematic  transverse  force  power  spectra  for  regular  waves  around  a 
smooth,  fixed  vertical  cylinder  for  various  Keulegan-Carpenter  numbers,  KC  (Refs. 
12,  73). 
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Figure  2.16:  Maximum  transverse  response  for  various  flow  configurations  (□  steady, 
planar  flow),  (o  oscillatory,  planar  flow),  (•  wave  flows). 


CHAPTER  3 

EXPERIMENTAL  INVESTIGATIONS 

3.1  Scaling  Parameters  and  Model  Selection 

The  modelling  of  fluid-induced  dynamic  forces  exerted  on  bodies  immersed  in 
a viscous  fluid  has  always  presented  difficulties  with  regard  to  similarity  and  scale. 
This  is  especially  true  when  the  structure  is  such  that  it  is  free  to  move.  If  dimen- 
sional analysis  is  applied  to  the  particular  problem  of  a transversely  unconstrained 
bottom-hinged  cylinder  in  nonuniform  oscillatory  fluid  flow,  the  amplitude  of  trans- 
verse motion  can  be  expressed  in  terms  of  the  lift  coefficient,  Cl , reduced  velocity, 
Vr,  stability  parameter,  Ka,  and  added  mass  to  cylinder  mass  ratio,  m„/m4,  (see 
Appendix  A for  details),  i.e., 

y/D  = Y(CL,VT,K3,ma/m3)  (3.1) 

where 

CL  = CL(Re,  KC,Ty).  (3.2) 

If  the  cylinder  is  smooth  and  the  external  generated  turbulence  “ Ty ” in  the  flow  is 
small,  Eqn.  3.2  reduces  to 

CL  = CL(Re,KC).  (3.3) 

Consequently,  Eqn.  3.1  reduces  to 

y/D  — Y(Re,KC,Vr,K3,ma/m3).  (3.4) 

For  complete  similarity  the  values  of  all  the  parameters  should  be  the  same  for 
both  model  and  prototype,  i.e., 
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where  m and  p refer  to  model  and  prototype,  respectively. 

Of  these  similarity  groups  the  most  important  parameters  are  the  Reynolds  num- 
ber ( Re ),  Keulegan-Carpenter  number  ( KC ) and  reduced  velocity  (K-).  Hydrody- 
namic force  coefficients  for  large  values  of  Re  are  needed  in  full-sized  offshore  struc- 
tures applications.  These  conditions  are  difficult  if  not  impossible  to  produce  in 
laboratory  experiments.  Wave  tank  testing  generally  produces  Reynolds  numbers 
(based  on  structure  diameter)  up  to  approximately  5 x 104,  whereas  most  prototype 
structures  experience  Reynolds  numbers  well  beyond  10s,  i.e.,  in  the  upper  subcrit- 
ical,  critical  and  supercritical  regions  (according  to  steady  flow  principal  Reynolds 
number  flow  regimes  around  smooth  cylinders).  Wave  tank  testing  has  also  the  addi- 
tional problem  of  not  being  able  to  precisely  control  Reynolds  number.  Ideally  both 
Reynolds  number  and  reduced  velocity  must  be  scaled,  ft  is  practically  impossible 
to  achieve  both  a desired  Reynolds  number  and  reduced  velocity  at  the  same  time 
unless  a full  scale  structure  is  tested. 

For  an  average  prototype  articulated  tower,  such  as  the  “ELF”  loading  tower 
operating  in  the  North  Sea,  (diameter  = 4.5  m,  length  = 150  m and  water  depth 
= 135  m),  the  reduced  velocity,  VT  varies  between  0 and  12  for  a current  velocity 
varying  between  0 and  2.5  knots,  ft  can  also  be  as  high  as  20  for  a 100  year  design 
wave  (having  a 17  sec.  period  and  30  m height  according  to  Kirk  and  Jain,  1977). 
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In  laboratory  testing  it  is  important  to  cover  the  range  mentioned  above  for  both 
prototype  Reynolds  number  and  reduced  velocity.  Using  a constant  diameter  cylinder 
this  can  be  done  by  adjusting  either  the  natural  frequency  /„  or  the  velocity  range 
or  both.  In  the  tests  carried-out  in  this  investigation  the  natural  frequency  was  held 
approximately  constant. 

The  intent  of  this  investigation  was  to  study  the  effect  of  flow  nonuniformity  and 
transverse  motion  on  the  vortex-induced  transverse  loading.  In  order  to  make  the 
results  of  this  study  as  useful  as  possible  to  the  design  of  offshore  structures,  the 
model  was  designed  to  produce  ranges  of  similarity  parameters  as  close  as  possible  to 
those  experienced  by  prototype  structures  where  vortex- induced  loading  is  important. 
Other  geometric  parameters  such  as  the  length/diameter  ratio  (£/ D ) should  also  fall 
within  the  prototype  range.  A survey  of  the  available  literature  revealed  that  vortex- 
shedding  is  considerable  in  the  following  ranges  of  similarity  parameters:  3 x 104  < 
Re  < 106,  5.0  < KC  < 12.0,  1.0  < Vr  < 20.0  and  i/ D between  20  and  30. 

The  requirements  for  achieving  similarity  (i.e.,  equivalent  values  for  model  and 
prototype)  for  Reynolds  number  conflicts  with  the  conditions  required  for  similarity 
for  reduced  velocity  and  other  parameters.  For  example,  to  produce  high  values  of 
Reynolds  number  requires  high  values  of  flow  velocity  since  v is  fixed  and  the  cylinder 
diameter  is  constrained  by  the  ratio  IfD.  To  satisfy  reduced  velocity  similarity  the 
higher  flow  velocities  imposed  by  Reynolds  number  similarity  require  high  system 
natural  frequency  which  means  a very  light  structure  which  is  limited  by  required 
rigidity.  Also  higher  flow  velocities  means  higher  amplitudes  of  oscillation  or  higher 
driving  frequencies  which  are  constrained  by  the  available  driving  mechanisms.  Due 
to  these  complications,  the  design  of  the  cylindrical  model  used  in  these  experiments 
was  accomplished  with  the  aid  of  a computer  program  that  optimized  the  range  of 
the  parameters  within  the  limitations  imposed  by  the  drive  mechanism,  tank,  budget, 
etc.  In  other  words,  given  the  values  of  the  facilities  available  such  as  maximum  water 
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depth,  maximum  in-line  velocity  and  others,  the  program  selects  the  dimensions  that 
cover  ranges  of  the  similarity  parameters  overlapping  with  those  of  prototypes.  In 
this  study,  model  dimensions  of  0.1524  m diameter  and  3.05  m length  were  selected, 
i.e,  a ratio  of  t/D  = 20.  Such  a model  covered  ranges  of  Re  between  6.1  x 103  and 
1.3  x 105,  KC  between  2.5  and  9.5  and  VT  between  3.15  and  46. 

For  an  elastically-mounted  cylinder,  experiments  should  be  carried-out  at  constant 
K,  to  facilitate  comparison  of  maximum  response  amplitudes  for  different  conditions. 
In  this  investigation,  despite  the  changes  in  temperature  during  the  course  of  testing, 
K,  was  approximately  constant. 

The  ratio  of  the  fluid  added  mass  to  the  structural  mass  ( ma/m ,)  should  be 
modelled  properly.  In  air  ma/m,  = 0(1O-3)  and  thus  is  of  little  importance.  In 
water,  however,  the  ratio  ma/m,  = 0(1)  and  affects  both  the  maximum  amplitude 
of  motion  due  to  vortex-shedding  and  the  velocity  range  over  which  lock-in  occurs. 

3.2  Experimental  Set-up 

Very  little  information  exists  on  the  effect  of  flow  nonuniformity  and  structure 
transverse  motion  on  vortex-induced  transverse  loading  and  vortex  span-wise  coher- 
ence length.  To  investigate  such  effects  two  types  of  experiments  were  designed  and 
performed.  Both  experiments  were  performed  in  a circular  water  tank  with  a rigid 
articulated  tower  model.  In  the  first  type  of  experiment  the  model  was  constrained 
from  transverse  motion  while  in  the  second  type  it  was  free  to  move  in  the  transverse 
direction. 

As  discussed  before  in  Chapter  2,  there  are  two  different  ways  to  produce  a relative 
oscillatory  motion  between  a cylindrical  model  and  the  surrounding  fluid.  One  method 
is  to  oscillate  the  fluid  past  a stationary  cylinder  while  the  other  is  to  oscillate  the 
cylinder  past  the  stationary  fluid.  Kinematically  there  is  no  difference  between  the 
two  situations  when  viewed  from  the  appropriate  reference  frame.  Experimentally, 
there  are  differences  in  implementation  between  the  two  methods. 
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In  the  experiments  of  this  investigation,  an  articulated  cylinder  was  oscillated  in 
an  otherwise  still  fluid.  This  approach  was  taken  because  it  was  relatively  easy  to 
implement  and  because  it  allowed  more  precise  control  over  the  parameters.  It  also 
allowed  the  fluid-induced  forces  to  be  determined  more  directly.  In  other  words  one 
can  determine,  after  subtracting  the  inertial  force  due  to  the  mass  of  the  oscillating 
cylinder,  the  coefficient  of  added  mass  ( Cm ) instead  of  the  inertia  coefficient  (Cm,  i.e., 
1 + C'm)  since  no  horizontal  pressure  gradient  exists  in  the  flow  field.  On  the  other 
hand,  this  approach  has  certain  disadvantages  which  must  be  overcome.  These  are: 
i)  waves  and  free  surface  disturbances  can  be  created  by  the  oscillating  cylinder,  ii) 
the  drive  mechanism  can  transmit  vibrations  to  the  cylinder  and  surrounding  fluid 
and  iii)  the  inertia  force  due  to  the  mass  of  the  oscillating  body  has  to  be  accounted 
for  in  the  measured  force  signal. 

In  this  experiment,  several  ideas  and  designs  were  employed  to  overcome  these 
difficulties.  Due  to  the  tank  to  cylinder  diameter  ratio  (=  60)  wall  effects  were 
negligible  (in  terms  of  a blockage  ratio  it  was  1.6%).  This  is  well  below  the  6.0% 
limit  given  by  West  and  Apilt  (1982)  for  no  influence  from  blockage  on  the  Strouhal 
number.  Soft  and  porous  packing  materials  were  also  placed  at  water  surface  half- 
way between  the  cylinder  and  the  tank  wall  to  absorb  surface  disturbances  and  reduce 
wave  reflections  from  the  wall  (see  Fig.  3.3).  Two  heavy  steel  I-beams  were  placed 
across  the  tank  to  support  the  linear  driving  motor.  The  feed  back  control  system  for 
the  linear  drive  motor  was  initially  a source  of  high  frequency  vibration,  but  this  was 
minimized  by  using  a lowpass  analog  filter  in  the  feed  back  circuit. 

The  two  types  of  experiments  used  in  this  study  were  designed  to  obtain  the  in- 
formation needed  to  quantify  the  dependence  of  the  vortex-shedding  process,  vortex 
span- wise  coherence  length  and  transverse  force  on  the  flow  and  structure  parameters. 
All  tests  were  performed  in  a deep  circular  water  tank  9.15  m in  diameter  and  3.1  m 
high.  The  tower  model  was  a right  circular  rigid  aluminum  cylinder  with  a diameter 
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of  0.1524  m,  a thickness  of  3.0  mm  and  a length  of  3.05  m.  It  was  instrumented  with 
miniature,  quick  response,  (2.15  mm  diameter  bead)  thermistors  embedded  in  its  sur- 
face to  measure  the  frequency  of  vortex  shedding  along  the  tower.  These  thermistors 
were  located  at  the  leading  and  trailing  sides  of  the  tower  at  20  different  locations 
(total  of  10  thermistors  along  each  side).  All  the  wires  of  the  thermistors  were  run 
through  the  inside  of  the  tower.  The  signal  processing  circuits  for  the  thermistors  were 
mounted  inside  the  cylinder  near  the  top.  The  cylinder  was  attached  to  the  bottom  at 
the  center  of  the  tank  through  a low  friction  hinge  designed  to  allow  X and  Y motion 
while  constraining  cylinder  rotation  about  its  axis.  Its  top  was  attached  through  a 
vertically  sliding  shaft  to  a linear  drive  electric  motor  mounted  horizontally  between 
the  I-beams.  The  vertical  sliding  shaft  was  designed  to  account  for  the  change  in  the 
vertical  position  of  the  top  of  the  tower  during  its  rotation  about  the  bottom  hinge. 
Ultra  low  friction  linear  ball  bushings  for  the  vertical  sliding  shaft  were  mounted  in 
the  upper  end  plate  of  the  cylinder.  An  X — Y force  transducer  was  inserted  between 
the  top  of  the  tower  and  the  table  of  the  linear  drive  motor.  To  negate  the  need  to 
measure  forces  at  the  base,  a pin  joint  was  placed  between  the  tower  and  the  X — Y 
force  transducer  at  the  top.  Two  Linear  Displacement  Transducers  (LDT)  were  used 
to  measure  in-line  and  transverse  displacements  of  the  tower.  The  active  strokes  of 
the  in-line  and  transverse  LDT  were  0.66  m and  0.46  m,  respectively.  All  alignments 
were  made  prior  to  adding  water  to  the  tank. 

In  the  first  type  of  experiments,  tests  were  performed  with  the  configuration  shown 
in  Fig.  3.1.  An  in-line  simple  harmonic  motion  was  imposed  at  the  top  of  the  tower 
while  constraining  the  transverse  motion.  Quantities  measured  during  these  tests 
were  i)  in-line  position,  ii)  in-line  force,  iii)  transverse  force  and  iv)  vortex-shedding 
frequencies  from  the  20  thermistors  mounted  along  the  tower. 

The  second  experiment  type  configuration  shown  in  Fig.  3.2  used  the  same  tower 
model  with  the  exception  of  the  mechanism  that  attaches  the  cylinder  to  the  linear 
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drive  motor.  The  tower  was  driven  with  the  same  in-line  motion  as  in  the  first 
type  experiment,  but  in  this  case  it  was  allowed  to  respond  freely  in  the  transverse 
direction.  Quantities  measured  in  these  tests  included  i)  in-line  position,  ii)  in-line 
force,  iii)  transverse  motion  and  iv)  vortex- shedding  frequencies  from  the  20  different 
thermistors  along  the  tower.  The  water  depth  for  all  the  tests  (type  I and  II)  was 
2.65  m.  Photographs  of  the  general  setup  and  transverse  motion  mechanism  are 
shown  in  Fig.  3.3. 

3.3  Instrumentation  and  Calibration 

As  discussed  in  the  previous  section,  the  following  quantities  were  measured 

1.  In-line  tower  position. 

2.  In-line  force. 

3.  Transverse  force  (type  I experiment  only). 

4.  Transverse  response  (type  II  experiment  only). 

5.  Thermistor  signals. 

Other  signals  such  as  the  input  driving  frequency  were  measured  and  held  constant 
during  each  test.  Quantities  such  as  water  temperature  and  surface  wave  activity 
were  also  monitored. 

A block  diagram  of  the  measurement  system  used  is  shown  in  Fig.  3.4.  The 
instrumentation  used  to  measure  these  quantities  are  listed  below.  The  calibration 
procedures  used,  where  needed,  are  also  discussed. 

3.3.1  Frequency  Generator 

The  sinusoidal  motion  imposed  at  the  top  of  the  tower  was  generated  using  a 
Hewlett  Packard  (HP)  frequency  generator.  A Hewlett  Packard  (HP)  frequency 
counter  in  parallel  with  the  frequency  generator  was  also  used  to  continuously  monitor 
the  input  driving  frequency  ( fd ) during  each  run. 
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1 Tank 

2 Articulated  Tower  (10*  x 6"<))) 

3 Linear  Motor  Drive 

4 Force  Transducer 

5 Trandsucer  Coupling 

6 Upper  Cylinder  Bearing  Guide 

7 Linear  Bushing  Bearing 

8 Vertical  Sliding  Shaft 

9 Universal  Joint 

10  Lumber  Frame  Cover 


Figure  3.1:  Schematic  diagram  of  transversely  constrained  experiment  set-up. 
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1 Tank 

2 Articulated  Tower  (10'  x 6"<{>) 

3 Linear  Motor  Drive 

4 Force  Transducer 

5 Trandsucer  Coupling 

6 Upper  Cylinder  Bearing  Guide 

7 Linear  Bushing  Bearing 

8 Vertical  Sliding  Shaft 

9 Spacer 

10  Universal  Joint 

11  Linear  Bushing  Bearing 

12  Lumber  Frame  Cover 


Figure  3.2:  Schematic  diagram  of  transversely  unconstrained  experiment  set-up. 
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Figure  3.3:  Photographs  of  experimental  set-up. 
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lln-llne  Position 


Analog/Digital  Processing 


In-line  Position 


Linear  Motor  & Controller 

Figure  3.4:  Block  diagram  of  measurement  system. 


3.3.2  Force  Transducers 

Due  to  the  high  cost  of  commercial  X — Y (horizontal  components)  load  cells, 
the  one  used  in  these  experiments  was  designed,  constructed  and  tested  in  the  de- 
partment laboratory.  The  transducer  was  a beam-type  load  cell  (see  Appendix  B, 
Fig.  B.l)  made  of  304  stainless  steel  with  overall  dimensions  3.8  x 3.8  x 7.5  cm.  It 
was  constructed  as  if  two  I-beam  were  placed  flange-to-flange  so  that  their  webs  make 
a “90°”  angle.  The  thickness  of  each  web  which  was  required  to  withstand  the  maxi- 
mum loads  anticipated  in  its  direction  was  calculated  using  Hooke’s  law,  with  a safety 
factor  of  1.5.  With  the  scantlings  selected,  the  transducer  had  a resolution  of  0.2  N 
in  each  direction.  Each  web  was  instrumented  with  four  (two  on  each  side)  active 
350  ohm  strain  gages  (with  a gage  factor  of  2.03  ± 1%).  These  gages  formed  the 
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elements  in  a standard  four-arm  wheatstone  bridge  circuit  and  provided  both  high 
sensitivity  and  temperature  compensation.  The  strain  gage  bridge  was  connected  to  a 
strain  gage  amplifier,  the  Omega  DMD  465  (see  Appendix  B,  Fig.  B.2  and  Table  B.l 
for  schematic  and  specifications).  A pulley  system  was  developed  to  calibrate  each 
channel  of  the  transducer.  The  calibration  was  carried  out  by  varying  the  loads  on 
one  channel  in  both  directions  (i.e.,  in  tension  and  compression).  Loads  were  added 
in  increments  then  removed  in  increments  to  test  for  hysterisis.  At  the  same  time  a 
constant  load  was  maintained  on  the  other  channel  to  examine  the  side  load  effect. 
This  was  repeated  for  different  constant  side  loads.  The  calibration  curves  are  shown 
in  Appendix  B,  Figs.  B.3,  to  B.5.  Note  that  the  two  channels  are  very  well  isolated 
(i.e.,  loads  in  one  direction  have  minimum  influence  on  the  transducer  output  in  the 
other  direction). 

3.3.3  Linear  Displacement  Transducers 

The  tower’s  in-line  and  transverse  motions  were  measured  using  two,  MTS  Tem- 
posonics  II,  LDT  (Linear  Displacement  Transducers)  systems  with  AOM  (Analog 
Output  Modules),  see  Appendix  B,  Fig.  B.6  for  a schematic  setup  and  Table  B-2  for 
the  specifications  of  the  AOM.  Both  LDTs  were  calibrated  by  giving  the  tower  known 
displacements  measured  with  a rule  while  the  output  from  the  A/D  data  acquisition 
board  was  recorded.  A statistical  analysis  (average  and  standard  deviation)  of  the 
recorded  calibration  data  at  each  position  was  then  carried-out.  This  procedure  was 
repeated  for  the  different  positions  around  the  vertical  position  of  the  tower  i.e.,  in 
the  positive  and  negative  directions.  The  resultant  data  gave  straight  line  calibration 
curves  for  each  transducer,  as  shown  in  Appendix  B,  Fig.  B.7  and  Fig.  B.8. 

3.3.4  Thermistors 

The  vortex  shedding  frequencies  were  measured  by  using  miniature,  quick  re- 
sponse, thermistors  embedded  in  the  surface  of  the  tower.  Dimensions,  thermal  and 
electrical  properties  of  these  thermistors  are  given  in  Appendix  B,  Table  B-3.  The 
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principal  is  similar  to  that  of  a hot  film  anemometer.  That  is,  the  overheated  ther- 
mistor is  cooled  by  the  flow  of  the  fluid  past  its  2.15  mm  diameter  bead  which  in  turn 
reduces  its  electrical  resistance.  The  circuit  designed  for  processing  the  signals  of 
these  thermistors  is  given  in  Appendix  B,  Fig.  B.9.  The  concept  of  using  thermistors 
to  measure  frequency  of  vortex-shedding  was  tested  using  a thermistor  embedded  in 
a 0.065  m diameter  and  0.61  m long  PVC  cylinder.  The  cylinder  was  then  placed  in 
a steady  flow  on  a flow  visualization  table  as  shown  in  Appendix  B,  Fig.  B.10.  Three 
different  flow  velocities  with  Reynolds  numbers  of  4.6  x 103 , 9.7  x 104  and  5 x 104  were 
used.  The  thermistor  output  signals  were  sampled  at  50  Hz  and  a spectral  analysis 
performed  on  each  signal.  The  spectra  of  the  signals  (see  Appendix  B,  Fig.  B.ll) 
showed  a single  large  spike  (relative  to  other  frequency  components  in  the  spectrum). 
The  Strouhal  number,  St  at  these  spikes  was  found  to  equal  0.23,  0.188  and  0.123, 
respectively  for  the  Reynolds  numbers  above. 

3.3.5  Lowpass  Filters 

Analog  filters  must  always  be  used  with  caution  since  they  can  produce  undesir- 
able as  well  as  desirable  effects.  The  undesirable  effects  are  in  the  form  of  signal  phase 
shift  and  amplitude  attenuation.  To  minimize  these  effects,  a second  order,  Butter- 
worth,  lowpass  filter  was  designed  for  the  in-line  position  signal  to  reduce  the  noise  in 
the  signal  to  the  servo  controller.  The  cut-off  frequency  was  carefully  selected  for  the 
filter  after  considering  the  useful  driving  frequency  range.  Since  the  highest  driving 
frequency  was  1.0  Hz,  a cut-off  frequency  of  3.0  Hz  for  the  feed  back  loop  filter  was 
deemed  sufficiently  high.  A schematic  diagram  of  this  filter  is  shown  in  Appendix  B, 
Fig.  B.12.  Two  other  lowpass  filters  were  used  for  the  in-line  and  transverse  forces. 
Their  cut-off  frequencies  were  selected  after  examination  of  the  unfiltered  signals  from 
exploratory  runs  for  all  the  proposed  driving  frequencies.  Spectral  analyses  were  then 
performed  on  the  measured  signals.  The  results  of  these  spectra  showed  no  strong 
frequency  components  (even  for  the  highest  driving  frequency,  1.0  Hz)  above  15.0  Hz 
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for  the  in-line  and  20.0  Hz  for  the  transverse  force  signals.  Accordingly  17.0  Hz  and 
22.0  Hz  cut-off  frequencies  for  the  in-line  force  and  transverse  force  signals  were  se- 
lected. A sample  output  of  the  transverse  force  time  series  and  its  power  spectrum 
before  and  after  filtering  is  given  in  Appendix  B,  Fig.  B.13. 

Other  more  standard  instruments  used  in  these  experiments  included:  A Normag 
DC-linear  drive  motor  (specifications  given  in  Appendix  B,  Table  B.4);  DC-servo 
controller  (Moog,  Model  82  — 300);  glass  bulb  mercury  thermometer  and  a 12— bit, 
multi-function,  high  speed  A/D  Metrabyte/Asyst/DAC  expansion  board  with  analog 
input  accuracy  of  0.01%  and  a variable  input  voltage  range.  This  board  was  installed 
in  an  IBM  compatible  personal  computer. 

3.4  Data  Reduction 

As  mentioned  before,  two  types  of  experiments  were  performed.  Type  I was  with 
the  transverse  motion  constrained  and  type  II  was  with  a cylinder  free  to  move  in  the 
transverse  direction. 

Tests  were  run  for  different  values  of  in-line  amplitude  and  driving  frequency.  In 
the  first  type  experiment,  the  cylinder  driving  frequency  (fd)  was  varied  from  0.4  Hz 
to  1.0  Hz.  At  most  frequencies  the  amplitude  of  oscillation  varied  from  0.1  m to 
0.28  m.  A total  of  86  runs  were  made  during  this  experiment.  Of  these  8 runs  were 
performed  to  test  the  repeatability.  In  each  run,  in-line  and  transverse  forces,  in-line 
position  and  the  signals  from  the  10  thermistors  were  recorded.  During  the  test  runs 
GLOBALLAB  software  was  used  to  perform  preliminary  analysis  of  the  measured 
data.  The  tower  driving  frequency  (/<*)  and  amplitude  (a)  were  also  monitored  during 
each  run  using  a frequency  counter  and  oscilloscope. 

In  the  type  II  experiment,  where  transverse  motion  was  allowed,  the  driving 
frequency  (/<;)  was  varied  from  0.15  Hz  to  0.8  Hz.  At  each  frequency  the  amplitude 
was  varied  from  0.076  m to  0.28  m.  A total  of  190  runs,  including  36  runs  to  test 
repeatability,  were  made  during  this  experiment.  In  each  run  the  in-line  force,  in-line 
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and  transverse  position  and  the  thermistor  signals  were  measured  and  recorded.  To 
reduce  the  scatter  in  the  measured  data  of  in-line  force  and  transverse  motion  some 
longer  test  runs,  where  only  in-line  force,  transverse  motion  and  the  signal  of  vortex 
shedding  from  the  top  thermistor,  were  made. 

Since  the  Direct  Memory  Access  (DMA)  page  registers  cannot  be  incremented  by 
the  controller,  the  maximum  data  area  available  (64K,  a page,  for  32,  767  conversions) 
was  used  to  acquire  data  for  each  run.  This  resulted  in  a duration  of  63  sec  for  each 
recorded  signal  in  the  first  type  experiment  where  13  data  channels  were  acquired  and 
91  sec  (273)  sec  for  the  second  type  experiment  when  9 or  (3)  channels  were  used. 

The  maximum  obtainable  in-line  amplitude  was  found  to  vary  with  the  driving 
frequency  (fj).  The  intent  was  to  obtain  as  wide  a range  of  Keulegan-Carpenter  ( KC ) 
and  Reynolds  ( Re ) numbers  as  possible.  The  range  covered  by  these  parameters  and 
others  during  both  experiments  together  with  a summary  of  the  test  conditions  are 
shown  in  Table  3.1.  The  natural  frequency  (/n)  and  damping  ratio  (£)  included  in 
Table  3.1  were  obtained  by  measuring  the  frequency  of  free  oscillation  in  still  water. 
The  measured  value  of  (/n)  compared  well  with  the  analytical  value  obtained  from 
the  free  vibration  response  predicted  by  a computer  program  developed  by  Omar 
and  Sheppard  (1991)  for  predicting  response  of  articulated  towers  under  the  action 
of  wind,  current  and  waves. 
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Table  3.1:  Test  conditions  for  experiments  I and  II. 


Item 

Exp.  I 

Exp.  II 

minimum  fd 

0.5  Hz 

0.15  Hz 

maximum  fd 

1.0  Hz 

0.8  Hz 

minimum  a 

0.085  m 

0.08  m 

maximum  a 

0.283  m 

0.29  m 

temp,  range 

83  - 88°  F 

60  - 74°  F 

minimum  Re 

6.1  x 103 

8.4  x 103 

maximum  Re 

2.15  x 104 

1.3  x 10s 

minimum  KC 

2.6 

2.4 

maximum  KC 

8.65 

9.35 

minimum  (3 

1792 

3151 

maximum  (3 

3692 

18336 

minimum  Vr 

2.66 

maximum  VT 

36.10 

fn 

0.15  Hz 

0.12 

water  depth 

2.625  m 

2.625  m 

No.  of  channels 

13  or  14 

3 or  9 

sampling  frequency 

40  Hz 

40  Hz 

CHAPTER  4 

MATHEMATICAL  MODELS 


4.1  In-line  Force 


The  following  equation  of  motion  was  used  to  reduce  the  in-line  force  data: 


4 nO(t)  — X!  -^total  = + Md  + Mfx, 


(4.1) 


where 


Io 

K 

M total 

Mg 

Md 

Mfx 

m 


= total  mass  moment  of  inertia  about  the  bottom  hinge, 

= cylinder  mass  moment  of  inertia, 

= added  mass  moment  of  inertia, 

= total  moment  about  bottom  hinge, 

= moment  due  to  tower  weight  and  buoyancy, 

= moment  due  to  drag, 

= moment  due  to  linear  drive  motor  and 
= in-line  deflection  angle. 

the  tower  shown  in  Fig.  4.1  the  moments  Mg , Md  and  Mfx  are  given  by 


t=i 


•d\  / co*  6 


(4.2) 


(4.3) 

(4.4) 


Mfx  = RL, 


where  Ac  — 7tD2/4  is  the  cross  sectional  area  of  the  tower  element.  The  added  mass 
moment  of  inertia,  /'  can  be  expressed  as 
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The  motion  of  the  tower  was  imposed  by  a horizontal  linear  drive  motor  mounted 
directly  above  the  tower  as  shown  in  Fig.  4.1.  The  position  of  the  moving  table  of 
the  motor  and  thus  the  upper  end  of  the  shaft  connected  to  the  tower  was  monitored 
by  a LDT  (Linear  Displacement  Transducer).  Due  to  the  relatively  small  angular 
movement  of  the  tower,  the  arc  L6(t)  can  be  approximated  by  the  measured  horizontal 
linear  displacement  Xm(t).  Substituting  Eqns.  4.2  - 4.5  into  Eqn.  4.*1  and  replacing 
LO  by  Xm  results  in 


T 12fc  L cos3[Xm{t)/L] 


Cm)Xn{t)  - 


d\ 


CDXm{t)\Xm(t)\ 


8 PwD L2  cos*[Xm{t)/L] 

+&<r.  -77 PwgD 2 - "y1m/n 

8 cos  2[Xm(t)/L] 


)sin  [Xm(t)/L)  = RL. 


(4.6) 


Since  the  measured  in-line  moment  Mm(t)  = RAri,  R — Mm(t)/Ar(.  Substituting 
this  expression  for  R into  Eqn  4.6  and  rearranging  results  in 
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Mn(t)  L/Art  + ^Xm(t)~ 

(£  D2  [Xm(t)/L]  = 

~12PwD  L cos3 [Xm (t)/L\ Cm  ~ 

8pw  DL2coS4[Xm(t)/L]CD  ^rn(t)\Xm(t)\.  (4.7) 


The  quantities  on  the  left  hand  side  in  Eqn.  4.7,  Mm(t),  Xm(t)  and  Xm(t)  are  known 
from  the  measured  values.  The  in-hne  motion  Xm(t)  was  monitored  by  the  LDT 
and  found  to  conform  with  the  input  signal,  asinu^f,  where, “a”  is  the  amplitude  of 
oscillation  and  ui^  the  oscillation  circular  frequency.  Thereupon,  the  velocity  Xm{t) 
and  acceleration  Xm(t)  were  computed  from  the  time  derivatives  of  Xm{t)  = a sin Udt. 

Equation  4.7  can  be  written  as 

fm(t)  = fi(t)Cm  + fo(t)CD,  (4.8) 


where 


fm(t) 

flit) 

fo(t) 


IQ  - ^ T 2 d\ 

-&r + t *-(t)  ~ (£ “kr*  ~ SM°  co.:put)/L] 


)sin[Xm(f)/L], 


/» 

~—pwD 


I2rw  L cos 3[Xm(t)/L] 


Xm[t)  and 


~o  Pv,D 


d\ 


8 L>cos4[Xm(t)/L] 


Xm(t)\Xm(t)\. 


Cd  and  Cm  are  the  unknown  quantities  in  Eqn.  4.8  and  can  be  obtained  by  minimizing 
the  squares  of  the  differences  between  the  computed  and  measured  values,  i.e.,  by 
minimizing  8 ? where 


1 N 

= J £[//(«.)  On  + MU)  Co  - MU) ?■  (4.9) 

N is  the  number  of  data  points  in  one  cycle.  From  this  equation  the  minimum  Cd 
and  Cm  were  obtained  by  solving  the  two  simultaneous  equations 
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dPx 

dCD 

dPx 

dCm 


0 and 
0. 


(4.10) 

(4.11) 


Equation  4.10  and  Eqn.  4.11  can  be  written  in  a matrix  form  as 


' El, 

E", Mt,)fD(t.) ' 

(Ml 

' E £, /.M/bM  ' 

. E,=i  U(U)1d{U) 

E.",  /?(<.) 

|cmj  I 

. E !,/.(<.)//(*.)  . 

(4.12) 


The  solution  of  Eqns.  4.12  results  in  the  following  equations  for  Cd  and  Cm 

r _ Ei/m(i.)/D(i,)i e /?(<,) - niMim i E[/p(i,)//(i,)i 

EWtoi  zm,)\  - e fo(t,  >/,(«.  w 1 ' 

r = af-wwi  e fhM  - )//(<■)]  ,,,,, 

E[/?W]E[/£W]-[E  *(*.)/,((.)?  • 1 ' 


where  the  summations  are  evaluated  from  * = 1 to  ./V  for  each  cycle.  To  put  more 
emphasis  on  the  large  data  values  and  thus  further  reduce  the  differences  between  the 
measured  and  predicted  forces  in  the  neighborhood  of  maximum  forces,  the  weighted 
least  squares  technique  was  also  applied  using  as  a weighing  factor.  The  suit- 

ability of  the  in-line  force  data  for  determining  Cp  and  Cm  was  also  evaluated  using 
Dean’s  (1976)  approach.  In  this  approach,  the  mean  square  error  fif  given  by  Eqn.  4.9 
defines  a quadratic  “error  surface”  which  is  a minimum  at  the  Cd  and  Cm  obtained 
by  solving  Eqn  4.12.  The  suitability  of  the  data  can  then  be  evaluated  from  the  lines 
of  constant  error  values,  S,(Cd , Cm),  which  are  ellipses.  For  example,  the  steeper  the 
slope  of  the  error  surface  with  Co,  the  better  the  data  are  suited  for  evaluating  the 
drag  coefficient,  see  Fig.  4.2. 


Cd  and  Cm  are  time  averages  over  several  cycles.  This  study  did  not  deal  with 
the  instantaneous  values  of  these  coefficients.  However,  since  the  measu rement  errors 
in  the  data  must  introduce  some  uncertainty,  the  variance  and  covariance  of  the 
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(cM> 


00  Data  Well-Conditioned  for  Deter- 
mining Drag  Coefficients 


(b)  Data  Well-Conditioned  for  Deter- 
mining Inertia  Coefficients 


Figure  4.2:  Contour  lines  defining  Error  surfaces  for  the  in-line  force  (Ref.  15). 


estimates  (Co  and  C,„)  were  computed.  The  variance,  aj,  in  the  value  of  any  function 
can  be  written  as 


(7 


2 

/ 


E(— )2- 

t 


(4.15) 


For  estimating  the  uncertainties  in  CD  and  Cm,  Eqn.  4.12  can  be  rewritten  as 


MM  = {/?}  (4.16) 

or 

2 

= (4.17) 

j=i 

where 

[a]  = 2 x 2matrix,  or 

akj  = Eh  X3(tt)Xk(tt),  k = 1, 2 and  j = 1,  2 and 

{/?}  = a vector  of  length  2,  or 

Pk  = Eh  fm(t,)Xk(t,)  k = 1,2. 

The  inverse  matrix  = [a]^1  is  closely  related  to  the  probable  (or  more  pre- 
cisely, the  standard)  uncertainties  of  the  estimated  parameters  a:  = (Co,Cm).  The 
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solution  to  Eqn.  4.12  (or  equivalently  Eqn.  4.16)  is 


N 


= EW  A = E c,„  E /»(«.)  **(<.)]• 

fc=i  fc=i  1=1 

Note  that  a3  corresponds  to  / and  /m(ft)  to  y,  in  Eqn.  4.15  thus 

d/  _ dQj  _ ^ v/i\ 

ay r SUt<)~h  * i(,)' 

Consequently,  the  variance  associated  with  the  estimate  parameters  a3  is 

= E E c*  [ E *<«]• 

k=\l=\  t=i 


(4.18) 


(4.19) 


(4.20) 


The  final  term  in  the  square  brackets  is  just  the  matrix  [a].  Since  this  is  the  matrix 
inverse  of  [C],  Eqn.  4.20  reduces  to 


<72(aj)  = C„. 


(4.21) 


In  other  words,  the  diagonal  elements  of  [ C ] are  the  variances  of  Crt  and  Cm,  while  the 
olf-diagonal  elements  C3k  are  the  covariances  between  Cp  and  Cm  (i.e.,  COV(Cd,  Cm )). 

An  estimate  of  the  goodness-of-fit  of  the  data  to  the  model  is  still  needed.  A 
simple  measure  of  the  goodness-of-fit  defined  as  the  average  percentage  error  between 
the  measured  and  predicted  signals  ( fm(t, ) and  was  used.  This  was  defined  as 


ERR  = 100 


\ 


nut.)  - /,(<.)]» 


i :fiM 

where  the  summation  is  evaluated  for  i = 1,2, ...,  N. 


% 


(4.22) 


Drag  Force 

Cq,  for  most  situations,  has  been  found  to  depend  on  Reynolds  number,  and  thus 
on  the  velocity.  In  an  attempt  to  minimize  this  dependence  and  maintain  the  drag 
coefficient  constant,  an  investigation  of  the  power  to  which  the  velocity  must  be  raised 
in  the  drag  force  equation  was  made.  The  drag  moment  in  Eqn.  4.3  was  replaced  by 
the  following  expression 

1 / rdi/cosO 

M-D  = — 2 PwDCd  J (r0)n  sgn(rd)  r dr, 


(4.23) 
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/ 

where  Cd  is  the  drag  coefficient  associated  with  this  expression  and  “ sgn ” is  a sign 
function  equal  to  ±1  depending  on  the  sign  of  the  argument,  rO.  Making  the  same 
assumption  Xm(t)  ~ LQ(t ) and  carrying  out  the  integration  of  Eqn.  4.23,  the  equation 
of  motion  (Eqn.  4.7)  reduces  to 


X„{t)  - 


(E^r.  “ 77 PwgD 


«=i 


p2 

12 Pw  Lcos*[Xm(t)/L\ 


cos  2[Xm(t)/L] 

Cm  Xm(t)  — 


)sin[Xm(f)/L]  = 


1 nr  dx 
2pw  lcos  [Xm(t)/L\ 


jn+2  j-  Cp  j | Xm(t)  |n 


5^n[Xm(<)]. 


(4.24) 


Ln  + 2J 1 L2 

In  this  equation,  the  quantities  on  the  left  hand  side,  Mm(t),Xm(t)  and  Xm(t)  are 
known  from  the  measurements.  Using  fm(t)  and  //(f)  defined  before,  and  defining 

fdr  = ~ pw  D sgn[Xm(t)], 


C{t)  = [ 

m = i 


dx 


cos  [Xn(t)/L\ 
Xm(t)t 


] and 


L2 


Equation  4.24  reduces  to 


/«(*)  = //(f)  Cm  + fdr  [C(t)]n+*  [D(t )]n  CD/(n  + 2).  (4.25) 

The  unknown  quantities  in  this  equation  are  Cm,  C/j  and  n which  can  also  be  obtained 
by  using  the  least  squares  equation: 

1 N 

6-  = N ELM'.)  + /*[C('.)]n+2  [»(',)]"  dD/(n  + 2)  - /m(f,)]2  (4.26) 

or  the  weighted  least  squares  equation 

? = JJ  E £('.)[/;('.)  Cm  + /*[C(f,)]"+2  [£>((,)]"  CW("  + 2)  - /m(i.)]2.  (4.27) 
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/ 

The  minimum  Cm,  Cp  and  n were  then  obtained  by  solving  the  coupled  nonlinear 
system  of  equations  given  by 


dPt 

dCm 

dsi 

dCD 

as? 

dn 


0, 


0 and 


0. 


(4.28) 

(4.29) 

(4.30) 


The  Newton-Raphson  method  for  nonlinear  systems  of  equations  was  applied  to  solve 
this  set  of  equations  (Eqns.  4.28  to  4.30).  The  method  is  given  in  many  text  books 
(e.g.  Press  et  al.,  1988). 

4.2  Transverse  Force 


For  several  reasons  the  transverse  force  has  been  the  subject  of  greatest  interest 
in  this  study.  First,  its  amplitude  can,  under  certain  conditions,  be  as  large  as  that  of 
the  in-line  force.  Second,  the  transverse  force  can  give  rise  to  fluid-induced  oscillations 
and  to  fatigue  failure.  Third,  even  small  transverse  motions  of  the  body  regularize  the 
wake  motion,  alter  the  span-wise  correlation  of  the  vortices,  and  drastically  change 
the  magnitude  of  both  the  in-line  and  transverse  forces. 

With  all  the  information  available  on  vortex-shedding,  there  is  no  simple,  explicit 
formula  to  predict  the  time  variation  of  the  vortex-induced  transverse  force.  This  is 
because  of  its  dependency  on  the  type  of  flow  and  structure  motion.  In  most  types  of 
flows  it  also  exhibits  some  degree  of  irregularity  that  is  usually  due  to  the  alternating 
eddies  behind  the  structure.  Therefore,  most  formulations  of  transverse  forces  are 
semi-empirical  and  based  heavily  on  experimental  data.  Of  course,  there  are  many 
numerical  models  available  (e.g.  Navier-Stokes  based-models,  discrete-vortex  models, 
wake  oscillator  models)  but  due  to  their  limitations  (such  as  laminar,  two-dimensional 
flows  at  low  Reynolds  number,  sensitivity  to  the  numerical  technique  used  and  their 
dependency  on  experimental  data  and  flow  visualization)  they  are  not  reviewed  here. 
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A number  of  authors,  including  Maull  and  Milliner  (1978),  Sarpkaya  and  ShoafF 
(1979)  and  Graham  (1980),  have  proposed  a method  that  uses  Blasius  equation  and 
the  “discrete- vortex  method” . Such  a method  assumes  that  the  induced-force  consists 
of  a component  due  to  the  attached  unsteady  irrotational  flow  past  the  body  and  a 
component  generated  by  the  vortices  shed  into  the  flow  a s a result  of  separation. 
The  application  of  this  model  requires  detailed  information  on  vortex  strengths  and 
trajectories.  Such  information  usually  is  difficult  to  obtain,  especially  for  the  cases 
involving  wave  flows. 

In  the  following  sections,  the  existing  semi-empirical  models  that  predict  trans- 
verse forces  are  briefly  discussed. 

4.2.1  Steady  Flow  Model 

This  model  is  similar  to  the  drag  force  in  Morrison’s  equation  and  is  written  in 
terms  of  the  lift  coefficient,  Cl  as 

Fl  = l-p„DtCLU2.  (4.31) 


The  model  was  originally  developed  for  steady  flow  where  it  has  been  confirmed  by 
many  investigators  that  the  transverse  force  spectra  has  a single  frequency,  predicted 
by  Strouhal  number,  St  = 0.2  for  103  < Re  < 2 x 10s.  Some  researchers  and  designers 
use  this  steady  flow  model  for  predicting  transverse  forces  for  oscillatory  and  wave 
flows.  In  order  to  account  for  the  time  variation  of  flow  velocity  and  transverse  force, 
a maximum  transverse  force  which  yields  a lift  coefficient  defined  by 


maximum  transverse  force 
jPtu  D i U^ax 

has  been  used. 

Others  use  an  rms  transverse  force  that  gives  a lift  coefficient  defined  by 


(4.32) 


C Lrms  — 


rms  value  of  the  transverse  force 

iPwDtUL*  ■ 


(4.33) 
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Other  methods  based  on  semi-peak-to-peak  values  of  the  transverse  force  or  different 
velocities,  such  as  those  corresponding  to  the  maximum  in-line  force  rather  than  to 
the  maximum  velocities,  will  not  be  discussed  here  in  order  to  avoid  confusion.  In  the 
author’s  opinion,  the  use  of  such  methods  for  oscillatory  flows,  especially  wave  flows, 
is  one  of  the  reasons  for  the  scatter  in  the  existing  Cl  data  as  illustrated  in  Fig.  2.14. 
This  is  because,  unlike  the  in-line  force  in  oscillatory  and  wave  flows,  the  transverse 
force  has  multiple  frequencies  and  exhibits  some  degree  of  irregularity  which  are  a 
result  of  the  different  shedding  mechanisms  in  these  types  of  flow. 

4.2.2  Quasi-Steady  Model 

Recently  the  steady  flow  model  has  been  extended  to  include  the  frequency  of 
vortex-shedding.  Verley  (1980),  followed  by  McConnell  and  Park  (1982b),  suggested 
the  following  simple  quasi-steady  model 

Fl  = -^Pw  DIClU 2 sin  uvt.  (4.34) 

The  model  is  based  on  the  instantaneous  values  of  flow  velocity  and  vortex-shedding 
frequency.  It  also  assumes  that  at  the  beginning  of  each  half  cycle  the  flow  starts 
from  rest  again  (i.e.,  does  not  account  for  the  flows  previous  history)  and  the  phase 
is  brought  to  either  0°  or  180°.  In  an  attempt  to  improve  this  model,  Bearman  et  al. 
(1984)  proposed  the  following  quasi-steady  model 

Fl  = ±-pw  D t Cl  cos  cj)  sin2  u ft,  (4.35) 

where  (j>  is  a function  of  time  given  by 

<p  = 0.2/C (7  [1  — cos  u)jt\  + ip,  (4.36) 

u)j  is  the  frequency  of  flow  oscillation  and  ip  is  a constant  to  be  adjusted  for  every 
half  cycle.  This  model  is  also  based  on  the  instantaneous  flow  velocity.  It  does  not 
predict  the  sign  of  the  force  which  depends  upon  the  sense  of  the  vortices  shed.  It 
is  based  on  the  assumption  that  the  Strouhal  number  is  constant  and  equal  to  0.2 
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(the  appropriate  value  for  steady  flow).  With  a suitable  choice  of  the  lift  coefficient 
Ci  for  each  half  cycle  of  the  incident  flow,  the  model  seems  to  work  reasonably  well 
for  K C > 20.  However,  the  variation  of  force  with  time  was  found  to  be  somewhat 
regular  and  the  vortices  tend  to  form  and  shed  in  a certain  prescribed  manner.  In 
general,  such  a model  may  be  helpful  in  understanding  the  flow  phenomenon  but  may 
not  be  suitable  for  design  purposes  as  it  is  sensitive  to  the  phase  angle  and  still  needs 
an  experimentally  determined  Cl. 

4.2.3  Series  Model 

Since  transverse  forces  exhibit  some  degree  of  irregularity,  Merrier  (1973)  found 
it  is  appropriate  to  express  the  transverse  force  in  the  form  of  a series.  Later  Isaacson 
(1974)  and  Chakrabarti  et  al.  (1976)  used  the  same  idea.  Their  results  showed  this 
model  to  be  superior  to  previous  models  with  regards  to  the  force-time  history.  In 
general,  the  series  form  of  the  transverse  force  can  be  written  as 

1 N 

FUt)  = Typw  DiUlZ  Cl(n)  cos(27rn/*  + <f>{  n)),  (4.37) 

Z n=l 

where  N is  the  number  of  harmonics,  Cx(n)  is  the  lift  coefficient  at  the  nth  harmonic 
and  <^(n)  is  the  phase  angle  associated  with  CL(n). 

In  1977  Sawaragi  et  al.  used  the  same  model,  but  introduced  the  assumption  that 
the  transverse  force  is  a random  variable.  Moreover,  the  spectrum  of  the  force  in  the 
region  of  the  dominant  harmonic  can  be  treated  as  a narrow-band  spectrum.  Based 
on  these  assumptions,  they  reduced  the  series  model  (Eqn.  4.37)  to 

pL(i)  = crL  \ cos(27T nft  + ^(n)),  (4.38) 

n=l  \ aL 

where  <r\  is  the  variance  of  the  transverse  force  defined  by 

= (4.39) 

and  S(nf)Af  is  the  spectral  energy  of  the  transverse  force  at  frequency  nf. 
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4.2.4  Proposed  Model 

Fixed  and  compliant  structures  in  deep  water  under  the  action  of  nonuniform 
oscillatory  or  wave  flows  have  varying  values  of  K C and  Re  starting  from  zero  at  the 
bottom  to  a maximum  at  the  surface.  Therefore,  using  the  maximum  velocity  in  the 
series  model  discussed  above  does  not  yield  a representative  total  transverse  force. 
Furthermore,  experimental  investigations  show  that,  in  general,  the  vortex-shedding 
process  for  this  type  of  flow  depends  on  KC  and  Re.  Thus,  the  lift  coefficient  is  not  a 
constant  over  the  structure.  In  an  attempt  to  improve  the  series  model  the  following 
assumptions  were  made.  The  structure  was  divided  into  N number  of  finite  elements 
with  length  A l.  Each  element  was  then  considered  to  be  subjected  to  uniform  flow 
(see  Fig.  4.3).  The  total  transverse  force  thus  could  be  expressed  by 

j M N 

Fl{1)  = ~PwD  J2'HCLmn  Ulmax  COs(2t TITlft  + <£(m)).  (4.40) 

* m=l n=l 

This  representation  allows  Cl  to  be  dependent  on  K C , Re  and  the  frequency  of  vortex- 
shedding. 


Figure  4.3:  Definition  sketch  showing  tower  elements  and  idealized  flow. 
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4.2.5  Fixed  Tower 


Transverse  force  data  measured  during  the  first  series  of  tests  were  reduced  us- 
ing the  proposed  model  presented  in  section  4.2.4.  Since  the  measured  data  was  in 
terms  of  the  moment  of  total  transverse  force  about  the  bottom  hinge  ( Mi(t )^),  the 


moment  form  of  Eqn.  4.40  was  used,  i.e. 


M N 


= £ CLmn  Unmax  r*  At  COs(2irm/f  + </’(„>), 


(4.41) 


m—1  n=l 

where  rn  is  the  distance  of  element  n from  the  bottom  hinge.  Substituting  Un  = 
anUd  sin  udt,  where  an  (=  arn/L ) is  the  amplitude  of  the  tower  oscillation  at  the 
element  n,  and  ud  (=  2irfd)  the  driving  frequency,  into  Eqn.  4.41  results  in 


L 1 a2u2  M N 

A 7t  = 2Pw  D ~ L ^ ^ ^ Clmn  r"  cos(<27rTnft  + 


(4.42) 


m=l  n=l 

Equation  4.42  is  analogous  to  the  Fourier  series  expression 

M M 

m =£c. eT"“  = Y.  Cm  cos(mwdt  + ^m)) 


m=l 


m=l 

M 


— ^(am  cos  TnuJdt  + bm  sin  mudt), 


(4.43) 


m=l 


where 

M = is  the  number  of  harmonics, 

f(t)  = ML(t)L/Aru 

= ^Em=l  f(t)  C°S  mUjf, 
bm  = sEm=i/W  sin  mudt, 

V’(rn)  = 4>(m)  = tan-1  bm/am  and 

Cm  = afil  + K,  = D (<J2  Ul/L2)  ZLl  CLmnrl  A 1. 

In  Eqn.  4.43  the  coefficients  Cm  were  obtained  by  taking  the  FFT  of  the  measured 
moment  of  transverse  force  Mi(t)L/ Art.  Then  knowing  Cm  (Fourier  components), 
the  minimum  lift  coefficients  at  each  element  along  the  tower  ( Cx„ , n = 1, 2, 3, ...,  TV.) 
were  obtained  by  using  the  method  of  least  squares: 

M 1 2,  ,2  N 

£ = £ {(nPn,  D Y‘ Clmn  rj  At)  - Cm]\ 

m=l  L L n 


(4.44) 
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where  the  minimum  Cimn  are  the  solutions  of 

m = and  n = l,2,...,N. 


a«-  = 0, 


(4.45) 


dCLmn 

In  other  words,  for  each  harmonic  the  minimum  lift  coefficients  along  the  tower  (Cl„) 
are  the  solutions  of  d62/dCLn  = 0 where,  n = 1, 2, ...,  N.  This  results  in  the  following 
N simultaneous  equations 

AX  = B,  (4.46) 


where 


A = 


•*6  ~,3  3 r3r3 

rl  rlr2  ' 1 r3 


r?r?  r6 


_3_3 


x = 


2' 1 '1  r2r3 


r3  »,3  r3  „3  r3  .,3 
' n' 1 rnr 2 rnr 3 

Clx 
CL2 


r3  3 
' 1 ' n 
r3  3 
'2 ' n 


■ and 


„3 


B = 


. 


The  solution  of  Eqn.  4.46  at  each  harmonic  is  then  given  by 


{X)  = {Cin}  = M]->{S}.  (4.47) 

4.2.6  Complaint  Tower 

When  the  articulated  tower  was  allowed  to  respond  freely  in  the  transverse  di- 
rection, the  transverse  force  was  not  measured.  Instead  the  transverse  motion  was 
measured.  Structure  motion  was  found  to  significantly  increase  the  transverse  force 
(see  McConnell  and  Park,  1982a).  One  way  to  quantify  the  effect  of  transverse  motion 
is  to  compute  the  change  in  the  transverse  force  due  t&  motion.  This  force  can  be 
determined  by  computing  the  force  needed  to  produce  the  measured  response.  This 
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method  was  used  in  this  study.  The  measured  data  for  the  tower  response  in  the 


transverse  direction  was  reduced  using  the  following  equation  of  motion 

AnVKO  = X/  Mt0tai  = Mg  + Md  + Mfy,  (4.48) 

where  Im , Mtotai,  Mg  and  Mr,  are  defined  in  Eqn.  4.1,  and  ip(t)  and  Mjy  are  the 
transverse  deflection  angle  and  applied  transverse  moment,  respectively.  For  the 
tower  shown  in  Fig.  4.4  the  moments  Mg,  Md  and  Mjy  were  defined  as 

J*  ^ rd\/  cos 

Mg  = (%2w,rt-  pw  Acgr  dr)  sin  rp(t)  (4.49) 

t=i  Jo 

Md  = -\pwD  *52  CD nrlAlri>2(t)sgn[ip(t)]  (4.50) 

^ n=l 

Y M N 

Mfy  = -pwD^2J2CLmnU^maxrnA£cos(2Trmft  + (f>(m))  (4.51) 

Z m=l  n=l 

where  Ac  = 7rZ?2/4  is  the  cross  sectional  area  of  the  tower  element.  The  added  mass 
moment  of  inertia  /'  is  given  by 

K=ElPu,D2Cmnr2nAL  (4.52) 

n=l  4 

Substituting  these  moments  into  Eqn.  4.48,  making  the  assumption  that  the  measured 
transverse  response  ym(t ) ~ Lip(t)  and  dx/ cos  ip(t)  ~ dx  gives 

(I°  + 'tjP»D*CmnrlAe)  ym(t)  + 

n=l  4 

1 N 

^Pv>DYl  CDnrn  A£yl(t)sgn[ym(t)]  + 

(? pwgD2dl  - J 2w'rt)ym(t ) = 

0 «=i 

1 a2u ;2  M N 

~2PwD~T~  ^ YjCLmnrlAicos(2'Kmft  + <f)(m)),  (4.53) 

m=l  n=l 

where  N is  the  number  of  elements  comprising  the  tower’s  underwater  portion,  M 
the  number  of  transverse  force  frequency  harmonics  and  k the  number  of  the  different 
components  making  up  the  tower. 
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Figure  4.4:  Definition  sketch  for  the  articulated  tower  showing  transverse  motion. 

The  left  hand  side  of  Eqn.  4.53  is  known  since  it  is  composed  of  the  measured 
quantities  ym(t),ilm(t)  and  ym(t).  The  quantities  ym(t ) and  ym(t)  were  obtained 
numerically  by  taking  the  first  and  second  derivatives  of  the  measured  response  ym(t). 
Equation  4.53  can  then  be  written  as 

1 a2a;2  M N 

/(*)  = ~Pw  D—J^  S ]£  CLmn  rjj  Sr  cos(27T mft  + (4.54) 

1 L m=ln=l 

where  f(t)  is  the  left  hand  side  of  Eqn.  4.53  but  of  opposite  sign.  Equation  4.54  is 
similar  to  Eqn.  4.40  thus  the  same  procedures  used  for  reducing  data  using  Eqn.  4.40 
were  applied  to  obtain  the  lift  coefficients  (Cl)  along  the  tower  for  each  harmonic. 


CHAPTER  5 

EXPERIMENTAL  DATA  ANALYSIS 


This  chapter  contains  the  data,  data  analysis  and  results  from  276  test  runs  for 
two  types  of  experiments.  The  sources  of  uncertainty  and  inaccuracy  associated  with 
the  data  are  also  discussed.  The  results  indicate  significant  effects  of  flow  nonuni- 
formity and  structure  transverse  motion  on  the  vortex-shedding  process,  loading  and 
thus  the  structure  response. 


5.1  In-line  Forces 

The  in-line  force  data  were  analyzed  using  Eqn.  4.8  from  which  the  inertia  and 
drag  coefficients  ( Cm  and  Co)  were  obtained  by  the  least  squares  curve  fit  technique. 
It  was  assumed  that  these  coefficients  are  constant  over  one  cycle.  The  mean  values 
of  Cm  and  Cq  were  then  computed  for  13  to  73  cycles  depending  on  the  driving 
frequency,  fa.  The  results  show  that  Cm  and  Co  do  not  vary  significantly  from  one 
cycle  to  the  next.  The  average  percentage  error  between  the  measured  and  predicted 
forces  was  also  computed.  For  most  of  the  results  the  error  for  Cm  and  Co  was  less 
than  5%.  A flow  chart  of  the  computer  program  (CDCM),  written  to  analyze  these 
in-line  force  data,  is  given  in  Appendix  C. 

For  the  range  of  Reynolds  numbers,  Re , and  Keulegan- Carpenter  numbers,  KC , 
tested,  the  Cm  and  Co  data  show  a clear  dependency  on  Re  when  KC  is  held  con- 
stant. This  is  consistent  with  the  results  in  the  lower  range  of  KC  for  a mechanically 
oscillated  cylinder  in  still  water,  see  Chakrabarti  et  al.  (1983).  The  results  for  Cm  and 
Co  versus  Re  for  different  values  of  KC  are  shown  in  Figs.  5.1  and  5.2.  These  figures 
show  that,  for  the  range  of  KC  and  Re  tested,  Cm  and  Co  are  strongly  dependent 
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on  A C and  Re  and  are  decreasing  as  Re  increases.  This  trend  is  also  consistent  with 
the  trend  reported  by  others  for  the  range  of  the  parameters  tested  (see  Sarpkaya, 
1976a).  Figure  5.1  also  suggests  that  the  fluid  added  mass  increases  as  the  cylinder 
amplitude  of  oscillation  increases,  i.e.,  as  KC  increases. 


Figure  5.1:  Cm  versus  Re  for  harmonically  oscillated  articulated  tower  (transversely 
unconstrained). 

It  was  observed  by  many  investigators  (e.g.  Sarpkaya,  1976a)  that,  for  lower  range 
of  KC,  the  in-line  force  is  usually  inertia  dominated.  The  data  from  these  tests  were 
also  more  suitable  for  extracting  Cm  than  Cd,  since  the  maximum  KC  attained  was 
9.35.  This  was  examined  by  calculating  the  “error  surfaces”  which  are  a minimum  at 
Comm  and  Cmmin • The  method  used  is  discussed  in  detail  by  Dean  (1976).  Figure  5.3 
shows  an  example  of  these  error  surface  contours  for  one  of  the  test  runs.  Most  of 
the  other  runs  show  the  same  shape  which  means  that  the  data  are  well- conditioned 
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Figure  5.2:  Cd  versus  Re  for  harmonically  oscillated  articulated  tower  (transversely 
unconstrained). 


for  determining  Cm.  The  scatter  in  Co  could  be  caused  by  the  errors  from  various 
sources  as  will  be  discussed  later  in  this  chapter.  In  summary,  the  number  of  useful 
data  points  for  establishing  the  Cp  — Re  relationship  is  small  compared  to  the  number 
used  for  Cm  especially  in  the  lower  range  of  KC.  The  phase  relationship  between  the 
measured  forces  and  the  corresponding  calculated  forces  was  also  examined.  The 
phase  differences  were  generally  found  to  be  small.  This  is  illustrated  in  Fig.  5.4, 
which  shows  a few  examples  of  the  force  time  history. 

In  an  effort  to  show  the  effect  of  flow  nonuniformity  on  these  data,  a comparison 
of  the  present  Cd  and  Cm  data  with  those  obtained  by  Sarpkaya  (1976a)  in  uniform 
oscillatory  flow  and  Bearman  et  al.  (1985)  in  waves  was  made.  The  results  are  shown 
in  Figs.  5.5  and  5.6.  Although  both  Cm  and  Co  show  the  same  trend  for  the  range 
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CD 


Figure  5.3:  Contour  lines  defining  error  surfaces  of  the  in-line  force  ( KC  = 8.4  and 
Re  = 8.4  x 104). 

of  KC  tested,  Cm  does  not  exhibit  a dependency  on  the  type  of  flow.  On  the  other 
hand,  Cd  seems  to  be  affected  by  the  flow  nonuniformity  in  that  its  magnitude  is 
less  for  nonuniform  and  wave  flows  than  for  uniform  oscillatory  flow.  The  Cd  values 
in  waves  and  nonuniform  oscillatory  flows  were  found  to  fall  within  the  same  range. 
The  differences  between  Cd  values  in  uniform  and  nonuniform  oscillatory  flows  could 
be  attributed  to  the  varying  Re  along  the  cylinder  in  the  nonuniform  oscillatory 
and  wave  flow  cases  since  Re  varies  from  a value  close  to  zero  at  the  bottom  to  a 
maximum  value  at  the  water  surface.  This  in  turn  could  cause  the  flow  to  change 
from  subcritical  to  critical  or  even  supercritical  (depending  on  the  water  depth  and 
the  cylinder  diameter)  along  the  cylinder.  It  may  also  cause  the  vortices  to  separate 
from  the  surface  at  different  moments  in  time  along  the  cylinder  leading  to  a phase 
difference  between  the  vortices  and  consequently  a smaller  correlation  length. 


in-line  force  (N.m)  in-line  force  (N.m) 
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Figure  5.4:  Phase  comparison  between  measured  and  predicted  in-line  forces  (a) 
h C = 6.6  and  Re  = 9.4  x 104  (b)  KC  = 8.4  and  Re  = 8.4  x 104  ( — measured, 

....  predicted). 
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Figure  5.5:  Comparison  of  Cm  for  different  flow  types. 
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Figure  5.6:  Comparison  of  Cp  for  different  flow  types. 
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Equation  4.25  was  also  used  instead  of  Eqn.  4.8  to  reduce  the  in-line  force  data 
and  to  examine  the  drag  force- velocity  relationship.  A flow  chart  of  the  computer 
program  (CDCMN)  to  compute  Cm,  Co  and  n is  given  in  Appendix  C.  Examples 
were  run  to  test  the  drag-velocity  relationship.  The  results  show  that  the  power  V 
to  which  the  velocity  is  raised  is  close  to  2 for  most  of  the  cases  tested;  see  Table  5.1 
for  samples  of  the  results. 


Table  5.1:  Samples  of  test  results  on  in-line  force  data 


Re 

I<C 

n 

CD 

Cm 

5.82  x 104 

5.9 

1.98 

0.35 

1.15 

8.35  x 104 

8.4 

1.81 

0.19 

1.16 

9.50  x 104 

5.8 

1.97 

0.55 

0.72 

8.30  x 104 

6.7 

1.86 

0.23 

0.90 

9.40  x 104 

6.6 

2.04 

0.13 

0.78 

1.20  x 10s 

6.6 

1.96 

0.15 

0.69 

4.20  x 104 

6.5 

1.88 

1.00 

1.77 

3.32  x 104 

6.5 

2.05 

0.88 

2.30 

5.2  Transverse  Forces 

In  this  section  the  analysis  of  the  data  pertaining  to  vortex-induced  transverse 
forces  is  presented  for  both  the  transversely  constrained  and  the  unconstrained  exper- 
iments. This  includes  the  analysis  of  measured  transverse  forces,  thermistor  signals 
and  transverse  motion. 

5.2.1  Constrained  Transverse  Motion 

When  the  tower  was  constrained  from  motion  in  the  transverse  direction,  the 
transverse  forces  were  analyzed  using  the  proposed  transverse  force  model,  Eqn.  4.42. 
The  lift  coefficients  Cl(„)  and  their  associated  phase  angles  for  the  first  eight 
harmonics  were  evaluated  using  Eqns.  4.47.  A flow  chart  of  the  computer  program 
(CLEF)  written  to  obtain  Ci(n)  and  <j)^  from  the  measured  transverse  force  data  is 
given  in  Appendix  C.  The  results  obtained  were  plotted  versus  different  parameters, 
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such  as,  KC , Re  and  (3.  The  values  of  the  lift  coefficients,  (?£(„)  were  found  to  be 
dependent  on  both  KC  and  Re.  However,  the  Cl(ti)  values  were  found  to  correlate 
better  when  plotted  versus  KC  for  constant  j3,  see  Figs.  5.7  and  5.8.  In  general,  Ci(n ) 
decreases  with  increasing  KC  as  shown  in  Figs.  5.7  and  5.8.  For  KC  > 7.0,  Cl( n) 
values  start  to  increase  slightly  indicating  the  possibility  of  a multi-peak  relationship 
with  KC.  This  behavior  has  also  been  observed  by  Bearman  et  al.  (1981)  and  Ikeda 
and  Yamamoto  (1981). 

The  dependency  of  the  CT(n)  values  on  (3  and  Re  was  also  found  to  be  strong, 
especially  for  the  first  four  harmonics  where  the  lift  coefficients  are  the  largest,  see 
Fig.  5.7.  Cnn)  were  also  plotted  versus  /3  for  constant  KC  as  shown  in  Figs.  5.9  and 
5.10.  The  trend  exhibited  in  Fig.  5.9  suggests  that  there  is  a critical  value  of  f3  or  Re 
at  which  the  lift  coefficient  is  a maximum.  Of  all  the  different  harmonics,  the  largest 
value  of  the  lift  coefficient  occurs  at  the  driving  frequency.  In  general,  the  magnitudes 
of  the  different  components  of  the  lift  coefficients  were  found  to  vary  between  0.0  and 
0.7  except  at  (3  = 2, 100  where,  Ci(\)  reached  a value  of  1.4. 

On  the  other  hand,  the  associated  phase  angles  show  no  orderly  dependency  on 
Re,  KC  or  /3.  As  can  be  seen  in  Figs.  5.11  and  5.12,  the  phase  angle  data  associated 
with  the  lift  coefficients  ((?£(„))  appear  to  be  random.  This  random  phase  could  be 
related  to  the  shedding  process  that  exhibits  modes  of  behavior  with  vortex  interaction 
between  newly  shed  and  earlier  formed  vorticies.  The  variability  in  the  modes  and 
thus  phasing  of  the  vorticies  shed  along  the  cylinder  could  also  be  a major  factor. 
The  phase  relationship  between  the  measured  transverse  forces  and  the  corresponding 
predicted  forces  using  the  proposed  model  for  computing  transverse  force  show  that 
while  the  maximum  forces  are  reasonably  correlated,  the  phase  differences  are  often 
large.  A few  examples  of  the  force  profiles  are  shown  in  Figs.  5.13  and  5.14. 
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Figure  5.7:  Lift  coefficient  harmonics  for  harmonically  oscillated  articulated  tower 
(transversely  constrained),  Cx(i)  to  Cx(4). 
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Figure  5.9:  Lift  coefficient  harmonics  versus  beta,  Cl(i)  and  Cu 2). 
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Figure  5.8:  Lift  coefficient  harmonics  for  harmonically  oscillated  articulated  tower 
(transversely  constrained),  Cl(5)  to  Cl(8). 
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Figure  5.10:  Lift  coefficient  harmonics  versus  beta,  Cx(3)  and  Cl^). 
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Figure  5.11:  Phase  angle  associated  with  the  lift  coefficients  harmonics,  (f>( i)  to  </>( 4). 
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Figure  5.12:  Phase  angle  associated  with  the  lift  coefficients  harmonics,  <^(5)  to  (f>(8). 
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measured  transverse  force  time  series 
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Figure  5.13:  Comparison  between  measured  and  predicted  transverse  forces  (Re  = 104 
and  KC  = 5.45). 
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power  spectrum  of  measured  transverse  force 
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Figure  5.14:  Comparison  between  measured  and  predicted  transverse  forces 

(Re  = 2 x 104  and  KC  = 5.9). 
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For  the  purpose  of  comparison  with  other  investigators’  data  the  transverse  force 
data  were  also  analyzed  in  terms  of  CLrms  which  is  defined  as 


C 


Lrms  — i 


Frm.it) 


PvDtUl, 


(5.1) 


where  FTma(t ) is  the  root-mean-square  ( rms ) value  of  the  measured  transverse  force. 
The  results  shown  in  Fig.  5.15  exhibit  the  same  trend  as  the  Cl  harmonics  and  are 
found  to  correlate  well  with  KC  for  constant  /?  with  magnitudes  varying  between  0.12 
and  0.6.  In  this  figure  the  results  show  that,  for  KC  < 5 (where  almost  no  published 
data  exist  for  nonuniform  oscillatory  flow)  Cxrma  increases  as  KC  decreases.  However, 
as  one  would  expect,  according  to  other  investigators  data  for  5 < KC  < 11,  CLrms 
increases  as  KC  increases  (see  Fig.  2.14). 


Figure  5.15:  RMS  lift  coefficient  versus  KC. 
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The  present  Cirms  data  were  also  plotted  with  other  investigators’  data  obtained 
for  waves,  see  Fig.  5.16.  As  shown  in  this  figure,  where  KC  overlaps,  the  lift  coeffi- 
cients, C Lrms , for  the  oscillating  tower  in  still  water  agrees  very  well  with  C/,rms  data 
for  waves.  For  the  range  of  Re  and  KC  tested,  Fig.  5.16  confirms  the  possibility  of  a 
secondary  peak  for  KC  between  2.0  and  3.0.  This  would  give  a.  KC  spacing  of  ~ 8.5 
from  the  main  peak  which  is  at  KC  ~ 11.0.  This  is  approximately  the  interval  of  KC 
found  by  Ikeda  and  Yamamoto  as  reported  by  Williamson  (1985).  Intervals  between 
the  Cl  peaks  were  found  to  be  approximately  7.5  for  KC  < 70.0. 
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Figure  5.16:  Comparison  of  Cirms  data  with  those  from  waves. 


Even  though  the  data  from  proprietary  studies  are  not  available,  the  ranges  of  Re 
and  KC  covered  by  these  data  sets  are  generally  known.  To  the  author’s  knowledge 
Cl  data  does  not  exist  for  most  of  the  range  of  the  parameters  covered  in  these  exper- 
iments. As  indicated  earlier  this  is  an  important  range  for  many  structural  elements 
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and  offshore  structures.  In  a paper  to  be  published  soon  by  Horton  et  al.  (1992)  will 
show  that  there  are  no  data  available  in  this  range,  see  Fig.  5.17.  This  figure  (taken 
from  Horton’s  to  be  published  paper)  is  reproduced  here  with  his  permisson. 


REYNOLDS  NUMBER 


Figure  5.17:  Range  of  data  from  wave  forces  hydrodynamic  experiments  (Ref.  27). 

There  are  several  ways  to  obtain  information  about  the  frequency  of  vortex- 
shedding.  One  common  way  is  to  perform  a spectral  analysis  on  measured  trans- 
verse force  data.  In  this  study  such  a method  was  used  and  the  assumption  that  the 
transverse  force  can  be  considered  a stationary  random  process  was  made.  Spectral 
analyses  were  carried  out  on  the  transverse  force  data  to  study  the  frequency  content 
of  the  signals.  This  was  done  by  using  a commercial  spectral  analysis  package  called 
GLOBALLAB  which  uses  an  FFT  algorithm.  Fractional  time  series  of  12.8  seconds 
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duration  and  0.025  seconds  time  intervals  (which  gives  a frequency  resolution  of  ap- 
proximately 0.078  Hz)  were  chosen.  The  final  spectrum  consisted  of  the  average  of 
5 spectral  estimates.  Since  the  frequency  of  vortex-shedding  in  oscillatory  flows  is 
known  (see  literature  review,  Chapter  2)  to  be  related  to  the  in-line  flow  oscillation 
frequency,  the  frequencies  in  the  power  spectra  were  normalized  by  the  driving  fre- 
quency. Also  to  illustrate  the  dependency  of  the  different  harmonics  in  the  transverse 
force  on  Re  and  AC,  the  power  spectra  were  plotted,  as  groups,  for  a range  of  KC 
tested  holding  Re  constant.  They  were  also  plotted  holding  0 constant.  Figure  5.18 
and  Fig  5.19  show  examples  of  these  spectra  for  constant  R^  and  0,  respectively. 

On  investigating  the  frequency  of  the  peaks  of  these  force  spectra,  it  was  ob- 
served that  the  frequency  of  the  dominant  peak  was  always  a multiple  of  the  driving 
frequency  (i.e.,  frequency  ratio  /„/ = N , where  N is  an  integer).  This  agrees  with 
the  findings  of  other  investigators  in  oscillatory  flows  and  waves  (e.g.  Chakrabarti 
et  a/.,  1976,  1983;  Bearman  and  Hall,  1987;  Justesen,  1989;  and  others).  The  plots 
in  Fig.  5.19,  where  0 is  the  correlating  parameter,  exhibit  the  same  trend  as  those 
of  other  investigators  shown  in  Fig.  2.9.  Because  the  frequency  of  vortex-shedding  is 
known  to  be  a very  important  parameter  in  this  type  of  fluid-structure  interaction 
problem,  the  measured  signals  from  the  thermistors  mounted  on  the  tower’s  surface 
were  also  analyzed  using  spectral  analysis  techniques.  The  spectra  of  the  thermistor 
signals  were  found  to  have  the  same  trend  with  KC,  when  0 was  held  constant,  as 
those  for  the  transverse  force  data  shown  in  Fig.  5.19.  Examples  of  thermistor  output 
spectra  are  shown  in  Figs.  5.20  and  5.21.  Each  figure  comprises  5 spectra  of  signals 
from  the  top  5 thermistors  (the  only  ones  that  worked  for  all  of  the  tests)  placed  one 
diameter  apart  along  the  length  of  the  cylinder,  see  Fig.  5.22. 

Similar  results  were  found  to  exist  for  all  values  of  0 tested.  This  supports  the 
finding  that  0 rather  than  Re  is  a better  correlating  parameter  for  the  vortex-shedding 
frequency.  In  general  these  spectra  show  that  the  frequency  of  the  dominant  peak  is 
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Figure  5.18:  Transverse  force  power  spectra  for  constant  Re(Re  = 1.15  x 104). 


Figure  5.19:  Transverse  force  power  spectra  for  constant  0(0  = 2,800). 


84 


Figure  5.20:  Thermistor  power  spectra  for  0 = 1,834. 


M 


FREQ.  RATIO,  f/fd 

Figure  5.21:  Thermistor  power  spectra  for  0 = 2, 161. 
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Figure  5.22:  Schematic  illustration  of  thermistors  locations. 

at  the  driving  frequency  for  lower  values  of  KC.  At  these  values  of  KC,  where  two 
asymmetric  vortices  remain  attached  to  the  cylinder,  the  ratio  fv/fd  = N = l.  As 
KC  increases  to  a certain  value  (depending  on  the  value  of  /3)  the  ratio  N shifts  to  2, 
where  two  vortices  with  opposite  signs  are  shed  in  one  cycle.  Figure  5.19  shows  that 
the  shift  to  N = 2 happens  for  a value  of  KC  between  7.1  and  7.9  for  (3  = 1,800. 
When  /?  is  2,800  this  shift  occurs  at  a lower  value  of  KC  ( KC  between  5.8  and  7). 
The  vortex  patterns  shown  in  Fig.  5.23  were  observed  in  these  tests.  This  observation 
helps  explain  the  behavior  shown  in  Fig.  5.19. 

Transverse  forces  are  not  always  well  correlated  over  the  length  of  the  cylinder. 
In  nonuniform  oscillatory  and  wave  flows,  forces  at  two  cross-sections  separated  by  a 
distance  as  little  as  two  diameters  may  be  in  phase  part  of  the  time  and  out  of  phase 
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KC  = 4 KC  = 9 

Figure  5.23.  Schematic  illustration  of  the  visualization  of  the  vortex-shedding  process. 

at  other  times.  Accordingly,  the  total  transverse  force  acting  on  a length  of  a few 
diameters  or  more  is  likely  to  be  considerably  smaller  than  that  suggested  by  a lift 
coefficient  derived  from  force  measurements  at  one  section.  The  vortex  correlation 
length  is,  therefore,  a very  important  parameter  in  interpreting  the  results  of  measured 
transverse  forces. 

In  this  study  correlation  lengths  were  computed  from  the  thermistor  signal  data. 
The  thermistors  were  mounted  at  10  different  locations  along  the  tower  model  (the 
top  7 thermistors  were  one  cylinder  diameter  apart  and  the  bottom  3 thermistors 
were  two  cylinder  diameters  apart).  The  correlation  was  calculated  between  the  top 
thermistor  signal  and  each  of  the  following  thermistors  using  the  correlation  coefficient 
7 xy(f),  defined  by 

,u  s„(f)-s„(/y  (5-2) 

where  Sxy(f)  is  the  cross  spectrum  between  the  two  signals  x(t)  and  y(t),  Sxx(f)  and 
Syyif)  are  the  auto  spectrums  of  the  signals  x[i)  and  y(t ),  respectively. 

Out  of  the  10  thermistors  on  each  side  of  the  tower,  only  the  top  5 thermistors 
were  found  to  work  during  all  tests.  They  are  located  up  to  6 diameters  from  the  water 


87 


surface.  The  computed  correlation  coefficients  of  these  thermistor  signals  were  plotted 
versus  the  span-wise  coherence  length,  z/D  (see  Fig.  5.24  for  definition)  for  all  values 
KC  and  (3  tested.  In  this  study,  the  thermistor  signals  were  considered  to  be  coherent 
if  their  correlation  coefficient  was  greater  0.9.  The  data  were  found  to  correlate  better 
with  (3  than  KC,  see  Fig.  5.25.  KC  has  been  found  to  correlate  coherence  length 
data  for  uniform  oscillatory  flows  (see,  Obasaju  et  al.,  1988  and  Kozakiewicz  et  al., 
1991).  Figure  5.25  indicates  that  the  vorticies  are  strongly  correlated  at  /?  = 2,161 
which  is  exactly  the  location  at  which  the  first  harmonic  lift  coefficient  peaks  (see 
Fig.  5.9).  This  phenomenon,  i.e.,  the  dependence  of  z/ D on  f3,  is  more  apparent  in 
a plot  of  z j D versus  /?,  as  shown  in  Fig  5.26.  The  general  behavior  exhibited  by 
this  figure  is  similar  to  that  reported  by  Overvik  (1982)  for  Scruton’s  (1967)  work  in 
steady  flow.  For  f3  = 2, 161,  is  near  unity  for  all  values  of  z/D  examined.  On 

the  other  hand  the  correlation  of  the  vortices  at  f3  = 3, 690  is  particularly  poor.  In 
general  the  results  presented  in  Fig.  5.26  suggest  that  the  vortex  correlation  length 
decreases  with  increasing  /?. 

D 


Figure  5.24:  Definition  sketch  of  span-wise  coherence  length. 
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Figure  5.25:  Span-wise  coherence  length  for  transversely  constrained  articulated 
tower. 


Figure  5.26:  Span-wise  coherence  length  for  transversely  constrained  articulated 
tower. 
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5.2.2  Unconstrained  Transverse  Motion 

When  the  tower  was  free  to  move  in  the  transverse  direction,  the  transverse  forces 
were  obtained  indirectly  by  measuring  the  transverse  response  and  computing  the 
forces  needed  to  produce  the  motion.  The  lift  coefficients  and  the  associated  phase 
angles  can  then  be  determined  in  a manner  similar  to  that  used  for  the  first  type 
of  experiment  (transversely  constrained  tower).  The  problem  with  this  approach  is 
that  reliable  values  for  Cm  and  Cd  for  the  range  of  KC  and  Re  (based  on  transverse 
motion)  experienced  are  needed.  Data  for  Cm  and  Cq  for  the  relatively  low  values  of 
KC  and  Re  in  these  experiments  are  limited  or  nonexistent.  The  values  of  KC  and 
Re  (based  on  transverse  motion)  for  the  tests  where  substantial  transverse  motion 
existed  are  shown  in  Table  5.2. 

To  quantify  the  effect  of  transverse  motion  on  lift  coefficients,  a computer  pro- 
gram (CLEUF)  using  Eqn.  4.54  was  written  for  computing  Cx(n)  and  the  associated 
phase  angles  from  the  measured  transverse  motion.  A flow  chart  of  this  program 
is  given  in  Appendix  C.  A few  cases  from  Table  5.2  were  tested  and  the  results  pre- 
sented in  Table  5.3.  For  these  cases  the  values  of  Cm  and  Co,  obtained  from  the 
measured  in-line  forces  (Figs.  5.1  and  5.2),  and  values  from  other  investigators  were 
extrapolated  to  the  transverse  motion  KC  values.  The  results  indicate  that  when 
transverse  motion  is  allowed,  the  lift  coefficients  are  increased  significantly,  at  least 
by  a factor  of  2.5.  McConnell  and  Park  (1982a)  oscillated  an  elastically-mounted 
cylinder  uniformly  in  still  water  and  found  that  response  in  the  transverse  direction 
increased  the  lift  coefficient  of  the  fixed  cylinder  by  a factor  of  up  to  2.0. 

The  Cl  data  contained  in  Table  5.3  were  then  used  to  evaluate  the  phase  relation- 
ship between  the  measured  and  predicted  transverse  motion.  This  requires  that  Eqn. 
4.53  be  solved  for  the  response  of  the  tower,  yv{t).  Because  the  equation  is  nonlinear, 
it  was  solved  numerically  using  the  Runge-Kutta  technique.  Details  of  this  technique 
are  given  in  Appendix  A and  a flow  chart  for  the  computer  program  (ATVSR)  used 
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to  calculate  the  tower  response  is  given  in  Appendix  C.  Results  of  the  program,  using 
the  Cl  data  given  in  Table  5.3,  show  that  fewer  harmonics,  than  for  the  constrained 
case,  are  needed  in  the  proposed  transverse  force  model  to  get  the  best  fit  between 
the  measured  and  predicted  response.  This  means  that  when  transverse  motion  is 
allowed  the  transverse  forces  due  to  vortex-shedding  are  more  regular.  Figure  5.27 
shows  a few  examples  in  which  only  2 harmonics  were  used.  The  lift  coefficients  were 
on  the  average  3 times  larger  than  those  corresponding  to  the  case  where  transverse 
motion  was  not  allowed. 


Table  5.2:  Test  conditions  for  cases  where  large  transverse  motion  were  measured. 


Ud 

in-line 

in-line 

in-line 

transverse 

transverse 

transverse 

(Hz) 

Re 

KC 

Re 

I<C 

0 

0.25 

2.72  x 104 

5.0 

5,432 

2.07  x 103 

4.12 

503 

0.25 

3.14  x 104 

5.8 

5,404 

2.28  x 103 

4.54 

503 

0.25 

2.07  x 104 

4.0 

5,170 

1.28  x 103 

2.68 

477 

0.25 

2.54  x 104 

4.9 

5,132 

2.28  x 103 

4.78 

477 

0.30 

3.19  x 104 

4.9 

6,501 

1.92  x 103 

3.18 

604 

0.30 

3.68  x 104 

5.7 

6,501 

2.49  x 103 

4.12 

604 

0.30 

4.21  x 104 

6.5 

6,500 

2.86  x 103 

4.74 

604 

0.30 

4.75  x 104 

7.3 

6,504 

3.24  x 103 

5.36 

604 

0.40 

5.20  x 104 

5.3 

9,890 

1.33  x 103 

1.44 

924 

0.40 

5.82  x 104 

5.9 

9,862 

1.90  x 103 

2.06 

923 

0.40 

6.75  x 104 

6.8 

9,923 

2.28  x 103 

2.47 

924 

0.40 

7.55  x 104 

7.6 

9,918 

3.04  x 103 

3.30 

922 

0.40 

8.35  x 104 

8.4 

9,926 

3.42  x 103 

3.71 

922 

0.40 

9.42  x 104 

9.3 

10,069 

3.80  x 103 

4.12 

922 

0.50 

7.22  x 104 

5.8 

12,424 

2.14  x 103 

1.86 

1,150 

0.50 

8.32  x 104 

6.7 

12,413 

2.38  x 103 

2.06 

1,154 

0.50 

9.23  x 104 

7.5 

12,396 

2.61  x 103 

2.27 

1,152 

0.50 

1.04  x 105 

8.4 

12,398 

2.85  x 103 

2.48 

1,150 

0.50 

1.13  x 10s 

9.1 

12,416 

3.56  x 103 

3.10 

1,150 

0.60 

8.04  x 104 

5.7 

14,159 

1.46  x 103 

1.11 

1,316 

0.60 

9.42  x 104 

6.6 

14,265 

1.95  x 103 

1.50 

1,302 

0.60 

1.04  x 10s 

7.3 

14,217 

2.33  x 103 

1.80 

1,296 

0.60 

1.18  x 10s 

8.2 

14,265 

2.71  x 103 

2.06 

1,317 

0.60 

1.28  x 10s 

9.0 

14,178 

3.26  x 103 

2.47 

1,318 

0.70 

1.07  x 10s 

6.6 

16,278 

2.49  x 103 

1.65 

1,512 
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Table  5.3:  Results  of  lift  coefficients  when  transverse  motion  existed. 


in-line 

Re 

in-line 

KC 

1/rmj/ D 
X100 

in-line 

Vr 

constrained 

O) 

unconstrained 

C't(i),(2)I(3) 

2.7  x 104 

5.0 

49% 

8.4 



0.57,  0.05,  0.11 

3.1  x 104 

5.8 

57% 

9.7 

_____ 

0.60,  0.05,  0.11 

4.7  x 104 

7.3 

60% 

14.6 

0.76,  0.06,  0.04 

5.2  x 104 

5.3 

16% 

14.0 

0.48,  0.024,  0.005 

5.8  x 104 

5.9 

28% 

15.6 



1.43,  0.053,  0.03 

6.8  x 104 

6.8 

35% 

18.1 



1.48,  0.05,  0.04 

7.2  x 104 

5.8 

21% 

19.4 

0.19,  0.05,  0.07 

0.68,  0.03,  0.007 

8.3  x 104 

6.7 

25% 

22.3 

0.25,  0.06,  0.15 

0.75,  0.03,  0.012 

9.2  x 104 

7.5 

28% 

24.8 

0.24,  0.12,  0.04 

0.77,  0.07,  0.002 

1.0  x 105 

6.2 

9% 

33.2 

0.31,  0.04,  0.03 

0.65,  0.03,  0.005 

1.0  x 105 

8.4 

33% 

27.8 

0.19,  0.20,  0.08 

0.61,  0.06,  0.02 

1.1  x 105 

6.6 

20% 

30.6 

0.10,  0.15,  0.07 

0.41,  0.01,  0.002 

1.3  x 10s 

9.0 

34% 

36.1 

0.14,  0.12,  0.03 

0.38,  0.02,  0.002 

Other  numerical  tests  were  also  performed  to  examine  the  effect  of  transverse 
motion  on  the  magnitude  of  transverse  forces.  In  these  tests  the  transverse  response 
was  predicted  using  the  measured  transverse  forces  on  the  constrained  tower  for  the 
same  in-line  conditions.  The  same  computer  program  “ATVSR”  was  used  to  predict 
the  tower  response  due  to  the  measured  transverse  forces.  This  was  carried  out  by 
replacing  the  expression  on  the  right  hand  side  of  Eqn.  4.53  by  the  the  measured 
transverse  force.  Many  cases  were  tested  and  most  show  that  the  amplitude  of  the 
predicted  response  are  less  than  that  measured  by  at  least  40%.  Results  also  show  that 
the  predicted  responses  are  more  irregular  than  those  measured.  Some  of  the  cases 
tested  are  shown  in  Figs.  5.28  to  5.31.  The  spectra  of  these  measured  and  predicted 
responses  (Figs.  5.29  and  5.31)  show  substantial  differences  in  energy.  These  results 
further  support  the  earlier  findings  that  the  transverse  motion  significantly  increases 
the  transverse  forces.  Therefore,  transverse  force  or  lift  coefficient  data  obtained  from 
measurements  on  stationary  (fixed)  structures  cannot  be  used  to  predict  the  response 
of  compliant  structures. 
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....  predicted). 
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Predicted  transverse  response 


Figure  5.28:  Comparison  of  measured  and  predicted  transverse  responses 

(Re  = 8.3  x 104  and  KC  = 6.7). 


Measured  spectrum 


Predicted  spectrum 


Figure  5.29:  Comparison  of  measured  and  predicted  spectra  of  transverse  responses 
(Re  = 8.3  x 104  and  KC  = 6.7). 
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Measured  transverse  response 


Predicted  transverse  response 


Figure  5.30:  Comparison  of  measured  and  predicted  transverse  responses  ( Re  = 105 
and  KC  = 8.4). 


Measured  spectrum 


Predicted  spectrum 


Figure  5.31:  Comparison  of  measured  and  predicted  spectra  of  transverse  responses 
(Re  = 10s  and  KC  = 8.4). 
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As  the  magnitude  of  the  transverse  force  was  found  to  be  dependent  on  the 
transverse  motion  and  in-line  flow  parameters,  the  frequency  of  this  force  is  like-wise 
believed  to  depend  on  these  quantities.  It  is  known  that  introducing  even  a small 
transverse  motion  will  cause  vortex-shedding  to  be  more  regular.  It  was  also  found 
that  when  the  tower  was  transversely  constrained,  the  vortex-shedding  frequency 
is  related  to  the  driving  frequency  and  0.  This  finding  is  believed  to  be  valid  for 
the  unconstrained  case  as  well.  The  frequency  of  vortex-shedding  was  analyzed  and 
presented  with  this  in  mind.  In  the  unconstrained  experiments  the  vortex-shedding 
frequency  was  obtained  from  spectral  analyses  of  the  measured  signals  from  the  ther- 
mistors mounted  along  the  cylinder.  These  spectra  show  vortices  being  shed  at  fre- 
quencies that  are  multiple  of  the  driving  frequency.  Out  of  the  several  peaks  in  the 
power  spectrum,  the  largest  peak  was  considered  to  be  the  dominant  frequency  of 
vortex-shedding. 

In  the  analysis  of  the  data  of  these  tests,  yTms/D  was  used  as  a parameter  to 
characterize  the  relative  magnitude  of  transverse  motion.  To  study  the  effect  of 
transverse  motion,  the  spectra  of  the  thermistor  signals  have  been  examined  in  many 
ways  i.e.,  for  constant  0 and  yrms/D , but  varying  KC  or  Re  or  VT.  It  was  found 
that  these  spectra  become  most  meaningful  when  grouped  together  for  constant  0 
and  yTma/D  but  varying  Vr.  The  test  runs  that  exhibited  transverse  motion  in  this 
second  set  of  experiments  covered  the  range  of  0 between  6,500  and  16,876.  For 
14,265  < 0 < 16,876,  (where  the  maximum  yTms/D  was  less  than  34%)  it  was  found 
that  the  dominant  vortex-shedding  frequencies  exhibited  the  same  trend  as  in  the 
first  set  of  experiments.  That  is,  at  small  values  of  VT  the  dominant  frequencies  of 
vortex-shedding  clustered  around  the  driving  frequency.  As  Vr  increased  to  a certain 
value  (depending  on  the  magnitude  of  yTTn3/D),  the  frequencies  clustered  around  the 
second  multiple  of  the  driving  frequency.  Figures.  5.32  to  5.36  show  a few  examples 
of  the  computed  spectra  for  the  thermistor  signals  at  different  0 and  yrms/D  values. 
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Vr  = 29.4 


Figure  5.32:  Thermistor  power  spectra  at  /?  = 14,265  and  yTms/D  = 29.5%. 


Figure  5.33:  Thermistor  power  spectra  at  0 = 14,265  and  yrm,/D  = 20.9%. 
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Figure  5.34:  Thermistor  power  spectra  at  = 14,265  and  yTms/D  = 8.2%. 
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Figure  5.35:  Thermistor  power  spectra  at  ft  = 16,260  and  yrmJD  = 16.8%. 
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Figure  5.36:  Thermistor  power  spectra  at  /?  = 16,260  and  yrms/D  = 8.8%. 

As  shown  in  Figs.  5.32  to  5.36,  the  transverse  motion  ( yrms/D ) does  affect  the 
conditions  under  which  the  frequency  of  vortex-shedding  shifts  from  the  first  multiple 
to  the  second  multiple  of  the  driving  frequency. 

To  show  the  effect  of  transverse  motion  on  the  transition  of  vortex-shedding  fre- 
quency, information  from  thermistor  spectra  was  plotted  (clustering  ranges  for  differ- 
ent transverse  motion  amplitudes  versus  reduced  velocity)  as  shown  in  Fig.  5.37.  This 
figure  shows  that  when  yrms/D  = 8.2%,  and  (3  = 14,265,  the  dominant  frequencies 
cluster  around  the  driving  frequency,  i.e.,  N = /„/ fd  = 1.  As  Vr  increases  to  17.5,  N 
shifts  to  2.  When  yrms/D  = 20.9%  , the  shift  to  N = 2 occurred  at  Vr  = 21.1.  As 
yrms/D  increases  further  to  29.5%  the  shift  to  N — 2 occurs  at  an  even  higher  value 
of  VT,  (Vr  = 26.4).  The  same  trend  was  also  observed  for  f3  = 16, 260  and  for  all  other 
values  of  (3  tested.  Consequently,  for  /?  > 14, 265,  the  smaller  the  amplitude  of  the 
transverse  motion  ( yrms/D ) the  smaller  the  value  of  Vr  at  which  the  shift  occurs. 


99 


Figure  5.37:  Schematic  of  the  effect  of  transverse  motion  on  the  behavior  of  vor- 
tex-shedding frequency. 

When  6, 500  < < 14, 265,  the  amplitudes  of  the  transverse  motion  were  larger 

than  when  ft  > 14,265  for  the  same  in-line  KC.  For  this  range  of  ft  it  was  found 
that  the  dominant  frequencies  of  vortex-shedding  did  not  cluster  around  the  driving 
frequency  at  small  values  of  Vr  as  they  did  for  ft  > 14,265.  Instead,  they  cluster 
around  a higher  multiple  of  the  driving  frequency,  the  value  of  which  depends  on  the 
magnitudes  of  the  transverse  motion  and  ft . For  example,  when  yTmal D was  56%  and 
ft  = 5,400,  N was  3 at  VT  = 4.6  and  remained  at  this  value  for  the  range  of  Vr  tested, 
see  Fig.  5.38.  On  the  other  hand  when  yrma/D  was  50%  and  ft  = 6,500,  N started 
at  3 at  Vr  = 6.6,  then  shifted  to  4 as  Vr  increased  to  7.7,  see  Fig.  5.39. 

Examples  of  thermistor  spectra  for  9, 900  < ft  < 14, 265  are  presented  in  Figs.  5.40 
to  5.44.  These  figures  show  that  N starts  at  4 then  shifts  to  higher  values  of  ./V  as  VT 
increases.  This  shift  depends  on  the  magnitude  of  the  transverse  motion  and  ft.  More 
data  are  needed  before  a complete  explanation  of  this  phenomenon  can  be  given. 
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Preferably  more  than  one  type  of  instrumentation  for  measuring  vortex-shedding 
frequency  should  be  used.  Perhaps  pressure  transducers  along  with  thermistors  could 
be  employed.  One  possible  explanation  of  the  phenomenon  described  above  might  be 
as  follows.  An  increase  in  Vr  can  be  thought  of  as  an  increase  in  KC  which  for  the 
articulated  tower  configuration  represents  an  increase  in  the  in-line  motion  amplitude. 
For  a transversely  constrained  cylinder,  the  number  of  vortices  shed  during  a cycle 
depends  on  the  in-line  motion  amplitude.  When  the  in-line  motion  is  symmetric,  as 
in  the  case  for  the  experiments  conducted  in  this  study,  the  increase  in  the  number 
of  vortices  shed  would  be  expected  to  occur  in  steps  of  two.  That  is,  as  the  in- 
line amplitude  increases,  one  additional  vortex  would  be  shed  per  half  cycle.  When 
transverse  motion  is  allowed,  the  vortex  shedding  process  is  altered  by  the  transverse 
motion.  In  general,  the  transverse  motion  is  not  regular  in  the  sense  that  many  half 
cycles  do  not  have  the  same  amplitude  as  their  adjacent  half  cycles.  This  could  result 
in  an  uneven  number  of  vortices  being  shed  during  a cycle  of  the  in-line  motion.  This 
might  explain  the  observed  incrementing  of  vortices  in  steps  of  one. 

Regardless  of  the  value  of  N at  which  the  dominant  vortex-shedding  frequencies 
cluster,  the  range  of  VT  over  which  these  frequencies  cluster  at  a particular  value  of 
N depends  on  the  magnitude  of  transverse  motion.  This  finding  is  similar  in  manner 
as  for  14,265  < 0 < 16,876.  For  example,  as  shown  in  Fig.  5.45  when  yrms/D  was 
16%,  N shifted  from  3 to  4 at  VT  — 9.95.  As  yrms/ D increased  to  28%,  N shifted  from 
4 to  5 at  VT  = 13.5,  then  to  6 at  Vr  = 14.  As  yTms/D  further  increased  to  35%,  the 
shift  from  N = 4 to  5 occurred  at  a higher  value  of  VT,  VT  — 15.6.  The  shift  to  N = 6 
then  occurred  at  VT  = 16.15,  which  is  higher  than  the  corresponding  VT  value  when 
yrms/D  = 28%.  At  yrms/D  = 40%  the  range  of  clustering  was  even  higher  and  so  on. 
This  general  behavior  was  exhibited  by  the  data  for  all  the  values  of  0 tested.  That 
is,  for  a constant  0 value,  the  smaller  the  transverse  motion  ( yrm3/D ),  the  earlier  the 
shift  of  vortex-shedding  frequency  to  a higher  multiple  with  respect  to  VT. 
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Figure  5.38:  Thermistor  power  spectra  at  /?  = 5,400  and  yrms/ D = 56%. 


Figure  5.39:  Thermistor  power  spectra  at  /?  = 6,500  and  yrms/ D = 50%. 
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Figure  5.40:  Thermistor  power  spectra  at  /?  = 9,900  and  yrma/D  = 16%. 


Figure  5.41:  Thermistor  power  spectra  at  (3  = 9,900  and  yrm3/ D = 28%. 
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Figure  5.42:  Thermistor  power  spectra  at  /3  = 9,900  and  t/rma/ D = 35%. 
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Figure  5.43:  Thermistor  power  spectra  at  /?  = 9,900  and  yTma/ D = 40%. 
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Vr  = 24.0 


Figure  5.44:  Thermistor  power  spectra  at  /?  = 12,500  and  yrms/D  = 33.5%. 
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Figure  5.45:  Schematic  of  vortex-shedding  frequency  behavior  (/ 3 = 9,900). 
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It  must  be  pointed  out  that  for  6,500  < /?  < 14,265,  when  yTms/D  exceeded 
35%,  the  dominant  vortex-shedding  frequencies  were  found  to  cluster  at  a value  of 
N = 4 for  a larger  range  of  Vr.  This  could  reflect  an  improved  correlation  of  the 
vortex-shedding  along  the  cylinder  due  to  the  larger  transverse  motion  which  causes 
the  vortices  to  be  more  regular.  To  gain  more  insight  into  this  clustering  behavior, 
the  correlation  coefficients  7xy(/)  (defined  in  Eqn.  5.2)  between  the  top  thermistor 
signal  and  each  of  the  lower  thermistors’  signals  were  calculated.  These  correlation 
coefficients  were  then  plotted  against  the  vortex  span-wise  coherence  length,  z/D,  as 
shown  in  Fig.  5.46  for  all  (3  values.  As  one  can  see  from  this  figure,  the  vortices  were 
highly  correlated  for  all  /?  values,  but  especially  for  /3  = 9, 900,  where  the  correlation 
coefficient  was  uniform  over  a longer  length  (z/D  > 4).  This  can  explain  why  the 
frequency  of  vortex-shedding  clustered  at  a value  of  TV  = 4 for  a larger  range  of  VT 
when  /3  = 9,900.  Compared  to  Fig.  5.25,  where  transverse  motion  was  not  allowed, 
Fig  5.46  shows  that  the  vortices  are  strongly  correlated  over  a wider  range  of  ft  and 
the  span-wise  coherence  length  is  larger  when  transverse  motion  is  allowed. 


Figure  5.46:  Span-wise  coherence  length  for  transversely  unconstrained  articulated 
tower. 
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The  frequency  of  tower  response,  /r,  was  obtained  from  the  frequency  analysis 
of  the  tower’s  transverse  motion.  This  frequency  (/r)  was  defined  as  the  dominant 
peak  in  the  response  spectrum.  The  frequencies  of  these  spectra  were  normalized  by 
both  the  natural  frequency  of  the  system,  /n,  and  the  driving  frequency,  fd.  The 
driving  frequency  was  found  to  be  the  more  meaningful  frequency  with  which  to 
normalize.  These  spectra  were  then  plotted  as  groups  with  constant  ft  but  varying 
VT.  Figures.  5.47  to  5.51  show  examples  of  these  computed  spectra  for  6,500  < 
ft  < 14,265.  As  shown  in  these  figures,  for  the  range  of  ft  values  examined,  the 
ratio  /r/ fd  (=  M)  always  has  a value  of  1.  The  vortex-shedding  frequency,  /„,  was 
also  higher  than  the  tower’s  response  frequency,  /r,  and  the  natural  frequency  in  still 
water,  /„.  At  the  lower  values  of  ft,  the  tower  response  exhibited  very  little  energy  at 
its  natural  frequency.  As  ft  increased,  the  energy  at  the  natural  frequency  increased 
until  it  reached  its  maximum  value  at  ft  = 12,415  and  Vr  = 30.3  then  decreased 
again.  Under  these  conditions  the  energy  at  the  natural  frequency  was  about  75% 
of  the  energy  at  fr/fd  = 1.0.  Such  behavior  can  be  interpreted  as  follows.  It  is 
believed  that  very  little  energy  is  exhibited  in  the  response  spectra  at  the  system’s 
natural  frequency  for  the  following  reason.  At  such  a low  value  of  /„  there  was  no 
vortex-shedding  observed.  It  was  also  found  from  the  first  set  of  tests  (constrained 
transverse  motion)  that  as  ft  increases,  the  transverse  force  increases  until  it  becomes 
a maximum  at  a certain  value  of  ft  called  ftcr  then  decreases  again,  see  Fig.  5.9.  By 
analogy  ft  = 12,415  is  ft„  for  the  case  when  transverse  motion  was  allowed.  Thus, 
this  could  explain  the  behavior  of  the  energy  exhibited  in  the  response  spectra  at 
the  natural  frequency  for  ft  < 12,415,  where  the  transverse  force  is  very  small.  At 
ft  — 12,415,  the  transverse  force  is  a maximum  and  therefore  the  energy  in  the 
response  spectra  is  a maximum.  Then,  as  ft  increases,  the  transverse  force  starts  to 
decrease  causing  less  energy  in  the  response  spectra. 
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Figure  5.47:  Transverse  motion  power  spectra  for  0 = 5,400. 


Figure  5.48:  Transverse  motion  power  spectra  for  /?  = 6,500. 
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Figure  5.49:  Transverse  motion  power  spectra  for  /?  = 9, 900. 


Figure  5.50:  Transverse  motion  power  spectra  for  /?  = 12,500. 
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Figure  5.51:  Transverse  motion  power  spectra  for  j3  = 14,265. 

The  response  of  a compliant  structure  is  very  much  dependent  on  its  natural 
frequency  (i.e.,  the  natural  frequency  of  the  fluid-structure  system).  Thus,  at  least 
indirectly,  transverse  forces  are  dependent  on  the  natural  frequency  of  the  system. 
There  was,  however,  no  attempt  made  to  examine  these  effects  in  this  study,  i.e.,  the 
system  still  water  natural  frequency  (as  defined  below)  remained  unchanged  during 
the  course  of  experiments. 

In  this  study,  the  system  natural  frequency  was  obtained  by  measuring  the  tower 
free  oscillation  in  still  water  (see  Fig.  5.52).  The  measured  value  compared  well  with 
the  analytical  expression  given  by 

= (5-3) 

where  k3  is  restoring  moment  coefficient  due  to  weight  and  buoyancy  of  the  tower, 
I0  and  /'  are  the  mass  and  added  mass  moment  of  inertia  about  the  bottom  hinge, 
respectively. 
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tower  in-line  free  oscillation 


tower  transverse  free  oscillation 


The  natural  frequency  of  the  system  appeared  to  change  with  the  operating  con- 
ditions of  the  experiments.  The  reason  for  this  apparent  change  was  attributed  to 
changes  in  the  added  mass.  The  “still  water”  natural  frequency  defined  above  was, 
however,  used  in  analyzing  the  response  and  thermistor  data. 

It  is  believed  that  there  should  be  a relationship  between  /„  and  the  other  frequen- 
cies involved  (/r,  /j  and  /n).  The  ratios  of  these  frequencies  were  carefully  examined 
using  the  information  from  the  thermistor  and  response  spectra  at  all  values  of  /? 
tested.  It  was  found  that  at  /3  = 5,400  and  t/rm4/D  = 56%  the  ratios  N and  M are 
integers,  more  specifically  N = 3 and  M = 1,  see  Figs.  5.38  and  5.47.  Under  these 
conditions  the  ratio  /„//„  is  also  found  to  be  an  integer,  fv/ fn  = 5. 
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Similar  situations  were  also  found  to  occur  at  /?  = 9,900  and  yrrna/ D = 16%  and 
at  /?  = 14,265  and  yrms/D  = 8.2%.  At  0 = 9,900,  it  was  found  that  N = 4,  M = 1 
and  fv/fn  = 10,  while  at  (3  = 14,265,  N = l,  M = \ and  /„/ fn  — 4. 

It  is  known  that  in  steady  flow,  when  the  natural  frequency  of  the  system  is 
considerably  greater  than  the  frequency  of  vortex-shedding  for  a stationary  cylinder, 
the  cylinder  has  a very  small  response  at  the  natural  frequency.  As  the  flow  velocity 
is  slowly  increased,  the  frequency  of  vortex-shedding  increases  as  does  the  cylinder 
amplitude  at  the  natural  frequency.  This  process  continues  until  the  amplitude  takes 
control  of  the  vortex-shedding  frequency  resulting  in  a maximum  amplitude  of  oscil- 
lation. A further  increase  in  the  flow  velocity  causes  the  amplitude  to  decrease  until 
eventually  control  of  the  shedding  frequency  by  the  transverse  motion  is  lost.  When 
the  amplitude  of  the  transverse  motion  is  large  enough  to  change  the  frequency  of 
vortex-shedding  to  that  of  the  response,  “lock-in”  occurs.  This  will  happen  when 
the  flow  speed  is  such  that  the  stationary  cylinder  vortex-shedding  frequency  is  in 
the  neighborhood  of  the  system  natural  frequency.  In  periodic  flows  the  concept 
of  lock-in  also  occurs.  However,  since  there  are  usually  more  than  a single  station- 
ary cylinder  vortex-shedding  frequency,  the  process  is  more  complex.  McConnell 
and  Park  (1982a,b)  found  that  lock-in  for  uniform  oscillatory  flow  occurs  when  the 
vortex-shedding  frequency  and  the  system’s  natural  frequency  are  related  to  the  fluid 
oscillating  frequency  in  an  integer  manner.  In  their  study  this  occurred  when  the 
fr/ fn  = fv/  fn  — 1 and  fv/ fd  = 6.22  for  4.4  < Vr  < 6.6.  This  means  that  the  natural 
frequency  of  their  system  was  higher  than  the  driving  frequency.  In  the  experiments 
conducted  in  this  study,  the  natural  frequency  of  the  system  was  much  lower  than 
the  driving  frequency.  This  is  usually  the  case  for  offshore  structures. 

To  examine  the  lock-in  condition  for  the  experiments  of  this  study,  all  combina- 
tions of  the  ratios  of  the  frequencies  involved  were  carefully  evaluated  for  the  cases 
where  considerable  transverse  motion  was  excited  at  the  natural  frequency  of  the 
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system.  In  this  study,  this  was  found  to  occur  at  j3  = 14,265  and  lock-in  was  ob- 
served when  fr/fn  = fv/fn  = 4 and  fv / fd  = l.  This  occured  over  a wide  range  of 
Vr  (22  < Vr  < 36),  see  Fig.  5.51  and  Figs  5.32  to  5.34.  In  other  words,  the  vortex- 
shedding  frequency,  /„,  changed  to  match  the  frequency  of  the  tower  response,  fr, 
which  was  one  of  the  multiples  of  the  system’s  natural  frequency,  /„.  Thus,  in  order 
to  predict  lock-in,  one  should  know  the  stationary  cylinder  vortex-shedding  frequency, 
system  natural  frequency,  driving  frequency,  reduced  velocity  and  the  frequency  pa- 
rameter. 

Transverse  motion  data  obtained  in  this  study  were  plotted  as  t/rm*/ D and  2/1/3/D 
versus  fd / fn  for  constant  KC.  The  results  presented  in  Fig.  5.53  show  that  the 
maximum  transverse  motion  amplitude  always  occurs  at  a ratio  of  fd / fn  = 1.67  for 
all  KC  values.  This  means  that  the  maximum  transverse  motion  occurred  at  a driving 
frequency  fd  = 0.25  Hz.  At  this  frequency,  the  ratios  fr/fd,  fv/fn , fv/ fd,  and  /„/ fr 
were  always  an  integer  for  all  KC  values  tested,  in  particular,  frf fd—  1>  fv/fn  = 5 
and  fjfd  = fv/fT  = 3. 

Knowing  the  range  of  the  governing  parameters  at  which  transverse  excitation 
is  likely  to  occur  is  vital  in  both  design  and  operation  of  offshore  structures.  The 
classification  societie’s  rules  for  the  design  and  construction  of  offshore  structures 
(such  as  Det  Norske  Veritas,  DNV,  1982),  recommend  the  use  of  VT  as  a parameter  to 
characterize  this  problem.  For  example,  for  steady  flow,  DNV  rules  recommend  the 
consideration  of  vortex-shedding  excitation  when  VT  > 3.5.  For  oscillatory  flows  and 
waves,  the  rules  recommend  that  cross-flow  excitations  be  considered  when  Vr  > 1.0. 
The  larger  range  of  VT  designated  for  oscillatory  flows  indicates  the  uncertainties  in 
this  phenomenon  and/or  the  lack  of  information  for  this  type  of  flow.  To  gain  more 
insight  into  this  problem,  the  tower  response  due  to  vortex-shedding  was  examined 
with  respect  to  the  parameters  controlling  the  excitation.  The  root-mean-square  and 
the  average  of  the  highest  1/3  values  of  the  transverse  motion  amplitudes  ( yrma  and 
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y\/z)i  normalized  by  the  tower’s  diameter  D , were  plotted  in  various  ways  such  as 
Urma/D  and  yi/3/D  versus  KC  while  holding  either  /?e,/?,  or  Vr  constant.  Plots  of 
yrma/D  and  2/1/3//}  versus  Vr,  Re  and  0 while  holding  KC  constant  were  also  made. 
It  was  found  that  plots  versus  VT  with  constant  KC , as  shown  in  Figs.  5.54  and  5.55, 
were  the  most  meaningful.  As  noted  from  these  figures,  for  any  KC  the  amplitude 
of  oscillation  increases  gradually  as  Vr  increases  and  reaches  a maximum  when  /„//„ 
approaches  an  integer  then  decreases  again  as  Vr  increases  further.  The  maximum 
rms  amplitude  of  oscillation  ( yrm, ) in  the  transverse  direction  due  to  vortex-shedding 
at  the  lowest  KC  {KC  = 4.2)  was  found  to  be  about  30%  of  the  cylinder’s  diameter. 
The  corresponding  1/1/3  was  about  50%/}.  The  results  shown  in  Figs.  5.54  and  5.55 
agree  with  the  results  of  Sumer  and  Fredspe  (1988,  1989)  for  uniform  oscillatory  flow. 
The  figures  also  demonstrate  that,  the  larger  the  KC , the  larger  the  value  of  VT  at 
which  the  maximum  transverse  motion  occurs. 


Figure  5.53:  RMS  transverse  motion  versus  frequency  ratio. 
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During  the  course  of  the  experiments  it  was  observed  that,  as  KC  increased,  the 
amplitude  of  transverse  oscillation  increased  until  the  cylinder  hit  the  ends  of  the 
apparatus.  This  happened  at  KC  = 5.8.  At  this  point  yTms  was  about  60%D  and 
the  corresponding  yi/3  was  about  90 %D.  Beyond  KC  = 5.8  and  over  a certain  range 
of  VT,  depending  on  the  value  of  KC,  the  amplitude  of  transverse  motion  was  larger 
than  that  allowed  by  the  apparatus.  The  dotted  lines  in  these  figures,  therefore,  show 
extrapolated  values  of  the  response  that  would  have  occurred  if  a larger  displacement 
had  been  allowed.  Thus,  the  response  amplitude,  shown  in  Fig.  5.56,  will  reach  a 
maximum  then  decrease  displaying  a peak  similar  to  that  obtained  by  Bearman  and 
Hall  (1987)  for  uniform  oscillatory  flow. 


Figure  5.54:  RMS  transverse  response  versus  reduced  velocity. 
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Figure  5.55:  Average  highest  1/3  transverse  response  versus  reduced  velocity. 
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Figure  5.56:  RMS  transverse  response  versus  Kuelegan- Carpenter  number. 
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5.3  Sources  of  Uncertainty  and  Inaccuracy 

The  purpose  of  the  experiments  carried  out  in  this  study  was  to  quantify  the 
effect  of  flow  nonuniformity  and  transverse  motion  on  vortex-induced  transverse  loads 
and  response.  Although  most  of  the  data  obtained  from  the  experiments  were  well 
behaved,  some  displayed  scatter.  For  example,  the  drag  coefficient  ( Cd ) exhibits 
some  scatter  at  lower  Re  and  KC,  particularly  for  Re  around  2 x 104  and  KC  < 6. 
There  is  also  some  scatter  in  the  lift  coefficient  data  for  1 , 300  < ft  < 1, 800  and 

0.4  x 104  <Re<  1.6  x 104. 

From  the  physical  insight  gained  from  the  results  presented  in  this  chapter,  it  is 
believed  that,  the  variation  of  the  velocity  and  thus  Re  along  the  tower  (that  varies 
form  zero  at  the  bottom  to  some  maximum  at  water  surface)  is  the  key  to  the  scatter 
in  some  of  the  data.  Some  of  the  possible  consequences  of  this  nonuniformity  are 
listed  below: 

1.  The  vortices  can  separate  from  the  surface  at  different  times  along  the  tower. 
This  results  in  phase  differences  between  the  vortices  shed  which  would  lead  to 
a smaller  correlation  length. 

2.  Different  flow  regimes  (subcritical,  critical  or  even  the  postcritical)  can  exist 
along  the  tower  at  any  given  time.  The  transition  from  one  regime  to  another 
is  not  only  dependent  on  Re  but  also  on  the  conditions  of  the  cylinder  surface 
and  motion  as  well.  In  these  experiments,  there  were  some  changes  in  cylinder 
surface  roughness  over  the  duration  of  the  experiments  (due  to  the  exposure  of 
the  aluminum  cylinder  to  the  chlorinated  water  for  a period  of  approximately 
three  months).  Small  vibrations  were  also  present  at  the  lower  frequencies  and 
amplitudes. 

The  scatter  observed  in  some  of  the  data  can  then  be  interpreted  as  resulting  from 
any  or  all  of  these  factors  listed  above. 
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Other  possible  sources  of  inaccuracies  in  the  data  include: 

1.  Errors  introduced  by  the  12-bit  A/D  convertor  which  only  had  a resolution  of 
0.0025  volts.  This  results  in  a resolution  of  ±1.9  mm  for  the  in-line  motion 
and  ±1.0  mm  for  the  transverse  motion  measured.  In  terms  of  the  moments 
measured,  this  is  about  ±0.004  N.m  for  the  in-line  and  ±0.005  N.m  for  the 
transverse  moment. 

2.  The  input  function  for  the  in-line  drive  mechanism  was  x(t)  = asinu;^-  This 
function  rather  than  the  measured  motion  was  used  in  the  analysis.  Measure- 
ments of  the  actual  motion  were  found  to  be  slightly  different  than  this  function. 
These  differences  resulted  from  inaccuracies  in  the  control  system.  The  errors 
introduced  by  this  procedure  are  estimated  to  be  less  than  2%. 

3.  In  order  to  compute  the  tower’s  transverse  velocity  and  acceleration,  the  re- 
sults of  measured  transverse  motion  were  numerically  differentiated.  After  each 
differentiation  the  results  had  to  be  filtered  to  remove  high  frequency  noise  in- 
troduced by  the  numerical  procedure.  Errors  were  introduced  in  the  results  by 
the  differentiation  and  filtering  techniques.  Higher  order  (4th  to  10th)  lowpass 
digital  filters  were  used  to  minimize  filter  errors  which  are  estimated  to  be  less 
than  6%. 


CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 


This  chapter  summarizes  the  most  important  observations  made  during  the  present 
investigation,  draws  conclusions  from  the  findings  and  makes  recommendations  for  fu- 
ture work. 

6.1  Summary  of  the  Results 

This  study  revealed  several  interesting  facts  about  the  effect  of  flow  nonunifor- 
mity on  fluid-induced  loading  on  fixed  and  compliant  structures  for  6.1  x 103  < Re  < 
1.3  x 105,  2.5  < KC  < 9.5  and  3.1  < VT  < 46.  Regarding  in-line  forces  on  compli- 
ant structures,  flow  nonuniformity  was  found  to  have  a greater  influence  on  the  drag 
coefficient  than  on  the  added  mass  coefficients.  The  drag  coefficient  was  less  than 
that  for  the  equivalent  uniform  oscillatory  flow.  Sarpkaya  (1976a)  and  others  showed 
that  the  vortex-induced  forces  in  the  in-line  direction  are  a major  factor  influencing 
the  drag  force.  The  effect  of  flow  nonuniformity  on  the  vortex-shedding  process  prob- 
ably accounts  for  the  difference  in  the  drag  coefficients.  However,  both  inertia  and 
drag  coefficients  show  clear  dependency  on  Re  and  KC.  Inertia  was  the  dominant 
component  of  the  force  for  the  range  of  Re  and  KC  attained  in  these  experiments. 

The  transverse  force  data  for  the  fixed  tower  identified  a number  of  differences 
in  the  vortex-shedding  behavior  between  uniform  and  nonuniform  oscillatory  flow. 
One  difference  is  nonuniform  oscillatory  flow  produces  multiple  frequencies  of  vortex- 
shedding.  This  is  most  likely  due  to  the  velocity  variation  along  the  tower  which 
can  cause  the  vortices  to  separate  from  the  surface  at  different  points  in  time.  The 
dominant  frequencies  cluster  around  the  driving  frequency  at  lower  values  of  KC.  As 
KC  increases  to  a certain  value  (depending  on  the  value  of  /?)  these  frequencies  shift 
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to  a value  twice  the  driving  frequency.  An  explanation  of  this  behavior  was  given  in 
the  previous  chapter.  Flow  nonuniformity  was  also  found  to  affect  the  magnitude  of 
the  lift  coefficient  and  the  span-wise  coherence  length.  Both  exhibited  a maximum  at 
a certain  value  of  (3  referred  to  here  as  (3  critical  or  (3cr.  In  these  experiments  (3^  was 
found  to  be  2, 160.  The  magnitude  of  the  lift  coefficient  was  found  to  depend  on  both 
KC  and  Re  rather  than  just  Re  as  it  is  for  steady  flow,  or  just  KC  as  it  is  for  uniform 
oscillatory  flow.  For  the  range  of  KC  tested,  the  lift  coefficient  data  of  these  tests 
showed  good  agreement  with  data  of  other  investigators  for  wave  flows.  The  data 
suggest  that  the  lift  coefficient  might  exhibit  a secondary  peak  around  KC  ~ 2.5. 

The  model  proposed  for  predicting  the  transverse  force  worked  well  for  both  the 
magnitude  and  frequency.  When  applied  to  the  data  from  these  experiments,  this 
model  showed  that  the  magnitude  of  the  transverse  force  (£/£,(„))  is  deterministic  and 
that  the  associated  phase  angle  <^(n)  is  random  in  nature. 

When  the  tower  was  allowed  to  respond  freely  in  the  transverse  direction,  interest- 
ing results  on  the  effect  of  transverse  motion  were  also  found.  The  most  pronounced 
effect  was  on  the  lift  coefficient,  frequency  of  vortex- shedding  and  vortex  span-wise 
coherence  length.  The  lift  coefficients  were  increased  by  a factor  of  at  least  2.5  due  to 
the  transverse  motion.  The  transverse  force  was  more  regular  than  when  transverse 
motion  did  not  exist  in  that  fewer  harmonics  (usually  2)  were  needed  in  the  proposed 
model  to  produce  the  best  fit  to  the  measured  data.  This  means  that  transverse 
motion  reduces  the  irregularities  in  the  transverse  force  by  synchronizing  the  shed- 
ding along  the  length  of  the  structure.  Spectra  for  the  thermistor  data  confirm  this 
explanation.  Regarding  the  frequency  of  vortex-shedding,  the  range  over  which  the 
dominant  frequencies  cluster  around  the  driving  frequency  (or  one  of  its  multiples) 
was  found  to  depend  on  (3,Vr,  and  the  amplitude  of  transverse  motion  rather  than 
KC  and  (3  as  for  the  cases  when  the  tower  is  not  allowed  to  have  transverse  motion. 
The  vortex  span-wise  coherence  length  was  found  to  be  much  larger  for  all  values  of 
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/3  tested  and  its  dependency  on  KC  and  /?  seems  to  not  be  influenced  by  the  trans- 
verse motion.  For  the  small  values  of  /3  the  energy  of  the  response  is  much  larger 
at  the  driving  frequency  than  at  the  system’s  natural  frequency.  As  /3  increased  the 
energy  at  the  natural  frequency  increased  until  it  reached  its  maximum  at  a value  of 
/?  = 12,415.  At  this  value  of  /?,  lock-in  occurred  over  a wide  range  of  reduced  ve- 
locity (22  < Vr  < 36)  when  the  dominant  frequency  of  vortex-shedding  matched  the 
response  frequency  and  a multiple  of  the  system  natural  frequency.  As  (3  increased 
further  the  amplitude  started  to  decrease. 

As  explained  before,  the  transverse  motion  increases  the  transverse  force  which  in 
turn  increases  the  transverse  motion.  This  process  continues  until  an  equilibrium  or 
ultimate  state  condition  is  reached  for  the  given  flow  conditions.  It  is  unfortunate  that 
the  apparatus  used  in  these  experiments  did  not  allow  the  largest  motion  that  could 
be  induced  due  to  vortex-shedding  for  the  range  of  KC  tested.  The  maximum  yrms 
measured  was  about  0.6 D,  with  a corresponding  1/1/3  of  about  0.95 D.  For  the  range 
of  KC  considered  the  results  show  that  as  KC  increases,  the  maximum  amplitude 
of  transverse  motion  increases.  Based  on  the  results  of  other  investigators,  it  is 
anticipated  that  these  amplitudes  will  go  through  a maximum  at  approximately  KC  = 
11. 

6.2  Conclusions 

This  study  started  with  an  extensive  literature  review  on  the  experimental  in- 
vestigations relevant  to  vortex-induced  transverse  loading  on  cylinders  in  water.  The 
review  revealed  a lack  of  information  and  understanding  of  the  effect  of  flow  nonunifor- 
mity and  three  dimensionality  on  the  vortex-shedding  process,  span-wise  coherence 
length  and  the  vortex-induced  loading.  Regardless  of  the  type  of  flow  there  were 
almost  no  quantitative  studies  on  the  effect  of  structure  transverse  motion  on  the 
quantities  mentioned  above.  The  objectives  of  this  study  were  thus  1)  to  gain  in- 
sight into  the  effects  of  flow  nonuniformity  and  transverse  motion  on  vortex-induced 


121 

transverse  loading  through  carefully  planned  laboratory  experiments,  2)  to  formulate 
a better  mathematical  model  for  predicting  transverse  loads  in  oscillatory  flows  and 
3)  to  produce  lift  coefficient  data  for  a meaningful  range  of  KC  and  Re.  A complete 
set  of  data  on  the  lift  coefficients  when  motion  exists  could  not  be  obtained  due  to  the 
lack  of  reliable  values  of  Cm  and  Cd  for  the  range  of  KC  imposed  by  the  transverse 
motion. 

The  study  increased  the  data  base  for  the  hydrodynamic  coefficients,  Cd  and 
Cm  for  ranges  of  KC  and  Re  where  very  little  has  been  reported  for  either  fixed 
or  compliant  structures.  The  results  obtained  in  this  study  should  be  useful  to  the 
design  of  articulated  mooring  towers,  risers,  TLP  tendons  and  other  components  of 
offshore  structures. 

During  the  course  of  this  study,  several  facts  emerged  upon  which  the  following 
conclusions  may  be  drawn: 

1.  Oscillating  an  articulated  cylinder  in  otherwise  still  water  appears  to  be  a rea- 
sonable technique  for  studying  the  effect  of  flow  nonuniformity  on  flow-induced 
loadings. 

2.  The  technique  used  for  measuring  frequency  of  vortex-shedding,  i.e.,  miniature 
thermistors,  was  successful.  The  data  obtained  provided  information  on  the 
vortex  span-wise  coherence  length  and  the  general  vortex-shedding  process  for 
both  transversely  constrained  and  unconstrained  cylinders. 

3.  The  proposed  model  for  predicting  transverse  forces  is  promising.  The  mag- 
nitude and  frequency  of  the  transverse  force  can  be  predicted  to  a reasonable 
degree  of  accuracy  using  this  model. 

4.  Even  though  more  data  are  needed  before  the  effect  of  flow  nonuniformity  and 
transverse  motion  can  be  quantified  for  a wide  range  of  KC , Re  and  VT,  progress 
on  this  problem  has  been  made  in  this  study.  The  results  of  this  work  cover 
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a range  of  KC  form  2.5  to  9.5  and  a range  of  Re  from  6.1  x 103  to  1.3  x 10s 
(ranges  where  very  little  data  exists  in  the  literature). 

5.  For  the  range  of  KC  and  Re  tested  the  in-line  forces  are  inertia  dominated.  CD 
is  more  dependent  on  flow  nonuniformity  than  Cm.  Both  coefficients  depend  on 
Re  and  KC. 

6.  The  vortex-shedding  process  and  the  transverse  forces  for  the  two  cases  tested 
(transversely  constrained  and  unconstrained)  are  remarkably  different.  For  the 
constrained  tower,  the  transverse  force  exhibits  multiple  frequencies.  Its  lift 
coefficients  are  deterministic  and  well  behaved,  while  the  associated  phase  angles 
appear  to  be  random.  The  magnitude  of  the  transverse  force  and  the  frequency 
of  vortex-shedding  depend  on  KC  and  ft.  The  largest  value  of  the  lift  coefficient 
(for  the  range  of  KC  tested)  occurs  at  the  in-line  driving  frequency. 

7.  For  the  transversely  constrained  tower  there  appears  to  be  a “critical”  frequency 
parameter  (Per)  where  the  transverse  force  peaks.  Thermistor  data  show  a 
maximum  vortex  span-wise  coherence  length  at  this  value.  In  this  study  the 
critical  value  of  /?  is  2, 160. 

8.  Transverse  motion  reduces  irregularities  in  the  transverse  force  and  provides 
a means  of  synchronizing  the  vortex-shedding  along  the  length  of  the  tower. 
Fewer  harmonics  in  the  proposed  model  are  needed  to  predict  the  transverse 
force  when  motion  is  allowed. 

9.  Transverse  motion  significantly  increases  the  value  of  the  lift  coefficients.  In 
this  study,  the  lift  coefficient  was  at  least  2.5  times  larger  than  that  for  the 
constrained  tower  under  the  same  environmental  conditions.  Lift  coefficients 
obtained  for  constrained  towers  underpredict  vortex-induced  response  by  at 
least  40%.  The  calculated  response  using  constrained  tower  lift  coefficients  is 
also  more  irregular  than  that  measured. 
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10.  Transverse  motion  influences  the  frequency  of  vortex-shedding  and  the  vortex 
span-wise  coherence  length.  The  frequency  of  vortex-shedding  and  the  range 
over  which  the  dominant  frequencies  cluster  about  a particular  frequency  de- 
pends onVr,  ft  and  the  amplitude  of  transverse  motion.  The  vortex  span- wise 
coherence  length  is  much  larger  when  transverse  motion  exists. 

11.  Lock-in  occurs  when  the  amplitude  of  the  transverse  response  at  the  system 
natural  frequency  is  large  enough  to  control  the  shedding  process.  In  this  case 
the  dominant  frequencies  of  vortex-shedding  are  close  to  the  natural  frequency 
(or  one  of  its  multiples). 

12.  The  amplitude  of  transverse  motion  depends  on  Vr  and  KC.  Maximum  trans- 
verse motion  due  to  vortex-shedding  occurs  when  KC  > 4 and  VT  > 6.5.  As 
KC  increases  the  amplitude  of  maximum  response  increases  and  the  value  of 
Vr  at  which  this  maximum  occurs  also  increases. 

13.  The  force  coefficients  obtained  from  oscillating  an  articulated  tower  in  still  water 
agree  well  with  the  results  of  other  investigators  for  waves  on  fixed  cylinders. 

6.3  Recommendations  for  Further  Work 

This  study  introduced  a new  method  for  measuring  frequency  of  vortex-shedding; 
proposed  an  improved  model  for  calculating  vortex-induced  transverse  forces;  and 
answered  some  of  the  questions  concerning  the  effects  of  flow  nonuniformity  and 
structure  transverse  motion  on  the  vortex-shedding  process,  the  transverse  loading 
and  the  vortex  span- wise  coherence  length.  The  problem;  however,  is  by  no  means 
solved  and  additional  work  is  needed. 

A major  uncertainty  is  in  the  mechanism  that  determines  the  phase  of  the  vortex- 
induced  transverse  force.  This  needs  to  be  investigated  further  and  a methodology 
developed  for  dealing  with  its  seemingly  random  behavior.  In  addition,  tests  similar 
to  those  conducted  here  need  to  be  performed  for  a wider  range  of  KC  and  Re.  More 
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specifically  the  range  of  KC  between  1.5  and  20  and  Re  between  104  and  2 x 106. 
The  influence  of  in-line  motion  on  transverse  forces  and  vortex  span-wise  coherence 
length  needs  to  be  included  for  the  same  range  of  KC  and  Re.  Drag  and  inertia 
coefficients  for  the  low  values  of  KC  and  Re  existing  for  transverse  motion  need 
to  be  measured.  The  effect  of  different  system  natural  frequencies  on  the  lock-in 
phenomenon  should  also  be  examined.  Finally,  the  effects  of  surface  roughness  on  all 
of  the  above  quantities  is  needed. 

Fluid-structure  interaction  problems  are  complex  and  difficult  to  solve.  In  spite 
of  the  research  conducted  and  progress  made  in  recent  years  in  vortex-induced  load- 
ing, many  questions  remain  unanswered.  Problems  associated  with  the  phase  of  the 
transverse  force  are  tied  to  flow  instabilities  which  are  very  sensitive  to  even  minor 
changes  in  the  flow  and  in  structure  surface  conditions.  The  phase  angle  may  thus 
have  to  be  treated  as  a random  phenomenon  for  some  time  to  come. 

It  is  hoped  that  the  results  of  this  study  will  provide  some  insight  into  the  complex 
mechanisms  involved  and  that  the  data  obtained  will  be  useful  to  design  engineers. 


APPENDIX  A 
ANALYSIS  TECHNIQUES 


A.l  Formulation  of  Dimensionless  Parameters 

When  an  articulated  tower  is  allowed  to  have  transverse  motion,  the  transverse 
force  and  response  are  interdependent.  To  study  the  effect  of  this  motion  and  other 
parameters  on  the  vortex- induced  transverse  loading  and  response,  the  equation  of 
motion  is  investigated.  The  general  form  of  this  equation  can  be  written  as 


Im  m + C,  6{t)  + fa  6{t)  |0(*)|  + k,  9{t)  = ML(t ) 


(A.1) 


where, 


m 

c. 

Pk 

Cd 

k, 

ML(t ) 


Im  = h + I'a, 

I0  = mass  moment  of  inertia  about  bottom  hinge, 

added  mass  moment  of  inertia  about  bottom  hinge, 
angle  of  transverse  oscillation, 
structural  damping  coefficient, 


hydrodynamic  damping  coefficient, 
r*‘  l 

o...  CJn  D r°  dr  = 

8 

drag  coefficient, 


= JQ  ^ Pw  CD  Dr3 dr  = ^ pwCD  D £4a , 


(A.2) 


restoring  moment  coefficient  due  to  buoyancy  and  weight, 
moment  due  to  lift  force  about  bottom  hinge, 

E M0  CLn  (A.3) 


n=l 


Mn  = 


components  of  the  moment  at  in-line  flow  frequencies  ncr  and 

1 


= j^\pu>DU2rdr=l-pxuDU2ea. 


(A.4) 
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If  the  nonlinear  drag  term  is  linearized  with  respect  to  time,  a solution  of  the  form 


N 

0(t)  = J2e»(n)  e,[*(n)-n<Tt], 

n=l 


(A.5) 


can  be  assumed.  Thus,  substituting  Eqn.  A.5  into  Eqn.  A.l  and  eliminating  the  time 
function,  result  in, 


0o(n)e‘*<n>  = 


n = 1,2, ...,  N.  (A.6) 


[fc,  - It(na)2]  - *(ncr)[Ae  + C3]  ’ 
where,  Ae  is  the  integral  over  the  length  £a  and  is  a function  of  Q0  and  <f>(n).  Further 
it  can  be  assumed,  without  loss  of  any  generality,  that  Ae  is  negligibly  small,  then  the 
magnitude  of  the  transverse  oscillation  becomes: 

M0/It 


^ {[u>2  - (ncr)2]2  + [2/3hu>(na)]2}2 

The  phase  angle  with  reference  to  the  in-line  flow  is 

*(■») = *<n> + taD''CA-(Z)^ 


(A.7) 


(A.8) 


where  u = 2irfn  = yj(ka/lt ) . 

If  It  = I0  + Ia  — (m,  + ma)£2/3,  then  substituting  this  expression  along  with 
Eqn.  A.4  into  Eqn.  A.7  gives 

. J A.  CLn  D U2 PJmJi  1 + m.lm,)P/ 3 


0o(n)  = . 

{[w2  - (tut)2]2  + [2 /3hLu(n<T)]2}i 
Also,  if  0o(n ) ~ y/£,  then  Eqn.  A.9  reduces  to 

\pwD(U/ Dfn)2  D2  £\  12n=i  CLn 


(A.9) 


y/l  = 


(4  7r2 f2  /3w2)  ma(l  + ma/mj){[w2  - (tut)2]2  + [2/?fcO>(n<r)]2}  J 
Since  U/Dfn  = Vr, 

y/D  = 3 pw£V2D2£]  ELiClu 

16  7T2  m.(l  + ma/ms)  l2  {[1  - (rur/a;2)2]2  + ^(tut/uj)]2}  J ' 

Further,  if  la  = d/  cos  6 is  approximated  by  £ then  Eqn.  A.ll  reduces  to 

3pw£V2D2T,n=iCLn 


(A- 10) 


(A.ll) 
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16  7r2m,(l  +ma/ma){[l  - (rur/u)2]2  + [2^(n<7/u)]2}5 


1 • 


(A.12) 
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Introducing  the  critical  damping  factor  £ in  water: 

C . 


£ = & = 


2It  u) 


(A.13) 


where  8 = 2tt£  is  the  logarithmic  decrement  of  free  oscillation  in  still  water.  Then, 
at  the  lock-in  condition  (where  na  = u ) Eqn.  A. 12  reduces  to 

3P„ev?D*  eL,  cLn 


U>ir8mt(l  + ma/m,)' 

The  stability  parameter  is  defined  as 

_ 2m,  8/1 

Then  with  the  dominant  lift  coefficient  Cl  ~ CTn>  Eqn.  A.  14  becomes 

3 V?CL 


y/D  ~ 


87tA",(1  + ma/ms) 


(A.14) 


(A.15) 


(A.16) 


A. 2 Runge-Kutta  Numerical  Solution 

In  quantifying  the  effect  of  transverse  motion  on  transverse  forces,  Runge-Kutta 
numerical  technique  was  used  to  solve  Eqn.  4.54.  For  this  technique,  the  dynamic 
equations  are  transformed  into  two  first  order  differential  equations: 


and 

~dt  = 


(A.17) 
(A. 18) 


where  Y = The  function  g = AfappUed  where  the  applied  moments  have 
been  described  previously  in  Eqns.  4.50  to  4.52.  The  function  / is  simply  f = Y. 
The  two  fundamental  dependent  variables  Y and  Y are  extrapolated  from  the  jth 
time  step  to  the  next  step  according  to 


'dY' 
dt  , 


8t,  and 


Yj+ 1 — Yj  + 


dY_ 

dt 


8t , 


(A- 19) 
(A. 20) 
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where 


dY 

dt 

XY_ 

dt 


— ( Ko  + 2 Ki  + 2K2  + K3 ) 1 and 
o 

^ (m0  + 2m!  + 2 m2  + m3) , 


where,  according  to  the  Runge-Kutta  theory, 

K0=  f(tj,  Yj,fy  , 

K\  = f (tj  + ^,  Vj  + A'o y, Vj  + m0 y)  > 

K2  = f [tj  + Y’ + mi  y)  ’ 

A3  = / (tj  + tit,  Vj  + K^St,  Yj  + m2£t^  , 
m0  = tf(tj,*j,X?)  , 

mi  = 9 (tj  + + ^of, ^j  + m0 y)  ' 

m2  = 5 ^tj  + Y'  Yj  + Ki  y'j  V}  + mi  Yy  , and 
m3  = g (tj  + 6t,  Vj  + /f2£t,  Yj  + m26t)  . 


(A.21) 

(A.22) 


APPENDIX  B 

INSTRUMENTATION  AND  CALIBRATION  DATA 


In  this  appendix,  the  instrumentations  used  in  these  experiments  are  presented 
along  with  pertinent  calibration  data. 
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^ |o.35" 

J 0.625 
1 0.4" 

ii  i — i.o^ — { !i 

Force  Transducer 


Wheatstone  Bridge 


Figure  B.l:  X-Y  force  transducer  and  wheatstone  bridge  circuit. 
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Table  B.l:  Specifications  of  strain  gage  amplifier/signal  conditioning  module. 


item 

gain  range 

gain  range  with  external  R 

gain  temp,  coefficient 

drift  (RTI) 

input  bias  current 

input  impedanceidifferential 

input  noise  (RTO)  at  gain  = 100 

10-200 
10-1000 
3 ppm/°C 
2 ± § pV/°C 
30  nA 
3000 

1 - 100  Hz  : \20pVRMS 

Figure  B.2:  Schematic  of  strain  gauge  amplifier/signal  conditioning  module. 
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output  voltage  (volt) 

Figure  B.3:  Calibration  curve  for  in-line  force  (for  test  runs  at  0.4  > fd  < 0.8  Hz.). 


output  voltage  voltage  (volt) 


Figure  B.4:  Calibration  curve  for  in-line  force  (for  test  runs  at  0.15  > fd  < 0.3  Hz.). 


transverse  moment  (N.m) 
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Boss  Detail  MS33514-8  (O-nng) 

3/4-16  UNF-3A  Thread 

0.375  in.  (9-52mm>  O.D.  3/32  Wall; 
304  Stainless  Steel 


t015  in.  (3.81mm) 
Reference 
_ 2.13  in.. 
(54.10mm) 


ZOO  in.(50mm) 
Null  Space 


1 .00  in.  to  300  in. 
(25.40mm  to  7620mmj 
Acuve  Stroke 


2.50  in. 
(63.50mm) 
_Dead 
Zone 


Hex,  0.375  in.  (9.52mm) 

1.75  in.  (44.45mm)  Across  Flats 


(15.7mm) 


Transducer  with  integral  cable  (pig tad) 


Figure  B.6:  Standard  rigid  linear  displacement  transducer. 
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Table  B.2:  Linear  displacement  transducer  specifications. 


Input  Voltage 

±12Vdc  to±15Vdc 

Current  Draw 

±15Vdc  @ 100mA  maximum,  25mA  minimum 
(current  draw  varies  with  magnet  position,  maximum 
draw  occurs  when  magnet  is  2 in.  (50  mm)  from 
the  flange  and  minimum  update  time  is  being  utilized). 

Null  Space 

Null  space  will  be  2 in.  if  the  stroke  length  is  stated  in  inches 
and  50  mm  if  the  stroke  length  is  stated  in  millimeters. 

Displacement 

Up  to  25  ft.  (7620  mm) 

Dead  Space 

2.5  in.  (63.5  mm) 

Head  Enclosure 

2.13  in.  (54.10  mm)  high  with  1.75  in.  (44.45  mm)  dia.; 
NEMA6 

Nonlinearity 

< 0.05%  of  full  stroke;  minimum  of  ±0.002  in.  £±0.05  mm) 

Repeatability 

< 0.001%  of  full  stroke  or  ±0.0001  in.  (±0.0025  mm), 
whichever  is  greater 

Hysteresis 

0.0008  in.  (0.02  mm)  maximum 

Frequency  Response 

Stroke  dependent,  200Hz  to  50Hz  is  typical  for  lengths 
of  12  in.  (304.8  mm)  to  100  in.  (2540  mm)  respectively. 
Wider  response  is  available. 

Temperature  Coefficient 

less  than  0.0001 1 inV°F  + [3ppm/°F  per  in.  of  full  stroke] 
(less  than  .005mm/°C  + [0.21ppm/°C  per  mm  of  full  stroke]) 

Operating  Temperature 
Transducer  Rod 
Head  Electronics 

-40  to  185°F  (40  to  85°Q 
40  to  150°F  (40  to  66°C) 

Operating  Pressure 

FM  certified  to  3000psi  continuous;  8000psi  static 

Output  Impedance 

10Q 
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displacement  (m) 

Figure  B.7:  Calibration  curve  of  in-line  linear  displacement  transducer. 


Figure  B.8:  Calibration  curve  of  transverse  linear  displacement  transducer. 


tn  ♦ 
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Table  B.3:  Thermistor  specifications. 


item 

FP07 

maximum  diameter  (D) 
length  (L) 

thermal  time  constant  (in  still  air) 
thermal  time  constant  (water  plunge) 
dissipation  constant  (in  still  air) 
dissipation  constant  (in  still  water) 
maximum  power  rating 
resistance  range  (ohms) 

2.15  mm 
12.5  mm 
0.1  sec 
7 msec 
0.05  mW/°C 
0.25  m\V/°C 
.006  watts 

1 K - 10  M 

c 

* Note  2 


NOTES: 

* 1.  Select  Resistors  to  establish 

a 2 Volt  D.C.  Level  at  Point  "A". 

* 2.  Select  C to  filter  high  frequency 

noise  if  required. 

* 3. 1C  is  LM358 


SIGNAL 

OUT 


Figure  B.9:  Schematic  of  thermistor  signal  processing  circuit. 
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spectral  density 
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Figure  B.ll:  Thermistor  signal  power  spectra  at  various  Reynolds  numbers. 
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Ri 

AAA/V 


R2 

AMAt 


V: 


Cirt: 


C2  rz 


v0 


Rt  = R2  = 10  K 
^1  = ^2  = ^ P F 


cr=  1.2  ESC 

»»  To*psr*turo=  27  Ca*e=  1 


Figure  B.12:  Schematic  diagram  and  frequency  response  for  second  order,  passive, 
Butterworth  lowpass  filter  circuit  (designed  for  the  feed  back  control  signal). 
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unfiltered  transverse  force  time  series 


filtered  transverse  force  time  series 


01  23456789  10 

Q frequency  [Hz] 


power  spectrum  of  filtered  transverse  force  signal 


£ 0.01 


o 

b 


S 6 

frequency  [Hz] 


Figure  B.13:  Samples  of  transverse  force  signals  and  their  corresponding  power  spectra 
(a)  before  filtering  (b)  after  filtering. 


Table  B.4:  Linear  drive  motor  specifications. 


motor  width  XY 

units 

45 

force  constant  ( Kp ) 

lbs/amp 

5.2 

continuous  force  ( Fc ) 

lbs 

26 

force  at  50%  D.C.  ( Fm ) 

lbs 

37 

force  at  10%  D.C.  ( Fp ) 

lbs 

82 

continuous  current  (7C) 

amps 

5.1 

current  at  50%  D.C.  (/m) 

amps 

7.1 

current  at  10%  D.C.  ( Iv ) 

amps 

16 

volts  for  Fc  at  10  inch /sec 

VDC 

28 

volts  for  Fc  at  20  inch /sec 

VDC 

34 

volts  for  Fc  at  30  inch /sec 

VDC 

40 

volts  for  7^  at  10  inch /sec 

VDC 

76 

volts  for  Fp  at  20  inch /sec 

VDC 

81 

volts  for  Fv  at  30  inch /sec 

VDC 

87 

back  EMF  Ke 

V/inch/sec 

.59 

resistance  at  20  C Rc 

ohms 

2.7 

resistance  at  180  C Rj 

ohms 

4.4 

inductance  L 

mhenry 

7.1 

time  constant  T at  20  C 

msec 

2.6 

magnetic  attraction 

lbs 

270 

APPENDIX  C 

FLOW  CHARTS  OF  THE  VARIOUS  COMPUTER  PROGRAMS 


In  this  appendix,  the  flow  charts  of  the  various  computer  programs  written  to 
analyze  the  experimental  data  are  presented. 
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Figure  C.l:  Computer  code  flow  cnari  101 
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input  data 


Figure  C.2:  Computer  code  flow  chart  for  CDCMN 
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input  data 


calculate  geometr 
matr i ces 
A and  8 


decompose  matrix 
[A] 


subroutine  SVDCMP 
decompses 
[A] 


calculate  f(t) 
from  measured 
force 


compute  Fourier 
coefficients  of 
f(t) 

subrout i ne 
FFT 

calculate  lift 
coef f i ci ents 
CL(n) 

subroutine 
SVBKSB  to  solve 
A-X  - B 

calculate 
associated  phase 
angles  <t>(n) 


print  out 
resul ts 


Figure  C.3:  Computer  code  flow  chart  for  CLEF. 
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decompose  matrix 
[A] 

subroutine  SVDCMP 
decompses 
[A] 

calculate  f(t) 
from  measured 
force 


calculate 
y and  y 

f rom  measured  y 


compute  Fourier 
coefficients  of 
f CO 


calculate  I i ft 
coef  f i c i ents 
CL(n) 


subrout i ne 
FFT 


subrout i ne 
SVBKS8  to  solve 
A-X  - B 


calculate 
associated  phase 
angles  <J>C n ) 


print  out 
resul ts 


Figure  C.4:  Computer  code  flow  chart  for  CLEUF. 
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